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Abstract 

This  thesis  is  motivated  by  a  series  of  biological  experiments  that  display  the  harbor 
seal’s  extraordinary  ability  to  track  the  wake  of  an  object  several  seconds  after  it  has 
swum  by.  They  do  so  despite  having  auditory  and  visual  cues  blocked,  pointing  to 
use  of  their  whiskers  as  sensors  of  minute  water  movements.  In  this  work,  I  elucidate 
the  basic  fluid  mechanisms  that  seals  may  employ  to  accomplish  this  detection.  Key 
are  the  unique  flow-induced  vibration  properties  resulting  from  the  geometry  of  the 
harbor  seal  whisker,  which  is  undulatory  and  elliptical  in  cross-section. 

First,  the  vortex-induced  vibration  (VIV)  characteristics  of  the  whisker  geometry 
are  tested.  Direct  force  measurements  and  flow  visualizations  on  a  rigid  whisker  model 
undergoing  a  range  of  1-D  imposed  oscillations  show  that  the  geometry  passively 
reduces  VIV  (factor  of  >  10),  despite  contributions  from  effective  added  mass  and 
damping.  Next,  a  biomimetic  whisker  sensor  is  designed  and  fabricated.  The  rigid 
whisker  model  is  mounted  on  a  four-armed  flexure,  allowing  it  to  freely  vibrate  in  both 
in-line  and  crossflow  directions.  Strain  gauges  on  the  flexure  measure  deflections  at 
the  base.  Finally,  this  device  is  tested  in  a  simplified  version  of  the  fish  wake  -  seal 
whisker  interaction  scenario.  The  whisker  is  towed  behind  an  upstream  cylinder  with 
larger  diameter.  Whereas  in  open  water  the  whisker  exhibits  very  low  vibration  when 
its  long  axis  is  aligned  with  the  incoming  flow,  once  it  enters  the  wake  it  oscillates 
with  large  amplitude  and  its  frequency  coincides  with  the  Strouhal  frequency  of  the 
upstream  cylinder.  This  makes  the  detection  of  an  upstream  wake  as  well  as  an 
estimation  of  the  size  of  the  wake-generating  body  possible.  A  slaloming  motion 
among  the  wake  vortices  causes  the  whisker  to  oscillate  in  this  manner.  The  same 
mechanism  has  been  previously  observed  in  energy-extracting  foils  and  trout  actively 
swimming  behind  bluff  cylinders  in  a  stream. 

Thesis  Supervisor:  Michael  S.  Triantafyllou 

Title:  Professor  of  Mechanical  and  Ocean  Engineering 
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blindfolded.  (B)  Same  setup,  but  the  second  seal  now  employs  a  more 


distinctively  undulatory  pattern  to  the  search.  Figure  from  [118].  .  .  54 

1-21  Different  shapes  that  a  seal  was  trained  to  distinguish  based  solely  on 

the  wake  left  while  they  were  held  fixed  and  towed.  Figure  from  [140].  55 

1-22  Ocean  structures,  such  as  offshore  platforms  (figure  from  [66])  and 

cables  (figure  from  [35]),  are  subject  to  VIV .  57 

1-23  A  bluff  body  in  a  uniform  flow  can  be  simply  represented  as  a  1-degree- 
of- freedom  linear  oscillator.  The  free  stream  approaches  with  velocity 

U  causing  vortex  shedding  at  a  frequency  of  fs .  59 

1-24  Examples  of  VIV  suppression  devices:  (a)  fairing,  (b)  splitter  plate, 

(c)  flag,  and  (d)  a  cylinder  with  helical  strakes.  Figures  from  [35].  .  .  66 
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1-25  The  drag  coefficient  of  various  VIV  suppression  devices  is  higher  than 
that  of  a  smooth  cylinder  at  certain  Reynolds  numbers.  This  is  espe¬ 
cially  true  for  a  cylinder  with  helical  strakes.  Figure  from  [35] .  67 

1-26  (a)  A  cylinder  with  bumpy  protrusions  placed  at  45°  angular  separation 
and  with  pitch  of  7  diameters.  Figure  from  [100].  (b)  The  wake  behind 
this  model  ranges  from  wide  to  narrow,  depending  on  whether  the 
model  has  a  peak  or  trough  in  the  disturbance  there.  Figure  from  [99].  67 

1-27  Maximum  transverse  vibration  amplitude  of  different  geometries.  As 
the  mass-damping  (m*£)  decreases,  all  geometries  vibrate  with  higher 
amplitude,  solid  circle :  smooth  circular  cylinder;  plus  sign :  sinuous 
cylinder;  others :  bumpy  cylinders  with  various  bump  sizes).  Figure 
from  [99] .  68 

1-28  Schematic  of  expected  response  of  the  whiskers  as  the  seal  is  swimming: 

the  whiskers  undergo  significant  VIV,  making  for  a  “noisy”  signal.  .  .  69 

1- 29  Schematic  of  the  hypothesis  behind  this  thesis.  In  open  flow,  the 

whiskers  produce  low  VIV.  This  enables  the  seal  to  more  easily  de¬ 
tect  unsteady  flow  features,  which  induce  large  response .  70 

2- 1  Close-up  view  of  the  whiskers  on  a  harbor  seal’s  muzzle.  The  “beaded” 

pattern  indicates  an  undulatory  geometry.  Photo  taken  courtesy  of  the 
New  England  Aquarium .  74 

2-2  (Top)  Side  view  and  (Bottom)  front  view  of  a  fabricated  replica  of 
the  whisker  test  section,  capturing  the  asymmetric  undulatory  geome¬ 
try.  Stereolithography  was  used  to  create  this  plastic  hollow  structure. 
Dimensions  were  taken  directly  from  [58]  and  enlarged  by  a  factor  of  50.  75 

2-3  Seals’s  whiskers  may  face  different  orientations.  Photo  taken  at  the 

New  England  Aquarium .  76 
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2-4  Schematic  of  the  forced  vibration  setup.  The  model  hangs  vertically 
in  the  water  tank.  It  is  directly  attached  to  a  force  sensor  and  a  linear 
stage  that  drives  a  sinusoidal  motion  (an  example  trajectory  is  depicted 
in  red)  while  the  carriage  moves  the  setup  down  the  length  of  the  tank 
at  velocity  U .  77 

2-5  Orientation  of  a  cross  section  of  the  model.  The  solid  line  cross  section 
represents  the  a  =  0°  case  and  the  dotted  line  cross  section  represents 
the  a  >  0°  cases.  The  projected  diameter  (d*)  is  the  dimension  used  for 
all  cases  and  is  the  length  perpendicular  to  the  flow,  which  approaches 
with  velocity  U .  79 

2-6  Schematic  of  the  PIV  setup.  The  model  again  hangs  vertically  in  the 
water  tank.  A  laser  sheet  illuminates  a  thin  cross-section  of  the  model 
and  the  fluid  in  that  plane.  A  camera  mounted  below  the  tank  captures 
the  motion  of  neutrally  buoyant  particles  in  the  water .  82 

2-7  Cd  is  plotted  as  a  function  of  amplitude  ratio  and  reduced  frequency  for 
the  circular  cylinder,  elliptical  cylinder,  and  whisker  (from  left  to  right). 

Re  =  5,300.  Asterisks  indicate  where  experiments  were  performed. 

The  contour  levels  are  in  increments  of  0.1  up  to  1,  then  in  increments 
of  0.25  to  1.5,  and  in  increments  of  0.5  beyond  that .  83 

2-8  Clv  is  plotted  as  a  function  of  amplitude  ratio  and  reduced  frequency. 

The  region  of  positive  values  indicates  where  VI V  occurs.  The  size  of 
this  region  decreases  from  the  circular  cylinder  to  the  elliptical  cylin¬ 
der  to  the  whisker  (from  left  to  right).  The  zero  contour  is  depicted 
in  bold.  Re  =  5,300.  Asterisks  indicate  where  experiments  were  per¬ 
formed.  The  contour  levels  are  in  increments  of  0.5  down  to  -0.5,  then 
in  increments  of  2  from  -2  on,  and  in  increments  of  -6  beyond  that.  .  84 
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2-9  Variation  seen  in  the  whisker  Ci)V  plot.  The  vertical  bars  indicate  the 
probability  that  that  point  on  the  plot  would  show  a  positive  C]  v  value. 

The  solid  black  dots  indicate  the  average  (over  10  trials)  Ci}V  value  at 
each  cycle  of  oscillation.  The  scatter  plots  (red)  show  the  spread  of  the 
10  individual  values  at  each  cycle.  The  standard  deviation  decreases 
over  time,  indicating  that  lock-in  is  reached  after  some  initial  start-up 
phase .  86 

2-10  Cm  contours  plotted  as  a  function  of  amplitude  ratio  and  reduced  fre¬ 
quency  for  the  circular  cylinder,  elliptical  cylinder,  and  whisker  (from 
left  to  right).  The  zero  Clv  contour  is  copied  over  and  shown  in  black. 

The  contour  levels  are  in  increments  of  1  and  the  3.5  contour  is  included 
as  well.  Asterisks  indicate  where  experiments  were  performed .  87 

2-11  Vorticity  contours  of  a  2-D  slice  of  the  wake  behind  the  elliptical  cylin¬ 
der  (left)  and  whisker  (right)  models  undergoing  steady  tow.  A  stan¬ 
dard  von  Karman  vortex  street  appears  behind  the  ellipse,  but  the 
whisker  decorrelates  the  vortex  shedding,  causing  the  vortex  street  to 
appear  much  farther  downstream,  where  the  forcing  on  the  body  is 
reduced .  88 

2-12  Vorticity  contours  of  a  2-D  slice  of  the  wake  behind  the  elliptical  cylin¬ 
der  (left)  and  whisker  (right)  models  undergoing  the  forced  vibration 
run  corresponding  to  their  maximum  VI V  amplitude.  A  standard  von 
Karman  vortex  street  appears  behind  both  models .  89 

2-13  Drag  coefficient  during  forced  vibrations  for  (left  column)  ellipse  and 
(right  column)  whisker  at  (a,b)  15;  (c,d)  30;  (e,f)  45;  (g,h)  60;  (i,j) 

75;  (k,l)  90°.  Contour  levels  for  15  and  30  °  are  in  increments  of  0.1. 
Otherwise  they  are  in  increments  of  0.5 .  92 

2-14  C[j,  contour  for  (left  column)  ellipse  and  (right  column)  whisker  at 
(a,b)  15;  (c,d)  30;  (e,f)  45;  (g,h)  60;  (i,j)  75;  (k,l)  90°.  Contour  levels 
are  in  increments  of  1 .  94 
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2-15  The  peak  amplitude  of  vibration  for  the  elliptical  cylinder  and  whisker 
increases  as  the  angle  of  attack  increases.  Re  =  5,300  —  7,030  from 
a  =  0°  —  45°  and  Re  =  2,530  -  3,210  from  a  =  60°  -  90° .  95 

2-16  Effective  added  mass  coefficient  during  forced  vibrations  for  (left  col¬ 
umn)  ellipse  and  (right  column)  whisker  at  (a,b)  15;  (c,d)  30;  (e,f)  45; 

(g,h)  60;  (i,j)  75;  (k,l)  90°.  Contour  levels  are  in  increments  of  0.5.  .  .  97 

2-17  PIV  of  the  wake  behind  an  elliptical  cylinder  (left)  and  whisker  (right) 

being  rigidly  towed  at  different  angles  of  attack.  Re  =  5,300 .  98 

2- 18  The  wake  shed  behind  the  freely  vibrating  whisker  for  the  (left)  45  and 

(right)  90°  cases.  (Experiments  conducted  by  B.  Seyed-Aghazadeh  and 
Y.  Modarres-Sadeghi) .  99 

3- 1  The  whisker  (yellow)  mounted  to  a  pressure  housing,  which  contains 

four  Bend  Sensors  (dark  orange) .  106 

3-2  (a)  A  square  plate  (in  black)  at  the  end  of  a  rod  connected  to  the 

whisker  pushes  on  the  four  Bend  Sensors  as  the  whisker  deflects,  (b) 

A  silicone  diaphragm  (in  orange)  between  the  whisker  and  a  retaining 
ring  allows  flexible  motion .  106 

3-3  The  electronics  that  are  stored  inside  a  separate  housing:  circuit  and 

data  logger.  The  battery  pack  is  not  pictured .  107 

3-4  Circuit  diagram.  The  Bend  Sensor  signal  lines  are  filtered  and  amplified.  108 

3-5  Setup  for  bench-top  calibration.  Deflection  of  the  whisker  was  con¬ 
ducted  in  increments,  using  the  ruler  on  the  left  as  a  reference .  108 

3-6  Voltage  output  versus  whisker  tip  displacement  for  each  of  the  four 

Bend  Sensors .  109 
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3-7  Towing  test  setup.  (Left)  The  sensor  is  mounted  to  the  carriage,  and 
a  data  logger  records  voltage  output.  Note:  The  wake  generator  (cir¬ 
cular  cylinder)  pictured  was  not  used  for  the  experiments  shown  in 
this  Chapter.  (Right)  Schematic  of  the  inside  of  the  sensor  housing 
(top  view).  The  Bend  Sensors  (purple)  are  numbered  with  this  con¬ 
vention  throughout.  Note:  For  a  =  0°,  Bend  Sensors  1  and  2  are  in 
the  crossflow  direction,  but  as  a  is  increased  to  90°,  Bend  Sensors  3 
and  4  become  fully  in  the  crossflow  direction .  Ill 

3-8  Time  traces  of  the  Bend  Sensor  outputs.  Black  vertical  lines  indicate 
start  and  end  of  the  carriage  motion.  (Top  row)  a  =  0°:  minimal 
change  occurs.  (Middle  row)  a  =  45°:  ramps  and  oscillations  in  voltage 
are  seen.  (Bottom  row):  a  =  90°:  larger  oscillations  occur  in  the 
crossflow  (Sensors  3,4)  and  larger  ramps  occur  inline  (Sensors  1,2). 
Sensors  1  and  4  are  shifted  by  a  constant  value  to  consolidate  plots.  .  112 

3-10  Field  test  setup.  The  sensor  is  mounted  to  the  CTD  rosette  in  a  manner 
that  allows  the  whisker  to  encounter  unobstructed  flow.  The  view  is 
from  the  port  side  of  the  ship.  (Photo  courtesy  of  Paul  Suprenand)  .  115 

3-11  (a)  Velocity  measurements  at  20m  depth  from  the  ADCP  and  (b)  volt¬ 
age  output  from  the  whisker-inspired  sensor .  116 

3-12  Orientation  of  the  ship  and  sensor.  The  ship  had  a  heading  of  -120°  to 
-140°  during  the  10mm  deployment.  The  arrow  indicates  the  general 
flow  direction  shown  by  the  ADCP.  The  numbers  indicate  the  orienta¬ 
tion  of  the  corresponding  Bend  Sensors .  116 

3-13  Isometric  view  of  the  wake  sensor,  which  consists  of  a  plastic  whisker 
model  mounted  on  a  four- armed  flexing  plate.  Strain  gauges  on  these 
arms  are  mounted  on  these  arms  for  measurement  of  the  whisker’s 
deflections .  119 
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3-14  (a)  A  Solidworks  simulation  of  the  strain  induced  on  the  flexure  as  the 
whisker  deflects,  (b)  Close-up  of  the  flexure  device.  A  spring  steel 
plate  made  using  the  water  jet  is  attached  to  a  bushing  (clear)  and  pin 
(gray)  at  the  end  of  each  arm.  This  design  reduces  coupling  between 
the  two  axes  by  allowing  rotation,  rather  than  deflection  in  the  off-axis.  119 

3-15  The  data-collection  electronics.  The  strain  gauges  are  conFigured  in 
a  Wheatstone  bridge,  using  additional  resistors.  A  PhidgetBridge  is 
used  to  filter  and  transmit  the  voltages  to  a  computer  through  USB.  .  120 

3-16  (Left)  A  linear  fit  was  placed  through  data  collected  for  different  tip 
deflections  of  the  whisker  in  x.  (Right)  A  pluck  test  in  air  was  used 
to  find  the  natural  frequency  of  the  whisker  in  y.  The  values  here 
(VynYxn)  are  centered  around  zero  by  subtracting  the  mean  of  the  signal.  121 

3-17  Picture  of  sensor  assembly  in  the  (left)  normal  configuration  and  (right) 
with  the  addition  of  set  screws  (circled  in  red)  to  restrict  bushing  rota¬ 
tion.  A  fourth  set  screw  through  the  remaining  bushing  is  not  visible 
from  this  view .  122 

3-18  Power  density  of  signals  from  y-direction  pluck  tests  in  water  for  the 
(left)  normal  configuration  and  (right)  configuration  with  no  bushing 
rotation .  122 

3-19  Vibrational  amplitude  in  both  axes  when  (left)  no  set  screws  were  used 

and  (right)  set  screws  were  used  to  restrict  bushing  rotation.  Ur  =  17.  123 


21 


4-1  (Left)  Top  view;  and  (Right)  side  view  of  the  experimental  setup.  A 
hydrodynamic  trail  was  generated  in  front  of  the  whisker  through  the 
use  of  a  cylinder,  which  shed  a  Karman  vortex  street.  Three  different 
downstream  distances  (dx)  were  tested  for  each  of  the  three  models, 
ranging  from  the  near  wake  to  the  far  wake.  Three  different  cylinder 
sizes  (dcyi)  were  also  tested.  Simultaneous  measurement  of  the  oncom¬ 
ing  wake  frequency  and  whisker  response  was  done  through  use  of  a 
pressure  sensor  placed  near  the  edge  of  the  wake.  A  mariotte  bottle 
supplied  dye  at  a  constant  flow  rate.  A  camera  underneath  the  model 
captured  a  bottom-up  view  of  the  wake-whisker  interaction . 129 

4-2  Two  variations  in  the  harbor  seal  whisker  geometry  are  captured  through 
the  depicted  parameters,  which  are  shown  here  at  the  scale  of  the  real 
seal.  The  inclination  angles  are  17.6°  and  15.27°  for  the  top  and  bot¬ 
tom  cross-sections  shown,  respectively.  The  geometry  from  Hanke  et 
al  (2010)  (“whisker  #1”)  has  a  A / Dm  ~  2,  whereas  the  new  geometry 
(“whisker  #2”)  has  a  X/Dm  ~  5.  The  models  used  in  this  chapter  are 
scaled  up  from  the  dimensions  shown  here  30  and  20  times,  respectively.  130 

4-3  The  model  whisker  (“whisker  #1”  pictured  here)  was  mounted  on  a 
set  of  flexing  plates  that  allow  free  vibration  in  in-line  and  crossflow 
directions.  Bushings  are  used  to  reduce  coupling  between  the  axes.  .  .  131 

4-4  VIV  response  of  each  model  positioned  in  the  streamlined  direction 
( a  =  0°).  Bars  depict  data  spread  over  three  repeat  runs.  (Left) 

Both  cross-flow  and  in-line  amplitudes  of  the  whisker  and  elliptical 
cylinder  are  significantly  smaller  than  those  of  the  circular  cylinder. 

The  frequencies  at  which  the  whisker  and  ellipse  vibrate  remain  close  to 
their  natural  frequency.  The  diagonal  dotted  line  indicates  a  frequency 
that  matches  the  Strouhal  relation.  (Right)  Both  whiskers  and  the 
elliptical  cylinder  vibrate  with  A*  <  0.2.  Whisker  #2  has  the  lowest 
VIV  amplitude .  135 
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4-5  Dye  visualizations  of  the  wake  behind  a  circular  cylinder,  elliptical 
cylinder,  and  whisker  (from  left  to  right)  at  a  =  0°.  The  first  two 
objects  shed  regular  2S  patterns,  whereas  the  whisker  #1  model  sheds 
an  incoherent  wake.  Images  were  taken  during  runs  corresponding  to 
each  model’s  peak  amplitude:  cylinder  at  U*  =  10,  ellipse  at  U*  =  28, 
and  whisker  #1  at  U*  —  28 .  136 

4-6  Dye  visualization  of  the  wake  behind  the  whisker  #2  model.  Similar 
to  the  whisker  #  1  model,  it  sheds  an  incoherent  wake.  The  image  was 
taken  during  the  run  corresponding  to  the  model’s  peak  amplitude, 
which  is  at  U*  =  32 .  136 

4-7  For  a  whisker  within  the  wake  of  an  upstream  cylinder  ( dcyi  =  2.5 dw) 
at  a  distance  dx  =  8dw\  (a)  The  whisker  vibrates  with  far  higher  ampli¬ 
tude  in  the  wake  than  it  does  on  its  own.  (b)  The  peak  frequency  of  the 
whisker’s  vibration  synchronizes  with  the  oncoming  wake  frequency. 

The  diamonds  and  “x”s  are  the  average  over  three  repeat  runs.  The 
bars  indicate  the  data  spread  of  the  measured  wake  frequency  over  the 
same  three  runs .  137 

Schematic  of  the  camera  perspective  of  the  dye  runs.  The  bottom-up 
view  shows  vortices  being  shed  by  the  wake  generator  on  the  right  and 
impinging  on  a  cross-section  of  the  whisker  on  the  left .  138 

Sequence  of  dye  visualizations  (left)  and  schematics  (right)  showing  the 
whisker’s  free  vibration  path  through  the  vortex  wake.  The  frequency 
of  the  whisker’s  motion  aligns  with  the  vortex  wake  frequency,  as  the 
whisker  slaloms  between  the  vortices.  (U*  =  16,  Re  =  3, 180) .  139 

Sample  result  of  the  cross-correlation  analysis.  (Top)  The  normalized 
signals  of  strain  and  pressure  measured  on  the  whisker  and  in  the 
wake,  respectively.  (Bottom)  The  cross-correlation  coefficient  at  each 
oscillation  cycle.  (dx  =  17 dw,  dcyi  =  4 dw,  U*  =  27) .  140 
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4-11  Cross-correlation  metrics  of  the  whisker  vibration  with  the  wake  fre¬ 
quency  for  three  different  size  cylinders  at  dx  =  15 dcyi.  (Left)  All  cases 
see  a  high  maximum  cross-correlation  coefficient.  (Right)  The  percent¬ 
age  of  the  run  where  the  whisker  is  “synchronized”  generally  increases 
with  cylinder  diameter .  141 

4-12  Average  of  percent  synchronization  across  all  speeds  tested  ( dcyi  = 

4 dw).  The  trends  are  displayed  versus  separation  distance  both  in 
the  downstream  (dx)  and  lateral  (dy)  directions.  Experiments  with  a 
lateral  offset  were  conducted  only  on  one  side  of  the  cylinder  but  the 
results  are  assumed  to  be  symmetric  and  are  therefore  mirrored  here.  142 

4-13  The  frequencies  of  response  nondimensionalized  by  the  theoretical  fre¬ 
quency.  At  the  close  distance  (8 dw),  all  four  models  vibrate  at  the 
theoretical  wake  frequency.  As  the  distance  increases  to  17 dw,  the  el¬ 
liptical  cylinder  and  whisker  continue  to  oscillate  at  the  wake  frequency, 
but  the  cylinder  deviates  from  it .  143 

4-14  Vibration  amplitudes  of  the  circular  cylinder,  elliptical  cylinder,  and 
whisker  (from  left  to  right,  respectively).  The  latter  two  vibrate  more 
when  they  are  in  the  wake  of  an  upstream  cylinder  (colored  lines)  than 
they  do  in  open  water.  The  circular  cylinder  vibrates  with  a  similar 
amplitude  for  all  cases.  The  amplitudes  of  response  decrease  as  the 
distance  from  the  upstream  cylinder  increases.  The  whisker  retains  an 
amplitude  greater  than  its  VIV  amplitude  even  at  the  largest  distance.  144 


4- 15  The  detectability  of  the  whisker  is  higher  than  the  other  models.  As 

the  separation  distance  increases,  the  wake  becomes  more  difficult  to 
detect,  but  the  whisker  remains  the  highest .  145 

5- 1  Pictures  of  the  whisker-inspired  flow  velocity  and  direction  sensor.  .  .  150 

5-2  Picture  of  the  re-designed  whisker  sensor,  which  was  made  for  system¬ 
atic  testing  in  a  vortex  wake .  150 
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5-3  Schematic  of  the  slaloming  trajectory  (dotted  line)  that  the  whisker 
takes  as  it  encounters  a  vortex  wake . 

5-4  (a)  Front  and  (b)  side  views  of  the  arrangement  of  whiskers  as  a  harbor 

seal  swims  forward.  Screenshots  taken  from  supplementary  video  in  [95].  157 

A-l  Digital  holographic  images  of  whiskers  from  a  harbor  seal  at  the  New 
England  Aquarium.  All  are  in  the  90°  configuration.  The  camera  field 
of  view  and  therefore  the  dimensions  of  each  DH  image  box  are  1.23  x 
1.23cm .  160 

A- 2  Digital  holographic  images  of  whiskers  from  a  harbor  seal  at  the  New 
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Chapter  1 


Introduction 


In  the  dark  and  murky  waters,  a  harbor  seal  swims  all  around,  searching  through  the 
turbid  waters  for  some  fish  to  eat  for  dinner.  Suddenly,  the  whiskers  on  her  right 
cheek  start  twitching.  And  she’s  off. 

Prey  had  passed  by  several  seconds  earlier,  and  biologists  at  the  University  of 
Rostock’s  Marine  Science  Center  [67]  were  clued  into  the  harbor  seals’  uncanny  ability 
to  chase  it.  These  scientists  placed  a  blindfold  and  earmuffs  on  a  seal  (Figure  1-1)  and 
found  that  the  animal  was  nonetheless  able  to  detect  and  track  moving  objects.  They 
found  evidence  that  the  whiskers  were  being  used  to  “feel”  what  was  happening  in  the 
surrounding  water,  based  on  how  they  vibrate.  It  was  already  known  that  seals  (and 
many  other  animals)  can  use  their  whiskers  in  a  tactile  manner  to  determine  the  size 
and  shape  of  objects.  But  discovery  of  their  potential  ability  to  detect  hydrodynamic 
information  opens  exciting  new  doors  on  our  understanding  of  the  biological  world. 

On  the  other  side  of  the  Atlantic,  research  was  being  conducted  in  my  lab,  the 
MIT  Towing  Tank  Laboratory,  on  the  vibrations  of  large  offshore  structures.  For 
the  last  couple  of  decades,  our  research  group  has  studied  the  physics  of  flow  across 
a  circular  cylinder.  This  situation  can  lead  to  vibrations  of  the  cylinder,  a  major 
concern  for  the  offshore  industry.  Simultaneous  to  these  studies,  our  group  has  been 
pioneering  the  field  of  marine  biomimetics.  This  is  where  we  use  features  observed  in 
nature  to  inform  the  design  of  new  marine  technologies. 
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Figure  1-1:  A  harbor  seal  wearing  a  blindfold  and  earmuffs,  ready  for  an  experiment. 
Photo  from  Marine  Science  Center,  Rostock  University. 

When  the  news  of  the  biologists’  work  reached  us,  we  jumped  at  the  opportunity 
to  marry  these  interests  in  vibrations  and  bio-inspired  design.  The  goal  on  our  minds, 
which  has  resulted  in  this  thesis,  was  to  uncover  how  the  seals  could  use  vibrations  as 
sensory  information.  Our  approach  is  to  study  the  fluid  mechanics  of  the  vibrations 
induced  by  the  whiskers’  interactions  with  the  surrounding  flow. 

1.0.1  How  might  seals  detect  a  fish  trail? 

Fish  leave  “footprints”  of  a  consistent  pattern  in  the  water.  Each  flick  of  their  fin 
causes  the  fluid  to  swirl  in  a  vortex,  which  lasts  for  many  seconds  after  the  fish  is 
gone.  Seals  could  reach  out  their  whiskers  to  feel  for  these  tracks  (Figure  1-2). 

It’s  surprising  that  this  would  actually  work  though.  As  our  research  group  knows 
well,  cylindrical  objects,  such  as  whiskers  or  offshore  riser  legs,  will  naturally  vibrate 
as  calm  water  flows  over  them.  So  how  is  it  possible  for  the  seal  to  notice  that  their 
whiskers  have  hit  the  fish  tracks,  if  their  whiskers  are  already  vibrating  a  lot  on  their 
own? 

The  biologists  in  Rostock  mounted  a  camera  to  a  harbor  seal’s  head  in  order  to 
get  a  close-up  view  of  the  whiskers.  What  they  noticed  was  that  the  whiskers  did 
not  appear  to  vibrate  as  the  seal  was  swimming  forward  normally.  This  observation 
suggests  that  the  whisker  is  tuned  for  heightened  sensing,  as  it  will  vibrate  minimally 
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Figure  1-2:  Schematic  of  a  seal  searching  for  a  fish  by  using  its  whiskers  to  detect  the 
hydrodynamic  trail. 

until  an  upstream  disturbance  is  encountered. 

In  this  chapter,  the  thesis  is  introduced  in  terms  of  two  topics:  first,  the  biological 
studies  that  inspire  the  work  and  second,  the  fundamentals  of  flow-induced  vibrations. 
The  remainder  of  the  thesis  takes  place  in  the  new  space  where  these  two  topics 
combine. 


1.1  Flow  Sensing 

Pinnipeds  often  feed  in  conditions  with  reduced  visibility  at  night  [116],  in  dark  waters 
[69] ,  even  diving  to  depths  >  300m  while  foraging  [77] .  Recent  behavioral  experiments 
suggest  that  the  reason  for  their  success  in  doing  so  is  the  use  of  their  highly  sensitive 
whiskers  to  detect  water  movements.  Their  whiskers  are  heavily  innervated  with  fibers 
that  detect  small  deflections  of  the  shaft.  Just  as  rats  and  cats  use  their  whiskers  as 
a  tactile  sensor  for  active  touch  discrimination  [31]  and  possibly  obstacle  avoidance 
[98] ,  seals  do  so  as  well  [23] .  The  study  of  their  potential  use  as  hydrodynamic  sensors, 
however,  has  just  begun. 

In  this  section,  I  outline  successive  background  context  which  leads  up  to  the  topic 
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of  the  whiskers  as  hydrodynamic  sensors:  1-  an  overview  of  all  of  the  seals’  sensory 
systems,  2-  the  structure  and  function  of  the  vibrissae  (whiskers),  3-  some  examples 
of  related  sensory  systems  in  other  animals,  and  4  -  the  types  of  flow  information  left 
by  fish.  Then  the  series  of  hydrodynamic  sensing  experiments  with  real  seals  that 
inspire  this  specific  study  is  described. 

1.1.1  Pinniped  Sensory  Systems 

As  with  many  animals,  the  sensory  systems  available  to  the  harbor  seal  are  rich  and 
diverse  [134,  110].  These  systems  may  enable  them  to  carry  out  impressive  feats 
observed  in  the  field,  such  as  precisely  navigating  across  long  distances  to  return  to 
their  breeding  area  yearly  [72]  and  relocating  breathing  holes  after  under-ice  excur¬ 
sions  [33]. 

The  following  section  presents  an  overview  of  the  known  functions  of  each  of  the 
seals’  sensory  systems  and  how  they  may  or  may  not  aid  hydrodynamic  tracking. 

Vision 

The  pinniped  vision  system  is  adapted  for  use  over  a  broad  range  of  light  intensities. 
The  iris  is  narrow  and  slit-like  in  daylight  but  dilates  to  a  circular  shape  as  ambient 
light  decreases  in  the  water  (Figure  1-3).  The  large  spherical  lens  allows  for  focusing 
of  images  on  the  retina  despite  the  loss  of  refractive  power  under  water.  Harbor  seals 
have  high  contrast  sensitivity,  which  enables  them  to  detect  movement  as  effectively 
as  cats  [53]. 

Behavioral  experiments  with  a  harbor  seal  and  California  sea  lion  show  them 
to  have  better  dark  adaptation  rates  and  relative  light  sensitivities  than  a  human 
subject.  And  an  elephant  seal,  which  is  a  deeper  diving  animal  than  the  others,  was 
measured  to  have  an  even  lower  threshold  estimate  [84], 

Conditions  of  high  turbidity  are  known  to  dramatically  reduce  visual  acuity  [136] . 
However,  a  recent  experiment  [47]  demonstrates  the  harbor  seal’s  potential  to  use 
particles  in  the  water  to  resolve  flow  direction  down  to  very  small  angles.  Despite 
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Figure  1-3:  A  harbor  seal’s  eye  covers  a  wide  pupillary  range  to  aid  visual  acuity  in 
both  water  and  air.  The  pupil  goes  from  (a)  circular  in  dark  water  to  (b)  a  vertical 
slit  on  land,  as  light  increases.  The  scale  bar  corresponds  to  1cm.  Figure  from  [53]. 

heightened  visual  perception,  seals  should  only  be  able  rely  on  this  to  locate  prey 
when  they  are  in  relatively  close  proximity.  The  question  still  remains  as  to  how 
these  animals  can  hunt  from  long  distances  away. 

Echolocation 

Pinnipeds  exhibit  a  range  of  vocalizations  [116].  The  pinniped  auditory  system, 
however,  does  not  contain  the  same  structural  adaptations  for  localizing  underwater 
sounds  as  seen  in  other  animals  (in  ringed  seals  the  auditory  bulla  is  not  acoustically 
isolated  from  the  skull  [69]). 

Other  animals,  such  as  dolphins,  are  known  to  have  an  active  sonar  system  [134], 
which  they  use  to  detect  prey,  communicate,  etc.  [105]  observed  sea  lions  to  emit 
sound  as  they  approached  pieces  of  fish  in  a  pool.  The  animals’  success  in  catching 
the  fish  regardless  of  the  light  conditions,  led  the  author  to  suggest  that  pinnipeds 
may  also  use  an  echo-ranging  sense  to  locate  prey. 

Over  the  next  few  decades,  multiple  studies  proceeded  to  provide  conflicting  evi¬ 
dence  for  pinniped  echolocation.  [Ill]  again  showed  pinnipeds  able  to  locate  objects 
in  total  darkness.  Moreover,  the  harbor  seals  here  were  able  to  discriminate  between 
two  objects,  which  were  visually  identical  but  differed  in  acoustic  impedance  one  air- 
filled  ring  and  one  water-filled  ring  (Figure  l-4a).  However,  since  neither  the  animals’ 
sound  emission  nor  reception  were  blocked  as  a  control,  reliance  on  echolocation  in 
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Figure  1-4:  Sketches  of  (a)  an  echolocation  experiment  and  (b)  a  magnetic  sensing 
experiment  conducted  with  harbor  seals.  Figure  from  [109]. 

this  case  could  not  be  unequivocally  stated. 

On  the  other  hand,  [98]  showed  a  negative  correlation  between  sound  emission  and 
navigational  ability,  [33]  recorded  no  vocalizations  during  a  search  for  ice  holes,  and 
[18]  recorded  very  few  vocalizations  by  Weddell  seals  during  under-ice  prey  attacks. 

[119]  argues  against  echolocation  in  pinnpeds,  pointing  out  that  these  animals  have 
not  developed  the  high-frequency  sound  production  or  reception  systems  required  for 
accurate  underwater  echolocation.  Most  recorded  sounds  in  the  wild  were  emitted 
at  low  frequencies  [112],  suggesting  they  are  related  to  social  interaction,  rather  than 
prey  detection.  [119]  also  notes  that  echolocation  is  a  result  of  adaptation,  rather 
than  a  generalized  skill.  Examples  of  opportunistic  echolocation  may  exist,  however, 
this  does  not  imply  that  they  have  evolved  an  accurate  echolocation  sensory  system. 

The  evidence  against  echolocation  now  surpasses  that  for  it.  The  frequency  and 
intensity  of  vocalizations  in  seals  are  low  compared  to  those  in  dolphins,  and  their 
ears  have  no  specific  adaptations  that  allow  for  better  functioning  in  the  two  media 
of  water  and  air.  It  is  therefore  unlikely  that  seals  would  rely  heavily  on  echolocation 
to  detect  prey. 
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Biomagnetism 


[109]  describes  a  test  of  a  harbor  seal’s  ability  to  detect  changes  in  a  magnetic  field. 
The  test  subject  was  asked  to  swim  through  a  Helmholtz  coil,  which  shifts  the  mag¬ 
netic  field  (Figure  l-4b).  In  these  test  conditions,  the  seal  exhibited  poor  ability  to 
detect  the  change. 


Olfaction 

There  is  evidence  for  the  olfactory  system  playing  a  role  in  social  interactions.  Its 
use  in  sensing  may  also  be  possible.  There  are  a  few  demonstrations  of  olfactory 
sensations  being  used  by  animals  underwater  [89],  including  harbor  seals  detecting 
dimethyl  sulfide,  which  can  indicate  profitable  foraging  areas  [78] .  However,  pinnipeds 
are  generally  considered  microsmatic  (having  a  poor  sense  of  smell)  [88]. 


Gustation 

The  number  of  taste  buds  in  pinnipeds  is  less  than  that  on  terrestrial  animals,  which 
suggests  a  lower  sense  of  taste.  A  study  by  [125],  nonetheless,  shows  harbor  seals  as 
being  able  to  use  gustation  to  discriminate  small  changes  in  salinity  levels  in  seawater. 

With  the  exception  of  echolocation,  none  of  the  sensory  systems  mentioned  up  to 
this  point  appear  to  be  capable  of  enabling  wake  detection  over  the  long  distances 
(tens  of  meters)  reported  in  multiple  studies  by  the  University  of  Rostock  research 
group.  The  seal  would  need  to  be  within  close  range  of  the  prey  in  order  to  successfully 
use  any  of  these. 

1.1.2  Vibrissae  (Whiskers) 

General  information  on  the  vibrissal  system  is  now  provided  to  help  determine  how 
it  may  or  may  not  be  able  to  play  a  role  in  sensing. 
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Why  they  might  be  important 


Vibrissae  are  found  in  almost  all  mammals  except  Homo  sapiens,  and  they  are  perhaps 
most  prominent  on  pinnipeds  [87] .  There  are  several  factors  that  point  to  the  vibrissal 
system  as  an  important  one  for  pinnipeds.  These  include  their  number,  size,  growth 
rates,  and  thermal  regulation  system. 

Adult  harbor  seals  have  about  40  vibrissae  next  to  either  nostril.  Their  lengths 
vary,  but  can  reach  up  to  about  125  mm  long  [117],  making  them  among  the  largest 
known  animal  structures  that  could  serve  as  flow  sensors. 

[62]  and  [51]  measured  the  growth  rates  of  harbor  seal  and  gray  seal  vibrissae, 
respectively.  The  vibrissae  of  both  species  were  found  to  not  shed  in  the  annual  molt 
with  the  rest  of  the  pelage.  They  are  shed  irregularly  and  individually.  This  suggests 
that  seals  have  developed  a  means  of  ensuring  continuous  presence  of  select  vibrissae 
for  sensing. 

Single  vibrissal  follicles  are  found  to  be  units  of  high  thermal  radiation.  It  is 
suggested  that  by  keeping  the  vibrissal  follicle  temperature  constant,  the  vibrissal 
system  can  maintain  high  tactile  sensitivity  regardless  of  any  decrease  in  ambient 
temperature  [90,  34],  Figure  1-5  shows  the  elevated  temperature  at  each  individual 
follicle. 

Known  functions:  tactile  sensing  and  others 

Vibrissae  are  known  to  serve  multiple  behavioral  roles,  such  as  communicating  socially 
(Figure  1-6),  monitoring  environment  parameters,  exhibiting  aggression,  etc.  [1].  It 
has  long  been  suggested  that  mammals  who  must  hunt  nocturnally  may  use  their 
whiskers  to  aid  hunting  [68],  though  the  means  by  which  they  do  this  is  unknown. 
It  has  also  been  speculated  that  the  whiskers  may  serve  as  a  filter  for  feeding  off  the 
ground  [116]. 

The  most  known  function  of  the  vibrissae  as  a  sensory  system  is  that  of  active 
touch  discrimination.  Several  studies  have  shown  vibrissae  to  be  sensitive  to  even 
slight  bending  as  they  contact  physical  objects.  [74]  and  [19]  were  the  first  to  demon- 


Figure  1-5:  (Left)  Thermogram  of  a  harbor  seal’s  face  immediately  after  leaving  the 
water.  The  warmest  places  are  around  the  eyes,  nostrils,  and  vibrissal  pads.  Figure 
from  [21].  (Right  top)  Photograph  and  (Right  bottom)  thermogram  of  the  vibrissal 
pad  of  a  harbor  seal  Figure.  The  follicles  have  higher  temperature  than  the  rest  of 
the  pad.  Figure  from  [90]. 


Figure  1-6:  A  naso-nasal  greeting  between  a  harbor  seal  mother  and  her  pup,  demon¬ 
strating  the  use  of  vibrissae  to  communicate  socially.  Figure  from  [110]. 
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Figure  1-7:  (a)  Different  shapes  that  a  sea  lion  was  asked  to  distinguish  using  vibrissal 
tactile  sensing,  (b)  Experimental  setup  for  the  blindfolded  sea  lion.  Figures  from  [19]. 


strate  this  tactile  sensitivity  in  pinnipeds.  Their  respective  experiments  showed  a 
Pacific  walrus  and  a  California  sea  lion  able  to  distinguish  objects  of  different  shapes 
despite  being  blindfolded  (Figure  1-7). 

Next,  tactile  sensing  was  demonstrated  to  enable  size  discrimination  of  objects. 
[20]  showed  a  blindfolded  California  sea  lion  able  to  identify  the  larger  of  two  cir¬ 
cular  discs.  The  detection  threshold  was  quantified  through  the  Weber  fraction 
(c  =  A  D/D,  D  =  disc  diameter),  which  indicates  the  relative  diameter  difference  de¬ 
tected  on  75%  of  occasions.  In  these  experiments,  the  mean  Weber  fraction  achieved 
was  c  =  0.26.  [23]  conducted  the  same  test  but  with  harbor  seals,  which  belong  to 
the  family  Phocidae,  rather  than  the  Otariidae  (California  sea  lion)  or  Odobenidae 
(Pacific  walrus)  families.  The  tested  harbor  seals  achieved  a  similar  Weber  fraction 
for  disc  diameters  of  the  same  size  used  in  [20].  For  larger  discs,  the  Weber  fraction 
reduced  to  as  low  as  c  =  0.08.  Those  tests  were  all  conducted  in  air,  so  [27]  tested 
the  harbor  seals’  tactual  resolution  ability  in  water.  No  significant  difference  was 
observed  between  water  and  air.  In  both  media,  the  seals  could  detect  size  differences 
as  low  as  2 mm. 

These  findings  reveal  the  tactile  acuity  of  harbor  seals  to  be  comparable  to  that 
of  human  subjects  measuring  the  diameter  of  objects  using  their  thumb  and  index 
finger  (c  ~  0.03)  [43]  and  is  nearly  as  good  as  their  own  visual  resolving  power  [26]. 
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Under  the  Skin:  Structure 


Vibrissae  are  different  from  pelage  hair  in  a  few  ways.  First,  the  base  of  the  vibrissae 
are  attached  to  striated  muscles  for  voluntary  control  [1],  Also,  each  vibrissa  is 
attached  to  an  innervated  follicle,  which  appears  to  have  similar  structure  across 
pinnipeds  and  cats  [70].  In  pinnipeds,  however,  the  follicle-sinus-complex  (FSC)  is 
the  largest  of  all  mammals.  They  can  reach  2 cm  in  length,  whereas  those  on  terrestrial 
mammals  of  similar  size  are  much  shorter  (for  example,  those  in  the  leopard  Panthera 
pardus  are  6.5 mm  long)  [22],  A  schematic  and  photo  of  the  FSC  are  shown  in  Figure 
1-8.  The  hair  shaft  is  encapsulated  in  three  sections:  the  upper  cavernous  sinus,  the 
ring  sinus,  and  the  lower  cavernous  sinus.  The  bottom  two  (RS  and  LCS)  contain 
mechanoreceptors.  The  top  part  (UCS),  which  is  significantly  shorter  or  nonexistent 
in  mammals  other  than  pinnipeds,  may  serve  to  distance  the  receptors  from  low 
ambient  water  temperatures  [69,  90].  No  nerves  are  present  inside  the  hair  shaft. 
Rather,  at  the  base  of  the  FSC,  a  nerve  bundle  penetrates  into  the  capsule.  This 
leads  away  to  a  specific  location  in  the  brain’s  sensory  cortex. 

In  [79],  different  parts  of  a  fur  seal  were  stimulated  while  the  activity  of  cortical 
neurons  was  recorded  via  microelectrodes  inserted  in  the  somatosensory  cortex.  Ex¬ 
citation  of  the  vibrissae  was  found  to  stimulate  a  large  portion  of  the  somatosensory 
cortex  (Figure  1-9).  The  area  occupied  by  projections  of  the  head  are  greater  than  the 
area  occupied  by  the  trunk  and  limbs  combined.  It  is  well  established  that  the  relative 
size  of  projections  indicates  the  importance  of  the  tactile  information  received  from 
those  zones.  At  the  time  of  that  study,  no  other  equal  to  this  vibrissal  representation 
had  been  seen  in  other  mammals.  Also,  despite  the  variation  in  size  of  vibrissae,  they 
all  take  up  approximately  the  same  area  in  the  cortex,  suggesting  that  each  is  equally 
useful. 

The  nerves  penetrating  the  FSC  in  pinnipeds  have  an  order  of  magnitude  higher 
innervation  than  those  in  terrestrial  mammals  (1,000-1,500  nerves  to  each  vibrissae 
in  ringed  seals  versus  100-120  in  beavers  and  rats)  [71].  The  same  comparison  can 
be  made  with  the  mechanoreceptor  density  in  the  ring  sinus  of  the  FSC,  where  the 
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(a)  (b) 


Figure  1-8:  (a)  A  diagram  of  the  follicle-sinus-complex  system,  which  holds  each  vib¬ 
rissa  and  connects  it  to  the  brain,  (b)  A  ringed  seal’s  FSC  partially  split  open.  1: 
upper  cavernous  sinus,  2:  ring  sinus,  3:  lower  cavernous  sinus,  C:  capsule,  Cm:  col¬ 
lagenous  membrane,  R:  liquid-filled  cleft,  Rw:  ringwulst,  N:  nerve  bundle  penetrating 
the  capsule,  G:  groove.  Figures  from  [69]. 


Figure  1-9:  (a)  The  somatosensory  cortex  portion  of  the  brain. 

The  points  on  the  body  that  stimulate  the  labeled  regions  on  the  cortex  from  (b) 
the  trunk,  (c)  the  head,  (d)  the  vibrissae  (marked  in  size  according  to  length  and 
diameter).  The  vibrissae  are  seen  to  take  up  a  sizable  portion  of  the  total  area. 
Figure  from  [79]. 
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number  of  Merkel  cells  in  a  ringed  seal’s  ring  sinus  is  2-5  times  greater  than  that  in 
cats  [69,  70].  The  increased  innervation  density  again  suggests  an  important  sensory 
role  for  pinniped  vibrissae. 

Under  the  Skin:  Nerve  Response 

The  term  vibrissa  comes  from  the  Latin  word  “vibrio” ,  which  means  to  vibrate.  Stud¬ 
ies  have  been  carried  out  to  determine  how  the  nerve  fibers  respond  to  vibrations,  as 
that  will  define  the  scope  within  which  vibrissae  can  serve  to  relay  sensory  informa¬ 
tion. 

These  nerves  are  composed  of  afferent  (indicating  direction  is  from  the  receptor  to 
the  central  nervous  system)  fibers  that  respond  to  deflections  of  the  hair  shaft.  Two 
classes  of  afferent  fibers  are  distinguished:  slowly  adapting  (SA)  and  rapidly  adapt¬ 
ing  (RA).  The  former  discharges  throughout  the  duration  of  a  stimulus.  The  latter 
provides  a  nerve  impulse  only  during  a  movement  and  therefore  does  not  respond  to 
steady  displacement.  It  is  assumed  that  SA  fibers  provide  position  information  and 
RA  fibers  provide  vibration  information. 

[30]  and  [31]  attempted  to  identify  the  range  of  stimuli  that  produce  neural  activity 
in  the  afferent  fibers.  The  author  stimulated  individual  afferent  fibers  on  vibrissae 
of  harbor  seals  and  cats,  and  he  surgically  exposed  the  infraorbital  branch  of  the 
trigeminal  nerve  to  measure  the  response.  Seals  were  found  to  have  a  somewhat 
higher  percentage  of  RA  fibers  than  cats.  Otherwise,  the  two  species  showed  minor 
functional  differences.  Their  results  suggest  that  harbor  seals  are  capable  of  detecting 
vibrations  between  5-7 00 Hz  and  with  a  minimum  threshold  of  about  10 urn.  For  some 
fibers,  the  discharge  frequency  was  unaffected  by  the  stimulus  amplitude.  In  both 
species,  the  response  of  the  SA  fibers  saturated  between  12-20°  of  vibrissae  deflection, 
indicating  their  sensitivity  to  minor  hair  deflections,  rather  than  large  ones.  The  RA 
fibers  showed  strong  ability  to  detect  small  vibratory  disturbances.  The  author  notes 
that  these  properties  should  aid  tactile  detection  of  texture  and  shape. 

Later  studies  [107,  92]  measure  the  smallest  amplitude  that  seals  may  perceive 
while  their  vibrissae  are  vibrating  (in  the  range  of  50-1000 H z) .  Harbor  seals  were 
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Figure  1-10:  Sketch  of  experiment,  in  which  a  harbor  seal  was  asked  to  identify 
whether  or  not  a  rod  in  contact  with  its  vibrissae  was  vibrating.  Figure  from  [109]. 

trained  in  air  to  touch  their  whiskers  to  a  rod  (Figure  1-10)  and  indicate  whether  or 
not  the  rod  was  vibrating.  Unlike  the  previous  studies,  these  results  indicate  higher 
sensitivities  at  higher  frequencies  and  an  ability  to  detect  displacements  <  Ijirn. 

The  natural  frequencies  of  whiskers  shed  by  fur  and  harbor  seals  were  measured 
[121,  59].  They  show  a  first  mode  of  vibration  at  around  A()-&)Hz.  This  is  on  the 
lower  end  of  the  spectra  tested  above.  The  question  of  interest  for  wake  following 
would  be  whether  or  not  the  vibrissae  are  sensitive  to  the  low-frequency  signals  that 
would  be  left  in  their  prey  trails  (0(1) Hz). 

Morphology 

Most  phocid  seals  exhibit  a  different  whisker  morphology  than  other  pinnipeds  [86]. 
Elephant  seals  (Mirounga  spp.),  gray  seals  (Halichoerus  grypus),  harbor  seals  (Phoca 
vitulina),  harp  seals  (Phoca  groenlandica) ,  hooded  seals  (Cystophora  cristata),  and 
ringed  seals  (Phoca  hispida),  have  an  undulatory  and  flattened  shape.  Other  pin¬ 
nipeds,  such  as  bearded  seals  (Erignathus  barbatus),  monk  seals  (Monachus  spp.), 
those  in  the  Odobenidae  family  (walruses),  and  those  in  the  Otariidae  family  (fur 
seals  and  sea  lions)  have  a  smooth  profile.  Sketches  of  the  different  morphologies  are 
shown  in  Figure  1-11  and  an  image  comparing  the  ringed  seal’s  and  sea  lion’s  whiskers 
is  shown  in  Figure  1-12. 

In  [46]  the  bumpy  morphology  of  the  vibrissae  on  harp,  hooded,  and  grey  seals  were 
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Gray  seal 


Figure  1-11:  A  few  of  the  different  vibrissal  morphologies  seen  in  pinnipeds.  Phocid 
seals  generally  exhibit  an  undulatory  profile  unlike  the  smooth  profile  seen  on  otariids. 
Figure  from  [45]. 


Crest  Width 


_ _  \ : 

Peak-to-Peak  Trough  Width 


Figure  1-12:  Pictures  of  a  vibrissa  from  a  ringed  seal  (top)  and  California  sea  lion 
(bottom).  Figure  from  [45]. 
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Figure  1-13:  The  number  of  bumps  is  seen  to  increases  somewhat  linearly  with  vibrissa 
length.  Figure  from  [46]. 

optically  measured.  The  bearded  seals’  smooth  vibrissae  were  used  for  comparison. 
The  wavelength  of  the  bumps  is  relatively  constant,  as  the  number  of  bumps  increases 
fairly  linearly  with  the  total  whisker  length  (Figure  1-13).  From  their  data,  the 
X/d  ratio  is  about  5  to  6,  where  A  is  the  undulation  period  and  d  is  the  average 
diameter  on  the  long  side  of  the  whisker.  They  note  that  the  shortest  vibrissae 
(length  <  25 mm)  had  no  bumps,  indicating  that  there  is  some  development  aspect 
to  the  bumpy  morphology.  Also,  the  width  of  the  vibrissae  tapers  towards  the  tip. 

In  [45],  morphological  differences  between  phocids  with  the  beaded  vibrissal  pro¬ 
file,  phocids  with  a  smooth  profile,  and  otariids  with  a  smooth  profile  were  compared. 
Using  a  combination  of  traditional  and  morphometric  analyses,  100%  of  vibrissae 
were  able  to  be  classified  according  to  the  species.  It  is  therefore  shown  that  these 
three  groups  occupy  distinct  morphospaces  (based  on  peak-to-peak  distance,  crest 
width,  and  trough  width). 

The  harbor  seal’s  vibrissal  geometry  is  described  by  [45]  as  having,  on  average,  a 
peak-to-peak  distance  of  3.3  mm,  crest  width  of  0.9  mm,  trough  width  of  0.73  mm, 
total  length  of  86.5  mm,  and  2.3  beads/cm.  The  three-dimensionality  is  described 
in  [58],  though  they  report  a  significantly  shorter  undulation  wavelength  than  that 
in  [45].  The  latter  matches  what  I  measured  from  vibrissae  provided  by  the  New 
England  Aquarium. 

A  variety  of  foraging  methods  are  employed  across  these  species.  It  is  therefore 
conceivable  that  their  vibrissal  morphologies  could  have  evolved  to  serve  slightly  dif- 
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ferent  functions.  Otariids  generally  do  not  hunt  in  as  deep  water  as  phocids  do,  and 
therefore  should  not  need  to  rely  as  much  on  non-visual  cues.  Likewise,  their  vibris- 
sae  perhaps  need  not  have  developed  the  extra  feature  of  undulations,  as  those  on 
phocids  have. 


1.1.3  Other  Animals 

Several  other  animals,  which  have  similar  systems  to  the  pinniped  vibrissae,  will  now 
be  briefly  over  viewed. 


With  hairlike  structures 

Many  marine  mammals  have  a  vibrissal  system  [87].  Baleen  whales  have  about  100 
vibrissae  around  their  blowholes  and  jaw.  Odontocetes  (dolphins,  etc.)  have  vibrissae, 
but  most  lose  them  postnatally.  Manatees  have  vibrissae  across  their  body  and  have 
been  shown  capable  of  detecting  hydrodynamic  dipole  stimuli  using  those  on  their 
face  [106,  42],  They  are  herbivorous,  however,  so  it  is  less  likely  that  they  would  need 
to  employ  wake  tracking.  Several  terrestrial  mammals  also  have  vibrissae  (cats,  rats, 
insects,  etc.).  Studies  on  some  of  these  have  shown  their  tactual  function  [38],  and 
have  inspired  creation  of  artificial  tactile  sensors  [122],  Figure  1-14  displays  a  few 
animals  with  vibrissae. 

There  are  also  animals  with  non- vibrissal  structures,  which  are  nonetheless  hair  like 
and  can  be  used  for  sensing  (Figure  1-15).  For  example,  moles,  shrews,  and  aquatic 
snakes  conduct  tactile  sensing  using  such  structures  [15,  14].  Also,  the  lateral  line 
organ  on  fish  consists  of  a  set  of  neuromasts.  These  are  a  group  of  thousands  of 
sensory  hair  cells  with  kinocillia  and  stereovilli  that  project  into  a  flexible  cupula 
(tens  to  hundreds  of  micrometers  in  size) .  Deflection  of  this  structure  provides  sensory 
information  to  the  fish. 
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Figure  1-14:  Examples  of  other  animals  with  vibrissae:  (a)  manatee  (figure  from 
[106]);  (b)  harbor  porpoise  (figure  from  [16]);  (c)  rat  (figure  from  [61]). 


Figure  1-15:  Examples  of  other  animals  with  hairlike  structures:  (a)  star-nosed  mole 
and  (b)  tentacled  snake  (both  figures  from  [14]). 


Other  types  of  sensing 

There  are  many  examples  of  non-vibrissal  based  sensing  in  animals.  A  few  are  listed 
here  to  give  a  sense  of  the  vast  range  of  sensory  capabilities  in  the  animal  kingdom. 
Bloodhounds  and  other  dogs  use  olfaction  to  detect  the  direction  of  a  track  [128]. 
Copepods  can  conduct  chemical  sensing.  Snakes  can  use  their  tongues  to  transport 
odors  and  follow  scent  trails  [120,  135].  Sharks  can  detect  electric  fields  [73].  The  star¬ 
nosed  mole  and  water  shrew  are  able  to  sniff  and  detect  semiochemicals  underwater 
[89].  Crocodiles  have  dome-shaped  pressure  receptors  distributed  across  their  face. 


Other  instances  of  tracking 

There  are  a  limited  number  of  studies  of  animals  conducting  wake  tracking.  One 
of  the  few  examples  is  that  of  a  catfish  tracking  a  small  guppy.  The  test  subject’s 
vision  was  obscured  by  shining  infrared  light  [104]  and  its  taste  sensation  was  ablated 
[103] .  Ablation  of  taste  was  not  seen  to  lower  the  effectiveness  of  guppy  capture,  but 
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Figure  1-16:  Sketch  of  the  reverse  von  Karman  vortex  wake  left  by  a  swimming  fish, 
which  produces  thrust.  Figure  from  [29]. 

restriction  of  the  lateral  line  did.  At  first  glance,  the  prominence  of  their  “whiskers” 
suggest  they  could  play  a  role  in  hydrodynamic  sensing.  However,  these  appendages 
have  taste  buds,  so  are  likely  used  in  chemical  sensing.  Also,  these  barbels  are  fleshy, 
unlike  the  stiff  vibrissal  hair  shafts,  which  sit  in  an  innervated  FSC.  It  should  be 
noted  that  this  tracking  was  done  over  relatively  short  distances  (50  guppy  lengths) 
and  times  (lOsec),  which  still  leaves  the  performance  of  the  harbor  seal  as  a  yet 
unrivaled  and  unexplained  feat. 

1.1.4  Hydrodynamic  Stimuli 

The  focus  of  this  thesis  revolves  around  detection  of  hydrodynamic  information.  This 
section  provide  details  on  what  information  may  be  left  in  the  wakes  of  swimming 
animals  for  the  seals  to  sense. 

The  hydrodynamic  trails  produced  by  a  swimming  fish  can  take  many  forms, 
depending  on  various  parameters  that  define  their  kinematics  (pitch  angle,  heave 
amplitude,  flapping  frequency,  etc.).  It  is  beyond  the  scope  of  this  thesis  to  describe 
multiple  types  of  wakes  and  study  seals’  responses  to  them.  Instead,  I  chose  to  focus 
on  the  type  of  wake  left  during  efficient  swimming,  which  is  that  of  the  reverse  von 
Karman  vortex  street  (Figure  1-16).  A  vortex  is  shed  with  each  flick  of  the  fin,  which 
results  in  an  array  of  counter-rotating  vortices  and  a  mean  flow  that  produces  thrust. 

Studies  by  [55]  and  [54]  provide  some  information  on  how  long  fish  wakes  can  last. 
In  [55],  even  the  wakes  of  small  (10cm  long)  goldfish  ( Carassius  auratus )  could  still  be 
measured  5 min  after  the  fish  swam  by  in  the  low  background  flow  lab  conditions  of  the 
study.  [54]  measured  the  trails  of  three  teleost  species  (38-8 6mm  long).  Typical  flow 
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velocities  measured  60sec  after  the  wake  generating  fish  had  passed  were  Ylmmjs. 
The  authors  state  that  this  implies  that  a  predator  could  detect  the  fish  at  a  distance 
of  tens  of  meters.  A  harbor  seal  itself  was  measured  to  leave  flow  velocities  of  at 
least  30 mm/ s  after  30s  [118].  Based  on  extrapolated  values,  [97]  proposed  that  water 
movements  caused  by  the  C-starts  of  teleost  fish  may  last  for  several  minutes,  even  if 
background  noise  is  present. 

An  efficiently-swimming  fish  will  flap  its  tail  at  a  nondimensional  frequency  de¬ 
scribed  as  the  Strouhal  number  (St  =  fw/U )  of  ~  0.3,  where  w  is  the  wake  width, 
U  is  the  swimming  speed,  and  /  is  the  flapping  frequency  [129].  The  size  and  speed 
of  fish  pursued  by  harbor  seals  should  therefore  produce  vortex  wakes  with  dominant 
frequencies  of  0(1  —  10 )Hz. 

1.1.5  Hydrodynamic  Perception  Experiments  with  Pinnipeds 

In  the  1950s,  and  perhaps  earlier,  it  was  suggested  that  seals  may  use  their  vibrissae 
to  detect  water  changes  as  a  cue  for  prey  detection  [60].  A  field  observation  of  a 
blind  grey  seal,  who  was  nevertheless  fat  and  healthy,  indicated  its  ability  to  catch 
food.  The  author  ruled  out  the  possibility  of  sight,  smell,  and  hearing  and  concluded 
that  its  vibrissae  must  have  been  used:  “One  may  suppose  that  by  means  of  them 
[vibrissae]  the  animal  is  made  aware  of  slight  turbulences  in  the  water  such  as  might 
be  produced  by  the  swimming  of  a  fish.” 

Since  the  1990s,  the  Marine  Science  Center  group  in  Rostock  has  led  a  series 
of  experiments  with  pinnipeds,  demonstrating,  in  particular,  harbor  seals’  ability  to 
perceive  hydrodynamic  information.  A  review  of  these  experiments  has  been  captured 
in  a  few  places  [57,  143,  56,  21]  and  an  overview  is  provided  here. 

Dipole  Stimuli  Detection 

[24]  was  the  first  behavioral  test  of  a  seal’s  ability  to  use  its  whiskers  to  detect  hy¬ 
drodynamic  stimuli.  A  similar  setup  as  has  often  been  used  to  test  the  lateral  line 
sensing  capabilities  of  a  fish  was  employed  here  (Figure  1-17).  An  oscillating  sphere 
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Figure  1-17:  Experiment  setup  for  dipole  stimuli  detection.  A  harbor  seal  wearing 
a  blindfold  and  earmuffs  was  stationed  at  a  set  distance  from  an  oscillating  sphere. 
Vibrissal  detection  threshold  of  the  resulting  water  movements  was  measured.  Figure 
from  [26]. 

generated  a  hydrodynamic  stimulus  over  a  range  of  lO-lOOiTz.  The  seal  was  trained 
to  rest  its  jaw  on  a  station  to  maintain  a  specific  distance  between  it  and  the  sphere. 
Optical  and  acoustic  cues  were  removed  by  using  eye  caps  and  headphones  with  pink 
noise  masking.  The  seal  was  trained  to  leave  the  hoop  station  upon  detection  of  the 
stimulus. 

The  distance  between  the  stimulus  and  the  whiskers  was  adjusted  and  used  to 
calculate  an  effective  stimulus  amplitude.  The  best  response  was  seen  at  50 Hz,  which 
corresponded  to  ~  250 fim/s.  At  that  frequency,  the  absolute  detection  threshold  of 
water  displacement  was  0.8/im.  The  seal  was  also  tested  while  wearing  a  wire  mesh 
muzzle,  which  blocked  whisker  movement.  In  those  cases,  he  never  responded  to  the 
stimuli.  These  results  indicate  the  harbor  seal’s  ability  to  detect  minute  movements 
using  its  whiskers,  though  the  detection  threshold  is  not  as  sensitive  as  that  seen  in 
some  fish  species. 

Toy  Submarine  Trail 

Next,  seals’  ability  to  follow  the  trail  left  by  a  moving  object  was  demonstrated 
[25].  The  authors  selected  a  propeller-driven  toy  submarine  of  a  size  that  should 
create  trails  with  similar  flow  velocity  and  spatial  extent  similar  to  those  wakes  of 
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Figure  1-18:  Schematic  of  the  submarine  trail  detection  experiment.  A  harbor  seal 
wearing  a  blindfold  is  asked  to  locate  a  toy  submarine,  which  swam  ahead  on  a  linear 
path  in  a  random  direction  several  seconds  earlier.  The  yellow  line  marks  the  actual 
trail,  the  red  line  indicates  a  potential  search  path  the  seal  could’ve  taken,  and  the 
blue  semicircle  marks  the  locations  where  the  sub  could  have  landed.  The  probability 
of  the  seal  finding  the  submarine  by  arbitrarily  choosing  a  direction  is  P  <  0.04. 
Figure  from  [25]. 

swimming  fish  measured  in  [55] .  Again  the  seal  was  blindfolded  by  wearing  a  stocking 
mask.  Headphones  were  used  for  acoustical  masking  while  the  submarine  swam  ahead 
(for  3-5sec).  About  two  seconds  after  the  submarine  stopped,  the  headphones  were 
removed  and  the  seal  began  its  search. 

The  sub  was  commanded  to  follow  a  linear  path  in  a  random  direction  (Figure 
1-18).  The  seal  began  its  search  by  protracting  its  vibrissae  to  the  forward  position 
and  slightly  moving  its  head  laterally  until  it  detected  the  trail,  at  which  time  it  swam 
towards  the  submarine.  Even  with  a  delay  (up  to  20sec)  between  starts  of  the  two 
objects,  the  seal  successfully  found  the  submarine  ~  80%  of  the  time.  The  authors 
note  that  this  20sec-old  wake  is  equivalent  to  the  seal  having  chased  the  object  from 
40  m  away. 

A  spontaneous  change  in  course  was  then  introduced  into  the  submarine’s  path. 
As  the  seal  encountered  this  change,  it  sharply  turned  to  stay  on  the  trail.  On  the 
other  hand,  if  the  engine  was  allowed  to  run  for  a  short  time  after  the  earmuffs 
had  been  removed,  the  seal  would  follow  a  direct  path  to  the  source  (Figure  1-19). 
Finally,  an  experimentally  naive  harbor  seal  was  able  to  carry  out  these  tasks  with 
similar  performance,  further  solidifying  evidence  for  the  harbor  seal’s  hydrodynamic 
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Figure  1-19:  The  two  paths  taken  by  the  harbor  seal  (red  arrow)  to  locate  the  toy 
submarine  (yellow  arrow).  In  the  absence  of  acoustical  cues,  the  seal  followed  the 
sharp  turn  that  the  sub  also  took  (solid  yellow  line).  When  the  sub  engine  was 
audible,  the  seal  followed  a  direct  path  to  it  (dotted  yellow  line).  Figure  from  [25]. 

perception. 

Biogenic  Trail 

[118]  showed  the  accuracy  and  success  of  trail  following  to  be  even  slightly  higher 
(success  rate  ~  90%)  when  the  seal  followed  a  second  seal  (biogenic  tracking).  A  sim¬ 
ilar  experimental  setup  was  used,  with  the  test  seal  wearing  a  blindfold  and  earmuffs 
until  the  wake  generator  reached  its  destination. 

Trail  generation  lasted  a  mean  of  13.2sec.  The  blindfolded  seal  was  able  to  closely 
follow  the  path  of  the  front  seal,  never  deviating  more  than  one  body  length  from 
it.  It  was  seen  to  employ  two  types  of  search  patterns:  in  ~  65%  of  the  cases,  it 
followed  a  linear  pattern  and  in  the  remaining  ~  35%  of  the  cases,  the  seal  repeatedly 
crisscrossed  the  path  (Figure  1-20).  The  authors  postulate  that  this  may  be  related 
to  the  angle  at  which  the  seal  approached  the  wake:  At  acute  angles,  it  could  stay 
head-on  with  the  trail;  from  more  obtuse  angles,  it  may  have  overshot  and  needed 
to  correct.  Control  experiments,  where  the  seal  wore  a  stocking  mask  muzzle,  were 
unsuccessful,  as  the  mask  irritated  the  seal  and  he  refused  to  perform  the  search. 

Artificial  Fin  Trail 

Then  [138]  tested  a  harbor  seal’s  ability  to  determine  trail  direction,  based  on  the 
wake  left  by  a  hnlike  paddle.  Two  paddles  were  placed  at  either  end  of  an  enclosed 
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Figure  1-20:  (A)  One  seal  swims  ahead  on  the  blue  path.  A  second  seal  swims  behind 
several  seconds  later  and  takes  the  red  path,  despite  being  blindfolded.  (B)  Same 
setup,  but  the  second  seal  now  employs  a  more  distinctively  undulatory  pattern  to 
the  search.  Figure  from  [118]. 

tank.  The  experimenter  moved  one  across  the  tank  (U  =  0.2 m/s),  and  when  the 
paddle  reached  the  other  end,  the  seal  had  to  identify  which  direction  the  movement 
had  been  in.  Based  on  a  metric  of  70%  correct  choices,  they  show  the  seal  is  able  to 
detect  this  motion  up  to  35sec  after  the  wake  was  created. 

It  should  be  noted  that  a  rigidly  clamped  fin  does  not  generate  a  reverse  von 
Karman  vortex  street,  as  a  swimming  fish  does.  So  the  manner  in  which  the  seal 
determined  direction  in  this  study  may  not  match  that  used  by  seals  chasing  fish. 
Based  on  the  information  presented  so  far,  it  is  difficult  to  say  how  the  specific 
flow  structures  may  have  aided  detection.  In  this  case,  and  in  the  previous  two 
tracking  studies,  it  is  possible  that  the  difference  in  the  mean  flow  velocity  of  the  wake 
compared  with  the  ambient  flow  was  sufficient  information  for  the  seal  to  determine 
the  direction  of  the  wake-generator. 

Impact  of  Glide  Phases 

Fish  often  do  not  swim  continuously,  rather  they  perform  alternating  burst  and  glide 
phases,  which  may  improve  swimming  energetics.  [139]  tested  the  impact  this  type 
of  swimming  would  have  on  seals’  detectability  of  the  wake.  A  submarine  was  com- 
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Figure  1-21:  Different  shapes  that  a  seal  was  trained  to  distinguish  based  solely  on 
the  wake  left  while  they  were  held  fixed  and  towed.  Figure  from  [140]. 

manded  to  run  continuously  for  3-5 m  then  glide  for  another  3-5 m.  The  seal  was  able 
to  track  this  wake  with  about  the  same  reliability  as  a  continuous  wake,  unless  the 
delay  was  longer  than  15sec. 

Object  Size  and  Shape 

[140]  showed  the  seals’  ability  to  distinguish  certain  sizes  and  shapes  of  objects  based 
on  their  hydrodynamic  signature.  Chemically  inert  paddles  of  different  shapes  (Figure 
1-21)  and  sizes  were  towed  in  front  of  a  blindfolded  seal.  There  was  a  three-second 
delay  between  the  models  reaching  their  endpoint  and  the  seal  being  allowed  to  enter 
the  tank. 

Based  on  the  wake  left  by  these  different  objects,  harbor  seals  are  able  to  detect  a 
difference  in  object  sizes  of  at  least  4cm.  This  is  not  as  sensitive  as  the  fish’s  lateral 
line  detection  in  terms  of  dimensional  value  [132] ,  however  it  would  likely  be  sufficient 
resolution  for  seals,  whose  normal  prey  size  is  at  least  10cm  in  length  [12]. 

Again,  it  is  difficult  to  say  what  specific  pieces  of  information  from  the  wake  the 
seal  may  have  used  in  this  case  to  conduct  detection. 

Comparison  with  Sea  Lion 

[48]  provides  the  first  behavioral  evidence  that  the  shape  of  the  whisker  may  play  a 
role  in  hydrodynamic  detection.  Using  similar  experimental  conditions  as  in  [25],  a 
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miniature  submarine  was  released  ahead  of  the  blindfolded  test  subject.  This  sub  was 
remotely  controlled  and  commanded  to  take  various  paths  with  curves. 

In  this  case,  the  experiment  was  conducted  on  a  California  sea  lion,  which  has  a 
different  vibrissal  morphology  than  that  of  the  harbor  seal.  The  former  has  a  smooth 
profile  and  the  latter  has  a  “beaded”  structure  as  described  in  Section  1.1.2. 

The  success  rate  (~  85%)  was  very  similar  to  that  of  the  harbor  seal  in  the  previous 
study,  and  the  test  subject  again  followed  any  curves  present  in  the  path.  However, 
unlike  the  harbor  seal,  if  the  sea  lion  had  to  wait  more  than  5sec  after  the  submarine 
had  left,  it  mostly  lost  the  trail. 

An  individual’s  performance  does  not  necessarily  represent  that  of  an  entire  species, 
however  these  results  may  suggest  that  sea  lions’  reduced  trail- following  ability  com¬ 
pared  with  harbor  seals  may  be  related  to  their  different  whisker  morphologies. 


1.1.6  What  is  the  mechanism  of  detection? 


These  demonstrations  naturally  raise  the  question  of  how  seals  can  do  this.  What 
enables  harbor  seals  to  detect  these  trails  using  their  whiskers  alone?  As  noted  in 
[84] ,  it  is  understandable  that  high  vibrissal  sensitivity  could  enable  heightened  tactile 
sensing,  which  can  be  used  for  “sit-and-wait”  foraging  strategies  or  burrowing  in  the 
ground.  But  how  can  the  vibrissae  be  used  while  the  seal  is  actively  chasing  prey? 
The  wakes  left  by  prey  should  be  difficult  to  detect  above  background  noise  while  the 
seal  is  swimming. 

This  is  the  question  explored  in  this  thesis.  A  key  hint  was  revealed  in  [58],  which 
showed  a  potential  connection  between  vibrations  and  sensing.  The  following  section 
provides  background  on  vortex-induced  vibrations,  which  are  the  underlying  physics 
that  will  help  reveal  the  answer. 
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Figure  1-22:  Ocean  structures,  such  as  offshore  platforms  (figure  from  [66])  and  cables 
(figure  from  [35]),  are  subject  to  VIV. 


1.2  Vortex-Induced  Vibrations  (VIV) 

As  a  cylindrical  object,  such  as  a  vibrissa,  moves  through  steady  flow,  it  can  develop 
significant  vibrations.  This  phenomenon,  known  as  vortex-induced,  vibrations  (VIV), 
has  been  studied  for  decades  and  several  reviews  have  been  written  [8,  114,  142], 
Despite  extensive  studies,  this  phenomena  continues  to  present  serious  challenges  in 
engineering  applications,  such  as  design  of  large-scale  ocean  structures  (Figure  1-22), 
which  will  experience  structural  fatigue  after  prolonged  vibrations. 


1.2.1  Why  does  VIV  occur? 

As  flow  goes  over  a  circular  cylinder,  vortices  are  shed.  These  exert  a  force  back 
onto  the  object.  If  the  natural  frequency  of  the  cylinder  is  sufficiently  close  to  the 
vortex  shedding  frequency  and  if  the  cylinder  damping  is  low  enough,  the  cylinder 
will  undergo  sustained  oscillations. 

Vortex  Shedding 

Bluff  bodies  are  geometries  that  have  a  cross-section  with  a  span-to-chord  ratio  of 
about  1.  A  circular  cylinder  is  the  canonical  example  bluff  body,  and  it  will  be  used 
throughout  this  chapter.  When  a  steady  flow  encounters  this  type  of  shape,  the  fluid 
is  forced  to  deviate  from  its  uniform  path  and  it  attempts  to  stay  attached  to  the 
cylinder  surface.  Reynolds  number  is  defined  in  equation  1.1,  where  U  is  the  free 
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stream  flow  velocity,  d  is  the  cylinder  diameter,  and  u  is  the  kinematic  viscosity  of 
water. 


Re  = 


Ud 

v 


(1.1) 


When  Re  is  >  40,  the  flow  separates  from  the  back  half  of  the  cylinder  and  peels 
off  into  swirling  vortices.  This  vortex  shedding  happens  in  an  alternating  pattern 
that  creates  a  turbulent  wake  known  as  the  von  Karman  vortex  street. 


This  street,  or  wake,  is  composed  of  periodically  shed  vortices  with  alternating 
rotation.  The  nondimensional  parameter  used  to  describe  this  flow  is  the  Strouhal 
number  (St)  and  is  defined  in  Equation  1.2  using  fs  as  the  vortex  shedding  frequency. 

St=f~§-  (L2) 

This  value  is  a  function  of  Reynolds  number  [11],  but  for  a  large  range  of  subcritical 
Reynolds  numbers,  it  remains  relatively  constant  at  ~  0.2  for  a  circular  cylinder. 


Fluid-Structure  Interaction 

Vortex  shedding  can  have  a  significant  effect  on  the  cylinder  itself.  If  the  cylinder  is 
either  flexible  or  flexibly  mounted,  its  natural  frequency  can  be  excited  by  the  vortex 
shedding  process.  When  these  two  frequencies  are  close  to  each  other,  the  cylinder 
may  undergo  large  amplitude  motions. 

This  is  because  a  vortex  is  a  low-pressure  region,  which  exerts  a  force  back  onto 
the  cylinder.  So  as  vortices  are  shed,  the  cylinder  experiences  an  oscillatory  transverse 
force.  This  transverse  force  corresponds  with  the  vortex  shedding  frequency,  which 
can  change  to  synchronize  with  the  motion  frequency.  This  is  called  lock-in.  The 
fluid  and  the  structure  operate  in  a  feedback  loop,  since  the  body  force  increases  the 
strength  of  the  vortices,  which  in  turn  increases  lift  and  drag  on  the  cylinder.  This 
results  in  high  oscillation  amplitudes. 
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Figure  1-23:  A  bluff  body  in  a  uniform  flow  can  be  simply  represented  as  a  1-degree-of- 
freedom  linear  oscillator.  The  free  stream  approaches  with  velocity  U  causing  vortex 
shedding  at  a  frequency  of  fs. 

1.2.2  Model 

An  elastically  mounted  bluff  body  in  a  free  stream  can  be  simply  represented  as  a 
1-D  harmonic  oscillator.  A  representative  mass- spring-damper  system  that  is  free  to 
move  in  the  transverse  (crossflow)  direction  with  position  y(t)  =  yo  sin  (cot)  is  shown 
in  Figure  1-23. 

The  cylinder  has  structural  mass  m,  structural  damping  b ,  and  spring  constant 
k.  The  motion  occurs  in  the  y  direction,  and  the  fluid  force  on  the  cylinder  in  this 
direction  is  given  by  L.  The  governing  equation  is  shown  in  Equation  1.3. 


my  +  by  +  ky  =  L  (1.3) 

This  can  be  expanded  into  Equation  1.4,  where  0  is  the  phase  angle  between 
motion  and  forcing  function,  and  u  =  2irf  is  the  frequency  at  which  the  body  is 
excited.  For  ease  of  analysis,  the  hydrodynamic  force  in  the  transverse  direction  is 
split  into  one  component  in  phase  with  velocity  and  one  in  phase  with  acceleration. 

my  +  by  +  ky  =  F0  sin  ( cot  f  o)  =  Fv  cos  (cot)  +  Fa  sin  (cot)  (1.4) 

This  free  vibration  system  has  mass  and  damping  that  are  dependent  on  the  fluid 
forces.  This  enables  an  effective  resonance  condition,  where  the  effective  natural 
frequency  varies  as  a  function  of  the  fluid  forcing.  In  VIV,  these  resonant  motions 
occur  over  a  much  wider  parametric  range  than  in  an  equivalent  mechanical  system, 
because  of  added  inertial  forces  from  vortex  shedding  and  fluid  acceleration. 
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The  specific  interplay  is  shown  using  the  following  equations  1.5  and  1.6  of  response 
amplitude  and  frequency.  These  can  be  derived  by  expanding  and  nondimensionaliz- 
ing  Equation  1.4  (as  shown  in  several  review  papers),  m*  is  the  mass  ratio,  Crn  is  the 
effective  added  mass  coefficient,  f*  =  f  /  fn  is  the  frequency  ratio,  ac  is  the  transverse 
force  coefficient  multiplied  by  a  constant,  and  (  is  the  damping  ratio.  The  two  key 
features  in  these  equations  are  the  damping  and  the  effective  added  mass. 

r  =  J  ,,r  77^  (1-5) 

V  m  +  Cm(/*) 


A* 


ac  sin  0 

+  Cm)C 


(1.6) 


From  equation  1.5,  it  is  seen  that  as  the  effective  added  mass  changes,  so  will  the 
excitation  frequency  (/).  The  feedback  between  the  fluid  forces  on  an  oscillating  bluff 
body  and  the  resulting  body  motion  can  cause  the  wake-capture  phenomenon.  This  is 
where  the  wake  frequency  fs  becomes  locked  into  the  structural  oscillation  frequency 
/.  The  effective  added  mass  is  thereby  altered,  bringing  the  natural  frequency  fn 
closer  to  the  wake  frequency  fs.  Also,  as  shown  in  equation  1.6,  for  low  mass-damping 
(m*£)  values,  the  amplitude  (A*)  is  greatly  affected  by  the  phase  angle  (0).  Even  if 
the  magnitude  of  the  forces  (ac)  are  low,  an  appropriate  phase  angle  will  make  the 
response  amplitude  large.  This  was  observed  by  [99]  for  low  mass-damping  cases  and 
is  shown  in  the  Griffin  plot  [142], 


Nondimensional  Parameters 

There  are  several  parameters  that  describe  the  dynamics  of  this  system.  Here  I 
pause  to  define  these  in  their  nondimensional  form  so  as  to  streamline  the  following 
discussions.  These  are  then  compiled  in  Table  1.1. 

Re : 

The  Reynolds  number  shows  the  relative  importance  of  inertia  forces  compared  to 
viscous  forces.  The  VI V  response  of  cylinders  varies  somewhat  as  this  value  changes 

[50]. 
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m*  and  m*(: 

The  mass  of  the  cylinder  can  be  normalized  with  respect  to  the  displaced  mass  of  the 
fluid  to  form  the  mass  ratio  (m*).  In  water,  this  value  is  often  low  since  the  density 
of  the  fluid  is  on  the  same  order  as  the  density  of  the  material.  The  damping  ratio 
of  the  system  is  written  as  £,  which  compares  the  structural  damping  to  the  critical 
damping.  The  complexity  of  the  VIV  phenomena  is  highlighted  in  the  fact  that  this 
damping  value  is  a  function  of  the  added  mass  The  effect  of  these  two  values  are 
often  combined  together  into  a  single  variable  (m*().  As  this  value  decreases,  the 
amplitude  response  increases  [50]. 

Ur  and  Urn : 

During  free-vibration  experiments,  the  flow  velocity  is  varied,  which  is  portrayed 
through  either  the  true  or  nominal  reduced  velocity  (Ur  and  Urn).  Similarly,  during 
forced-vibration  experiments,  the  motion  frequency  (/)  and  motion  amplitude  (A) 
are  varied,  which  affects  the  reduced  frequency  (fr).  This  is  equivalent  to  the  inverse 
of  Ur. 

Others: 

The  remaining  parameters  vary  as  a  function  of  the  reduced  velocity:  amplitude 
ratio,  frequency  ratio,  added  mass  ratio,  drag  coefficient,  and  lift  coefficient.  The 
analytical  solutions  to  these  are  derived  in  the  next  section.  Note  that  drag  (force 
in  the  direction  of  the  free  stream  flow)  and  lift  (force  perpendicular  to  the  free 
stream  flow)  are  written  in  dimensionless  form  by  dividing  by  the  dynamic  pressure: 
Pd  =  \pU2ds ,  where  s  is  the  cylinder  span. 

Derivations 

The  motion  profile  of  the  body  can  be  described  with  a  sinusoidal  profile,  as  can  its 
derivatives: 

y  =Asin(aA) 

y  =Au>  cos(cot)  (l-^) 

y  =—Aoj 2  sin  (cut) 
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name 

equation 

description 

reynolds  number 

Re  =  5^ 

inertial  force  /  viscous  force 

mass  ratio 

m*  —  m 

m  —  d2 

PTT^-S 

cylinder  mass  /  displaced  fluid  mass 

damping  ratio 

c-  b 

^  47rmfn 

system  damping  /  critical  damping 

nominal  reduced 
velocity 

u  —  — 

Urn  _  fnd 

flow  velocity  /  (natural  frequency*diameter) 

true  reduced  ve¬ 
locity 

u  —  — 

Ur  ~  fd 

flow  velocity  /  cylinder  velocity 

reduced  frequency 

f  fd 

Ar  —  U 

cylinder  velocity  /  flow  velocity 

amplitude  ratio 

motion  amplitude  /  diameter 

frequency  ratio 

f*  =  f 

In 

motion  frequency  /  natural  frequency 

added  mass  coeffi¬ 
cient 

m  a 

—  d2 
PTT^S 

effective  added  mass  /  displaced  fluid  mass 

drag  coefficient 

r* ,  —  D 

drag  force  /  (stagnation  pressure  *  projected 
area) 

“  §pU2ds 

lift  coefficient 

^  “  §pU2ds 

lift  force  /  (stagnation  pressure  *  projected  area) 

Table  1.1:  Nondimensional  parameters  relevant  to  VIV. 
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The  fluid  force  can  also  be  written  using  a  sinusoidal  profile  with  a  phase  shift 
(/>  between  the  force  and  the  motion.  This  can  be  expanded  in  terms  of  sines  and 
cosines,  as  in  Equation  1.8 


Fy(t)  =Lsm(u>t  +  <fi) 

=L  sin ((f))  cos(a jt)  +  L  cos (0)  sin(o jt) 


(1.8) 


The  equation  of  motion  (Equation  1.3)  can  be  expanded  by  plugging  in  the  profiles 
from  Equations  1.7  and  1.8.  Then,  two  different  equations  can  be  written  that  group 
terms  together  based  on  either  having  a  cos  (cut)  or  a  sin(t vt)  term.  These  are  shown 
in  Equation  1.9. 


b(Au)  cos(ujt))  =  Lsm(<f)  cos(cnt) 
(m  +  ma)(— Auj2  +  kA)(sm(ojt))  =  Lcos((f>)  sin(o;t) 


(1.9) 


The  amplitude  of  motion  and  the  fluid  forces  on  the  body  can  be  solved  by  ma¬ 
nipulating  these  equations.  The  method  to  do  so  is  simply  to  nondimensionalize  all 
the  terms  (using  the  definitions  in  Table  1.1)  and  isolate  the  desired  variable. 


Before  going  through  those  steps,  two  more  variables  should  be  defined  to  help 
simplify  the  equations.  The  lift  coefficient  in  phase  with  velocity  (Qjt))  and  the  lift 
coefficient  in  phase  with  acceleration  (Qj0)  are  the  dimensionless  versions  of  the  two 
components  of  lift.  The  lift  force  can  be  split  into  a  part  that  aligns  with  sin  (cut)  and 
one  that  aligns  with  cos(c at).  There  are  shown  in  nondimensional  form  in  Equation 
1.10. 


Q,v=  Ci  sin {(f>) 
Cl:a=  Cl  COS (0) 


(1.10) 


Now  all  the  components  are  in  place  to  solve  for  the  amplitude  of  motion  and  the 
fluid  forces  on  the  body. 
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A/d: 

Here  the  steps  to  do  so  are  shown  for  the  equation  that  shares  the  cos(c ut)  term 
(Equation  1.9),  which  leads  to  the  solution  of  A/d  given  in  Equation  1.11.  Again, 
this  is  done  by  nondimensionalizing  all  terms  and  isolating  the  desired  variable. 


(C47T  fnm)(Auj) 
(C47r/n(m>7r^s))(A(27r/)) 


A 

1 


Ci-pdsU 2  sin  (0) 
Ci,v  2  pds  ( Urn  fn  d) 

47 T3(j-)(m* 


(1.11) 


Similarly,  the  sin  (tut)  equation  can  be  rearranged  into  a  form  that  shows  the  fluid 
forces  on  the  cylinder,  given  by  the  added  mass  coefficient  (Equation  1.12). 


Cm  = 


-c,a 


jTYd2sAu!2 


(1.12) 


Ci,v: 

The  lift  coefficient  in  phase  with  velocity  is  a  particularly  important  variable  in  this 
context.  It  describes  the  energy  transfer  between  the  fluid  and  the  body.  The  lift  force 
is  also  assumed  to  have  a  sinusoidal  shape  with  a  phase  lag  (0).  The  nondimensional 
lift  is  then  split  into  two  terms:  one  in  phase  with  acceleration  and  one  in  phase  with 
velocity,  as  seen  in  Equation  1.13. 


Ci{t)  =  =  Ci  sin(ut  +  <p)  =  Ci  a  sin(a;t)  +  Ci  v  cos  (cot)  (1.13) 

Pd 

The  term  in  phase  with  acceleration  is  related  to  added  mass,  and  the  term  in 
phase  with  velocity  is  related  to  damping,  or  energy  transfer.  For  this  reason,  Ci>v 
can  be  used  to  determine  whether  energy  is  being  transferred  from  the  fluid  to  the 
model  (positive  sign,  which  indicates  VIV)  or  vice  versa. 
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Other  Vibrations 


Ct,v  = 


2  <  Q(t),y{t)  > 

T  y/<  WhW)  > 


(1.14) 


Vortex- induced  vibrations  fit  within  the  larger  class  of  flow- induced  vibrations.  In 
the  case  of  VIV,  the  vibrations  are  a  result  of  the  large-scale  vortex  shedding,  which 
causes  the  motions  to  be  stable  and  self-limiting  (to  about  A/d  =  1)  [131].  Another 
phenomenon  within  this  class  include  unstable  oscillations,  such  as  galloping,  which 
is  caused  by  geometrical  asymmetry.  These  induce  large  loads,  which  can  lead  to 
catastrophic  failures.  A  third  phenomenon  is  wake-induced  vibrations  (WIV),  where 
a  body  vibrates  in  response  to  an  external  vortex  wake.  In  that  case,  the  body  can 
still  undergo  a  lock-in  condition  with  wide  bandwidth  of  response.  The  vibrations  are 
not  self-limiting  in  that  case.  This  phenomenon  will  be  studied  later  in  the  thesis, 
when  the  response  of  the  whisker  to  an  externally  generated  wake  is  explored. 

VIV  oscillations  are  dominantly  in  the  crossflow  direction  (perpendicular  to  the 
free  stream).  However,  there  are  also  oscillations  that  arise  in  the  in-line  direction 
(parallel  to  the  free  stream),  which  are  usually  0.5  diameters  or  less.  These  in-line 
oscillations  occur  at  double  the  frequency  of  vortex  shedding,  since  one  pair  of  vortices 
is  shed  per  cycle  on  the  downstream  side  of  the  cylinder. 

1.2.3  VIV  Suppression  Mechanisms 

The  negative  impact  of  VIV  on  structures,  including  unsteady  forcing  as  well  as  drag 
amplification,  has  led  to  the  development  of  multiple  means  of  VIV  suppression,  as 
reviewed  in  [144,  133,  35]. 

Wake  Stabilizers 

One  method  of  VIV  reduction  is  to  stabilize  the  wake  by  placing  items  such  as  fairings, 
splitter  plates,  and  flags  (Figure  l-24(a-c))  into  the  wake.  These  streamline  the  body, 
which  serves  to  reduce  the  interaction  of  the  entrainment  layers.  Two  downsides  of 
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Figure  1-24:  Examples  of  VIV  suppression  devices:  (a)  fairing,  (b)  splitter  plate,  (c) 
flag,  and  (d)  a  cylinder  with  helical  strakes.  Figures  from  [35]. 


these  methods  are  that  they  require  extra  material  and  that  their  effectiveness  is 
dependent  on  a  unidirectional  flow. 


Strakes 

Another  means  of  reducing  VIV  is  to  change  the  surface  of  the  cylinder  itself.  The 
most  widely  used  suppression  device  in  the  offshore  oil  industry  is  that  of  helical 
strakes  (Figure  l-24(d)).  These  surface  protrusions  (with  height  up  to  25%  of  the 
diameter)  spiral  along  the  length  of  the  cylinder  in  order  to  disrupt  the  coherent 
formation  of  vortices. 

Such  surface  modifications  can  incur  the  expense  of  increased  drag.  In  particular, 
helical  strakes  experience  a  higher  drag  force  than  a  smooth  cylinder.  This  is  shown 
in  the  log-scale  plot  of  Figure  1-25. 


Bumpy  surfaces 

Bumps  on  various  objects  in  nature  may  serve  the  purpose  of  drag  reduction  [13].  The 
effectiveness  of  bumps  on  a  cylinder  to  reduce  VIV  was  tested  [99].  Hemispherical 
protrusions  were  attached  in  a  spiral  arrangement  to  a  cylinder,  as  shown  in  Figure 
l-26(a).  The  wake  patterns  were  seen  to  differ  along  the  span  (Figure  l-26(b)),  which 
allows  the  drag  on  the  model  to  be  reduced  up  to  25%  (based  on  the  diameter  of  the 
bare  cylinder). 

It  is  important  to  note  that  as  the  experiment  was  changed  from  fixed  boundary 
conditions  to  an  increasingly  flexible  mounting  (indicated  by  a  lower  mass-damping 
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Figure  1-25:  The  drag  coefficient  of  various  VIV  suppression  devices  is  higher  than 
that  of  a  smooth  cylinder  at  certain  Reynolds  numbers.  This  is  especially  true  for  a 
cylinder  with  helical  strakes.  Figure  from  [35]. 


Figure  1-26:  (a)  A  cylinder  with  bumpy  protrusions  placed  at  45°  angular  separation 
and  with  pitch  of  7  diameters.  Figure  from  [100].  (b)  The  wake  behind  this  model 
ranges  from  wide  to  narrow,  depending  on  whether  the  model  has  a  peak  or  trough 
in  the  disturbance  there.  Figure  from  [99]. 
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Figure  1-27:  Maximum  transverse  vibration  amplitude  of  different  geometries.  As  the 
mass-damping  {rn*Q  decreases,  all  geometries  vibrate  with  higher  amplitude,  solid 
circle :  smooth  circular  cylinder;  plus  sign :  sinuous  cylinder;  others :  bumpy  cylinders 
with  various  bump  sizes).  Figure  from  [99]. 

value) ,  the  transverse  response  of  the  model  increased  substantially.  This  is  the  same 
trend  as  seen  on  a  smooth  cylinder  (Figure  1-27),  indicating  that  suppression  of  vortex 
shedding  under  fixed  conditions  does  not  necessarily  equate  to  reduction  of  VIV. 

Wavy  surfaces 

A  final  method  is  seen  in  the  use  of  wavy  surfaces  on  cylinders.  Periodic  variations 
along  the  span  of  a  cylinder  decorrelate  vortex  shedding.  Studies  by  [2,  75,  81,  80, 
83,  145]  have  shown  cylinders  of  varying  waviness  effective  in  reducing  the  mean  drag 
force  by  about  20%  and  the  unsteady  lift  force  significantly,  under  certain  conditions. 
The  effectiveness  is  attributed  to  the  formation  of  an  undulating  wake  which  de¬ 
correlates  the  induced  forces,  as  well  as  the  formation  of  streamwise  vortices.  [83] 
found  a  maximum  effectiveness  when  the  wavelength  to  mean  diameter  ratio,  A / Dm, 
is  close  to  2.  [80]  found  a  second  optimal  value  of  A /Dm  =  6,  which  provides  even 
better  performance.  The  latter  value  agrees  with  the  reported  optimal  wavelength  to 
diameter  ratio  for  strakes,  which  is  about  5  for  single  strakes,  and  between  15  and  17 
for  triple  start  strakes,  i.e.,  effectively,  between  5  and  5.6. 
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whisker  "noise" 


Figure  1-28:  Schematic  of  expected  response  of  the  whiskers  as  the  seal  is  swimming: 
the  whiskers  undergo  significant  VIV,  making  for  a  “noisy”  signal. 

1.3  Synthesis  of  Sensing  and  Vibrations 

These  two  topics  of  hydrodynamic  perception  and  flow-induced  vibrations  are  inves¬ 
tigated  together  in  the  remainder  of  this  thesis. 

The  reason  for  their  synthesis  is  triggered  by  the  following  (depicted  in  Figure 
1-28):  Given  that  vortex- induced  vibrations  can  result  in  significant  amplitudes  of 
motion  on  cylindrical  objects,  we  expect  that  the  seals’  whiskers  should  vibrate  sig¬ 
nificantly  as  the  seals  are  swimming.  If  the  whiskers  are  sensors,  shouldn’t  they 
experience  strong  background  noise  from  VIV? 

1.3.1  Previous  Work 

The  study  by  [58]  was  the  first  to  attempt  to  synthesize  these  two  topics.  The  authors 
made  a  claim  that  the  whisker  geometry  was  able  to  suppress  VIV.  They  postulated 
that  this  would  enable  heightened  sensing  capabilities. 

In  that  study,  three  tests  were  conducted.  First,  a  camera  was  mounted  on  a 
harbor  seal’s  head  to  gain  a  close  view  of  the  whiskers.  It  was  observed  that  the 
whiskers  did  not  vibrate  (at  least  beyond  the  0.2 mm  resolution  of  the  camera)  as  the 
seal  swam  normally.  Second,  whisker  specimens  were  placed  in  a  flume.  Forces  at  the 
base  of  the  undulatory  harbor  seal  whisker  were  measured  to  be  6.2  times  lower  than 
that  of  a  sea  lion  whisker,  which  has  a  smooth  profile.  Finally,  CFD  was  conducted  on 
the  numerical  representation  of  the  harbor  seal  whisker.  The  whisker  was  held  rigidly 
and  forces  were  calculated.  Crossflow  forces  were  seen  to  reduce  90%  and  inline  forces 
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Figure  1-29:  Schematic  of  the  hypothesis  behind  this  thesis.  In  open  flow,  the  whiskers 
produce  low  VIV.  This  enables  the  seal  to  more  easily  detect  unsteady  flow  features, 
which  induce  large  response. 


40%  compared  with  a  smooth  circular  cylinder. 

Their  result  suggests  that  the  whiskers  can  move  smoothly  in  a  uniform  flow,  while 
experiencing  few  vibrations  and  therefore  accurately  measure  characteristics  of  the 
oncoming  flow. 


Other  whisker  vibration  studies 

Since  then,  a  few  other  parties  have  carried  out  studies  on  the  vibration  properties 
of  seal  whiskers  on  both  real  and  model  specimens  [91,  143,  95].  The  important 
thing  to  note  is  that  none  of  the  experiments  that  provided  quantitative  results  were 
conducted  under  conditions  that  allowed  the  whisker  to  vibrate  freely  in  response  to 
the  flow.  The  experiments  described  in  Chapter  2  are  designed  to  address  this  gap. 


1.3.2  Hypothesis 

The  evidence  of  reduced  fluid  forcing  felt  by  whiskers  in  open  water  [58]  leads  to  the 
hypothesis  of  this  thesis  (depicted  in  Figure  1-29):  The  whiskers  experience  low 
VIV  in  open  flow  so  that  the  seal  can  more  easily  detect  unsteady  flow 
features. 
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1.4  Outline  of  Current  Study 


This  thesis  investigates  whether  or  not  the  harbor  seal  whisker  morphology  can,  at 
least  in  part,  be  responsible  for  the  heightened  sensory  ability  of  the  harbor  seals.  The 
fluid  dynamics  of  the  whiskers’  response  to  different  flows  is  investigated,  in  order  to 
uncover  mechanisms  by  which  vibrational  response  could  aid  sensing. 

Chapter  1  provided  background  information  on  flow  sensing  in  marine  animals. 
Also,  the  basics  of  flow- induced  vibrations  were  introduced  along  with  the  relevant 
nomenclature  and  parameters. 

Chapter  2  quantifies  the  harbor  seal  whisker  geometry’s  VIV  properties.  Experi¬ 
ments  are  conducted  using  forced  vibrations  on  a  scaled  whisker  model,  allowing  for 
dynamic  investigation  of  the  vibrational  response.  Explanations  for  the  whisker’s  re¬ 
sponse  are  provided  through  wake  visualizations  and  measurement  of  the  added  mass 
and  damping  properties.  The  effect  of  the  angle  with  which  the  geometry  faces  the 
flow  is  also  studied. 

Chapter  3  describes  the  design  and  characterization  of  a  flow  sensor,  which  incor¬ 
porates  the  unique  vibration  properties  discovered  in  the  harbor  seal  whisker  geometry 
into  a  free  vibration  configuration. 

Chapter  4  presents  experiments  of  the  freely  vibrating  whisker  model  encounter¬ 
ing  a  vortex  wake.  This  simple  representation  of  the  seal  interacting  with  a  fish  wake 
is  used  to  elucidate  potential  mechanisms  employed  for  wake  detection.  Strain  mea¬ 
surements  and  dye  visualizations  reveal  the  whiskers’  unique  mode  of  interaction  with 
external  vortices. 

Chapter  5  concludes  the  thesis  by  summarizing  the  principal  contributions.  Sug¬ 
gestions  for  future  work  and  further  hypotheses  of  the  real  seals’  behavior  are  pro¬ 
vided. 
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Chapter  2 


Whisker  in  Open  Flow:  Forced 
Vibrations 


Can  harbor  seal  whiskers  passively  reduce  VIV? 

How  sensitive  is  the  vibration  to  changes  in  flow  angle? 

2.1  Motivation 

In  contrast  to  other  animals  with  whiskers,  most  phocid  seals  exhibit  a  unique  whisker 
geometry  consisting  of  undulating  radii  in  both  major  and  minor  axes  and  an  elliptical 
cross  section  (average  aspect  ratio  =  2).  A  close-up  view  of  the  whiskers  in  situ 
reveals  this  “beaded”  pattern  (Figure  2-1).  Observations  of  the  behavior  of  harbor 
seals  suggest  that  this  geometry  enables  enhanced  underwater  sensing  capabilities. 
These  aquatic  animals  are  known  to  employ  these  sensing  techniques  underwater 
while  swimming  with  relatively  high  forward  motion,  hence  the  sensing  mechanism 
needs  to  be  sensitive  enough  to  detect  minute  changes  in  the  local  flow  field,  while 
robust  enough  to  noise  induced  through  self-excitation  from  the  animal’s  forward 
motion. 

Few  studies  have  investigated  the  hydrodynamic  interaction  of  the  whisker  geom¬ 
etry  with  its  wake  and  hydrodynamic  scaling  effects  on  the  geometry,  particularly  in 
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Figure  2-1:  Close-up  view  of  the  whiskers  on  a  harbor  seal’s  muzzle.  The  “beaded” 
pattern  indicates  an  undulatory  geometry.  Photo  taken  courtesy  of  the  New  England 
Aquarium. 

explaining  the  effective  use  of  this  geometry  as  a  sensor  or  as  a  vortex  suppression 
device. 

2.1.1  Whisker  in  Frontal  Flows  (ct  =  0°) 

[58]  posited  that  the  unique  geometry  of  the  harbor  seal  whisker  suppresses  the  natural 
flow-induced  vibrations  experienced  by  the  whiskers  because  of  the  forward  swimming 
motion  of  the  animal.  This  statement  was  based  on  evidence  that  vortex  shedding  is 
suppressed  in  the  wake  of  the  whiskers  when  aligned  with  the  direction  of  swimming 
for  the  seal,  which  is  what  is  referred  to  here  as  “frontal  flows” .  This  study,  however, 
was  limited  to  a  static  investigation  of  the  forces  exerted  on  the  whisker  geometry.  For 
systems  near  a  resonant  condition,  a  reduction  in  the  magnitude  of  dynamic  forces 
acting  on  a  fixed  body  does  not  preclude  a  reduction  in  the  dynamic  response  of  the 
body  (as  seen  in  [99]).  A  body’s  dynamic  response  is  dependent  on  damping  and 
effective  added  mass,  which  influence  the  phase  between  the  body  motion  and  the 
fluid  forces,  potentially  leading  to  lock-in  over  a  wide  parametric  range. 

This  chapter  examines  the  full  VIV  characteristics  of  the  whisker  geometry.  A 
scaled  whisker  model  (Figure  2-2)  is  tested  under  dynamic  conditions  and  a  range  of 
mass-damping  values.  The  goal  is  to  determine  if  the  harbor  seal  whisker  geometry 
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Figure  2-2:  (Top)  Side  view  and  (Bottom)  front  view  of  a  fabricated  replica  of  the 
whisker  test  section,  capturing  the  asymmetric  undulatory  geometry.  Stereolithogra¬ 
phy  was  used  to  create  this  plastic  hollow  structure.  Dimensions  were  taken  directly 
from  [58]  and  enlarged  by  a  factor  of  50. 

is  able  to  suppress  VI V  under  conditions  that  allow  for  self-excited  motions.  If  it  is 
able,  it  would  imply  that  the  geometry  is  specialized  to  reduce  the  background  noise 
caused  by  the  whiskers  own  wake  in  the  detection  of  the  upstream  target  a  desirable 
feature  that  would  motivate  a  seal  whisker-inspired  flow  sensor. 

Experimental  force  measurements  are  employed  over  a  range  of  oscillations  to 
find  the  parametric  regime  in  which  VIV  occurs.  Flow  visualizations  are  conducted 
to  elucidate  the  wake  patterns  associated  with  these  cases.  The  geometry  is  first 
tested  in  the  streamlined  direction  (a  =  0°). 

2.1.2  Whisker  at  Angles  of  Attack  (a  >  0°) 

We  lack  solid  description  of  the  whisker  array’s  position  during  prey  pursuit  in  the 
wild,  however  seals  often  keep  their  whiskers  erect  during  tactile  sensing  [27]  and  may 
do  so  as  well  during  hydrodynamic  tracking  [18].  This  could  serve  to  cover  a  wide 
area  of  the  environment.  Given  that  this  may  happen  and  that  the  whiskers  may  sit 
in  the  cheek  at  different  angles  within  the  seals’  cheek,  it  is  desired  to  know  what 
changes  the  vibrational  response  may  have  as  the  whiskers  see  the  flow  from  different 
orientations. 

The  same  methods  of  experimental  force  measurement  and  flow  visualizations  are 
used  on  the  whisker  geometry  undergoing  forced  vibrations  while  placed  in  different 
orientations  to  the  flow. 
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Figure  2-3:  Seals’s  whiskers  may  face  different  orientations.  Photo  taken  at  the  New 
England  Aquarium. 


2.2  Methodologies 


2.2.1  Force  Measurements:  Forced  Vibrations 

Forced  vibration  of  a  cylinder  has  been  shown  to  accurately  predict  the  fluid  forces 
felt  on  a  freely  vibrating  cylinder  [64,  94,  124,  93].  This  is  the  methodology  used  to 
investigate  the  whisker  VI V  properties. 


Apparatus 

The  geometry  provided  in  [58]  was  scaled  up  50  times  and  fabricated  into  a  model 
using  stereolithography.  The  model  was  vertically  submerged  in  a  towing  tank.  It 
was  mounted  to  a  force  transducer  at  its  base,  which  is  above  the  water.  A  simple 
harmonic  motion  profile,  given  in  Equation  2.1,  was  imposed  on  the  model  in  the 
crossflow  direction  (y)  while  it  was  towed  (at  velocity  U )  down  the  length  of  the  water 
tank.  The  expected  vibration  was  thus  simulated.  Figure  2-4  shows  a  schematic  of 
this  setup. 


y(t )  =  Asin(ciA)  (2.1) 

The  same  process  was  repeated  for  multiple  combinations  of  amplitudes  ( A )  and 
frequencies  (/  =  oj/('2tt)).  These  are  written  in  dimensionless  form  as  amplitude 
ratio  ( A/d *)  and  reduced  frequency  (/r  =  fd*/U).  The  lift  coefficient,  where  lift  is 
the  hydrodynamic  force  in  the  direction  of  crossflow  motion  (L  =  Fy),  was  recorded 
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Figure  2-4:  Schematic  of  the  forced  vibration  setup.  The  model  hangs  vertically  in 
the  water  tank.  It  is  directly  attached  to  a  force  sensor  and  a  linear  stage  that  drives  a 
sinusoidal  motion  (an  example  trajectory  is  depicted  in  red)  while  the  carriage  moves 
the  setup  down  the  length  of  the  tank  at  velocity  U. 


simultaneously  with  the  position  of  the  model.  The  phase  between  these  two  signals 
was  determined  through  the  calculation  of  C^v,  the  lift  coefficient  in  phase  with 
velocity. 

The  sign  of  this  variable  indicates  if  the  model  underwent  VI V  or  not.  If  the  value 
is  positive,  it  indicates  an  energy  transfer  from  the  fluid  to  the  body,  which  is  what 
occurs  during  VIV.  If  the  variable  is  negative,  the  energy  transfer  goes  from  the  body 
to  the  fluid,  so  that  case  is  not  physical.  So  those  runs  that  resulted  in  a  value  of 
C\v  >  0  are  cases  that  could  occur  if  the  body  was  allowed  to  freely  vibrate.  For 
more  background  on  the  CitV  variable,  refer  back  to  Section  1.2.2. 

This  forced  vibration  technique  is  preferable  to  the  traditional  free  vibration  ex¬ 
periment  in  a  few  ways.  One  is  the  ease  of  experimental  setup.  A  single  linear  stage 
is  needed  to  run  the  present  technique,  rather  than  a  full  spring-dashpot  system  that 
is  correctly  tuned.  Also,  with  free  vibrations,  a  minor  change  in  flow  speed  can  cor¬ 
respond  to  a  large  difference  in  amplitude  and  a  change  in  vortex  shedding  mode. 
These  transitions  are  challenging  to  study  using  free  vibration,  since  the  conditions 
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cannot  be  closely  controlled.  Forced  vibrations  allow  for  parameters  to  be  specified 
with  greater  precision  and  to  be  covered  over  a  wider  range.  However,  the  number  of 
experiments  required  to  cover  the  appropriate  parameter  space  is  larger  for  forced  vi¬ 
brations.  There  are  two  variables  to  be  varied  (frequency  and  amplitude),  as  opposed 
to  one  (tow  velocity),  when  using  free  vibration. 


Parameters 

A  description  of  the  experimental  parameters  first  start  from  a  reference  cross  section 
of  the  whisker  geometry.  The  reference  (Equation  2.2)  is  taken  as  an  ellipse  with  in¬ 
line  diameter  (/)  and  crossflow  diameter  (d)  found  from  the  average  of  the  maximum 
and  minimum  diameters  from  the  two  extreme  cross  sections.  The  values  shown  here 
correspond  to  the  size  of  the  model  used  in  the  first  set  of  experiments  (50x-scale 
model) . 


d  ( dmax  T  j 2  2.65cm 


l  (.^niax  +  ^min)/2  5.35  CTTl 


(2.2) 


The  model  was  tested  at  different  angles  of  attack  (a)  with  respect  to  the  oncoming 
flow  direction.  This  angle  is  depicted  in  Figure  2-5.  First,  tests  were  conducted  at 
a  =  0°,  such  that  the  orientation  is  streamlined.  Then,  different  values  of  a  up  to 
90°  were  tested. 

The  reference  diameter  used  in  all  cases  is  the  projected  diameter  ( d *).  This  is 
defined  as  the  length  of  the  cross  section  (the  reference  cross  section  given  in  Equation 
2.2)  that  is  perpendicular  to  the  oncoming  flow  direction.  This  was  calculated  for  each 
angle  using  Equation  2.3.  That  can  be  derived  by  first  considering  the  tangent  line 
to  the  ellipse  that  passes  through  the  furthest  point  of  the  ellipse  on  the  x-axis,  and 
considering  the  line  perpendicular  to  that,  which  also  runs  through  the  origin.  Then 
the  distance  from  the  origin  to  the  intersection  of  those  two  perpendicular  lines  can 
be  calculated. 
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Figure  2-5:  Orientation  of  a  cross  section  of  the  model.  The  solid  line  cross  section 
represents  the  a  =  0°  case  and  the  dotted  line  cross  section  represents  the  a  >  0° 
cases.  The  projected  diameter  (d*)  is  the  dimension  used  for  all  cases  and  is  the 
length  perpendicular  to  the  flow,  which  approaches  with  velocity  U . 


d*  =  2\fb2  sin2  ct  +  a2  cos2  a  (2.3) 

Forced  vibrations  were  employed  here  on  three  different  geometries:  the  whisker 
model  pictured  above  (uses  the  geometry  as  described  in  [58]),  an  elliptical  cylinder, 
and  a  circular  cylinder.  All  have  the  same  material  properties  and  the  same  average 
cross-stream  diameter  d  =  2.65 cm  (at  a  =  0°),  which  is  based  on  a  50x  scale  of  the 
real  whisker  size.  For  the  low  angle  of  attack  cases  (a  <  45°),  a  towing  speed  of 
U  =  0.2m/ s  was  used,  which  corresponds  to  a  Reynolds  number  (Re  =  Ud*/u )  of 
5,300  —  7,030.  For  the  tests  at  a  >=  45°,  smaller  scale  models  were  used  in  order 
to  fit  the  vibration  amplitude  limits  of  the  experimental  setup.  This  model  is  a  30x- 
scale  version  of  the  real  whisker,  giving  it  d  =  1.59 cm  (at  a  =  0°).  A  towing  speed 
of  U  —  0.1  m/s  was  used,  resulting  in  Re  =  2,530  —  3,210.  The  nondimensional 
length  ( L/d *)  was  the  same  in  all  cases,  in  order  to  capture  the  same  number  of 
undulation  wavelengths.  The  mass  ratio,  which  is  found  using  the  model  volume 
given  in  Solidworks,  is  similar  for  all  of  the  models,  ranging  from  m*  =  0.5  to  1.1. 
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Calculations 


Forces  were  measured  on  the  models  while  they  underwent  imposed  transverse  (y) 
oscillations  of  a  certain  frequency  (written  in  nondimensional  form  as  reduced  fre¬ 
quency:  fr  =  fd*/U)  and  amplitude  (written  in  nondimensional  form  as  amplitude 
ratio:  A/d *)  and  were  towed  down  the  tank  at  U  =  0.2  m/s.  This  was  repeated 
for  different  combinations  of  frequency  and  amplitude.  These  combinations  differed 
for  the  circular  cylinder  versus  the  elliptical  cylinder  and  whisker,  as  their  excitation 
ranges  were  found  to  differ.  In  the  case  of  the  whisker,  each  of  these  runs  was  re¬ 
peated  10  times  and  averaged.  Signals  were  acquired  with  a  ATI  Gamma  6-axis  force 
sensor  (a  ATI  Nano  6-axis  sensor  was  used  for  low- amplitude  cases  on  whisker)  and 
collected  by  a  NI  USB-6218  DAQ  board  at  a  sampling  rate  of  1000Hz.  The  raw  data 
was  processed  through  a  two-pole,  low-pass  Butterworth  filter.  Inertial  forces  were 
calculated  by  measuring  the  model’s  mass  and  acceleration,  and  they  were  removed 
from  the  total  force  signal. 

The  remaining  hydrodynamic  force  component  was  used  to  calculate  four  variables 
that  describe  the  vibrational  response:  the  drag  coefficient  C-h,  the  lift  coefficient  in 
phase  with  velocity  Clv,  the  effective  added  mass  coefficient  Cm ,  and  the  damping 
ratio  (• 


The  drag  coefficient  was  calculated  from  the  in-line  force  ( D )  during  the  forced 
vibration  experiments.  It  is  normalized  with  respect  to  the  projected  diameter  ( d *): 


Equation  2.4. 


(2.4) 


The  lift  coefficient  in  phase  with  velocity,  shown  in  Equation  2.5,  depicts  the 
structure’s  excitation,  as  the  body  response  amplitude  is  proportional  to  the  force 
component  that  is  in  phase  with  velocity  (as  seen  in  Equation  1.6).  A  positive  Clv 
value  indicates  momentum  transfer  from  the  fluid  to  the  body  and  therefore  existence 
of  VIV.  Equivalently,  a  negative  value  indicates  transfer  in  the  opposite  direction  and 
no  VIV. 
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Clv  =  CL  sin  (<p) 


(2.5) 


Vortex  shedding  and  acceleration  of  the  surrounding  fluid  exert  forces  on  the  body 
that  are  in  phase  with  the  body  acceleration.  Together,  these  are  described  by  the 
effective  added  mass  coefficient.  The  effective  added  mass  is  determined  in  Equation 
2.6  by  first  dividing  the  lift  force  in  phase  with  acceleration  ( La  =  L  cos(d))  by  the 
magnitude  of  acceleration  ( y  term  from  Equation  1.4).  This  is  then  nondimensional- 
ized  with  respect  to  displaced  fluid  mass  in  Equation  2.7,  where  V  is  volume. 


ma  _  CLa 
pV  -  2r>(£)/? 


(2.7) 


(  was  found  using  the  same  method  used  to  derive  Cm ,  but  now  equating  the 
magnitude  of  the  damping  term  ( y  term  in  Equation  1.4)  to  the  lift  in  phase  with 
velocity  ( Lv  =  Lsin(0)),  as  in  Equation  2.8.  The  system  damping  ( b )  is  then  nondi- 
mensionalized  with  respect  to  critical  damping  ( bcr )  in  Equation  2.9.  Mass  ratio  (rn* ) 
is  used  to  describe  the  model  mass  with  respect  to  the  displaced  fluid  mass. 


Lv  Lv 

lull  Auj 


(2.8) 


2.2.2  Particle  Image  Velocimetry  (PIV):  Steady  Flow 

In  order  to  first  verify  the  results  shown  in  [58]  for  the  static  case,  flow  visualizations 
were  conducted  using  Particle  Image  Velocimetry  (PIV). 
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Figure  2-6:  Schematic  of  the  PIV  setup.  The  model  again  hangs  vertically  in  the 
water  tank.  A  laser  sheet  illuminates  a  thin  cross-section  of  the  model  and  the  fluid 
in  that  plane.  A  camera  mounted  below  the  tank  captures  the  motion  of  neutrally 
buoyant  particles  in  the  water. 

Apparatus 

The  model  was  mounted  as  before  (schematic  shown  in  Figure  2-6),  but  now  in  the 
towing  tank  at  our  CENSAM  facility  in  Singapore,  and  it  was  towed  down  the  length 
at  U  =  0.1  m/s.  A  high-speed  camera  (Phantom  V10)  operating  at  100  frames/sec 
remained  stationary  below  the  tank.  A  continuous  wave  laser  (A  =  532 nm)  located 
behind  the  model  illuminated  the  particle-laden  flow  with  a  horizontal  sheet  of  light 
at  40%  span.  Cross-correlation  of  the  captured  images  was  used  to  determine  the 
motion  of  the  particles  within  a  128xl28-pixel  sized  interrogation  window  followed 
by  a  64x64-sized  window,  both  with  50%  overlap.  Light  smoothing  was  carried  out 
through  use  of  a  3x3  filter.  The  same  thing  was  repeated  for  an  elliptical  cylinder  of 
the  same  average  diameter  as  the  whisker  in  both  directions  (2:1  ratio). 

2.2.3  Particle  Image  Velocimetry  (PIV):  Forced  Vibrations 

The  same  setup  and  processing  steps  used  in  the  steady  flow  PIV  experiments  was 
used  again,  but  now  with  imposed  motion  profiles  to  view  the  wake  under  dynamic 
conditions.  Based  on  the  force  measurements,  the  runs  corresponding  to  the  maxi¬ 
mum  VIV  amplitude  were  chosen  and  used  for  visualization  here.  For  the  whisker 
model,  this  corresponds  to  A/d*  =  0.075,  fd*/U  =  0.33,  and  for  the  elliptical  cylinder, 
this  corresponds  to  A/d*  =  0.175,  fd*/U  =  0.27. 
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Figure  2-7:  Cd  is  plotted  as  a  function  of  amplitude  ratio  and  reduced  frequency 
for  the  circular  cylinder,  elliptical  cylinder,  and  whisker  (from  left  to  right).  Re  = 
5,300.  Asterisks  indicate  where  experiments  were  performed.  The  contour  levels  are 
in  increments  of  0.1  up  to  1,  then  in  increments  of  0.25  to  1.5,  and  in  increments  of 
0.5  beyond  that. 


2.3  Results  (ce  =  0°) 

2.3.1  Force  Measurements 

Drag  ( Cd ) 

Figure  2-7  shows  a  comparison  of  the  models’  responses  during  the  forced  vibration 
tests.  The  drag  coefficient  goes  from  1.25  to  1  to  0.7  for  the  circular  cylinder,  el¬ 
liptical  cylinder,  and  whisker  in  the  near-zero  amplitude  case.  This  shows  that  the 
whisker  geometry  reduces  the  mean  drag  one  of  the  two  desired  components  for  V1V 
suppression. 

The  drag  is  defined  as  the  hydrodynamic  force  acting  on  the  body  in  the  in-line 
direction  ( D ).  The  formula  was  given  in  nondimensional  form  in  Equation  2.4,  where 
s  is  the  submerged  span. 

If  one  compares  the  drag  during  a  forced  vibration  (crossflow)  of  the  same  reduced 
frequency  and  amplitude  ratio  (for  example,  fr  =  0.2  and  A/d*  =  0.4),  the  drag  is  seen 
to  decrease  from  circular  cylinder  to  elliptical  cylinder  to  whisker.  In  that  example, 
it  goes  from  1.9  to  1.3  to  1.1.  Comparison  of  drag  during  zero  vibration  gives  values 
of  1.5  to  1  to  0.8,  in  the  same  order. 
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Figure  2-8:  Clv  is  plotted  as  a  function  of  amplitude  ratio  and  reduced  frequency. 
The  region  of  positive  values  indicates  where  VI V  occurs.  The  size  of  this  region 
decreases  from  the  circular  cylinder  to  the  elliptical  cylinder  to  the  whisker  (from  left 
to  right).  The  zero  contour  is  depicted  in  bold.  Re  =  5,300.  Asterisks  indicate  where 
experiments  were  performed.  The  contour  levels  are  in  increments  of  0.5  down  to 
-0.5,  then  in  increments  of  2  from  -2  on,  and  in  increments  of  -6  beyond  that. 


Energy  Transfer 

The  second  desired  component  is  the  reduction  in  vibration  response  (amplitude  and 
bandwidth).  This  is  calculated  via  C^v ,  which  is  plotted  in  contours  against  amplitude 
ratio  ( A/d *)  and  reduced  frequency  (/r)  in  Figure  2-8.  Presented  in  this  manner,  the 
full  VIV  characteristics,  including  potential  self-excited  motions,  are  captured.  The 
asterisks  indicate  the  frequency  and  amplitude  values  of  the  multiple  runs  conducted. 
In  the  whisker  case,  the  Clv  values  displayed  are  the  average  of  10  trials  for  each 
asterisk.  The  area  of  positive  value  (the  zero  contour  peaks  at  A/d*  ~  0.075)  for  the 
whisker  is  drastically  smaller  than  that  of  a  circular  cylinder.  Thus,  the  regime  in 
which  VIV  is  experienced  by  the  whisker  is  largely  suppressed.  The  baseline  case  of 
a  cylinder  with  the  same  cross-stream  diameter  has  peak  vibration  at  A/d*  ~  0.7, 
which  aligns  well  with  the  value  found  by  [49].  The  elliptical  cross  section  plays  a 
significant  role  in  reducing  the  vibrations,  since  the  elliptical  cylinder  exhibits  a  peak 
vibration  of  A/d*  ~  0.18. 


Variation  in  CLtV 

Another  interesting  point  is  that  there  is  variance  in  certain  Clv  values  across  the  10 
trials.  While  some  runs  appear  with  negative  Clv  on  average,  they  still  experience 
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a  positive  value  some  of  the  time.  The  probability  that  those  runs  will  experience 
VIV  based  on  the  ten  trials,  is  shown  in  Figure  2-9  as  a  third  dimension  on  top  of 
the  average  contour  plot.  It  is  seen  that  values  of  Clv  fluctuate  in  time,  so  while 
they  may  show  up  negative  on  average,  the  actual  value  is  positive  once  lock-in  has 
been  reached.  The  whisker,  therefore,  is  shown  to  require  a  longer  time  to  reach 
lock-in  than  a  cylinder.  The  length  of  tank  used  limited  the  length  of  time  available 
to  capture  data.  This  variability  also  suggests  that  a  greater  variation  in  shedding 
pattern  along  the  whisker  span  versus  the  cylinder. 

An  explanation  for  this  variation  can  be  unveiled  if  we  look  at  the  Cl,v  value  versus 
time.  The  value  actually  fluctuates  over  time  (refer  to  subplots  in  Figure  2-9).  For 
these  low-amplitude  cases,  lock-in  is  not  reached  until  multiple  cycles  of  oscillation 
have  occurred.  The  scatter  plots  shown  for  each  cycle  depict  the  spread  of  Cl,v  values 
over  10  trials.  The  middle  black  line  with  dots  shows  the  ensemble  average  for  each 
cycle.  With  the  exception  of  Run  c,  these  cases  have  a  decreasing  standard  deviation 
of  Cl,v  over  time,  and  these  values  eventually  settle  into  a  region  that  is  consistently 
either  positive  or  negative.  Run  c  does  not  appear  to  have  settled  into  a  constant 
value,  but  the  sign  is  not  in  question  it  is  clearly  positive  after  the  initial  phase. 

This  extended  time  needed  for  the  whisker’s  wake  to  synchronize  (compared  to 
a  cylinder)  means  that  only  the  latter  portion  of  data  should  be  considered  in  the 
analysis.  The  initial  start-up  phase  can  be  discarded  from  the  calculation  of  average 
Cl,v-  R  is  chosen  here  to  only  consider  the  segment  of  data  where  the  standard 
deviation  has  reached  half  its  initial  value.  Using  that  shortened  data  set,  the  average 
of  the  ensemble  averages  is  then  calculated  and  is  stated  as  the  representative  Cl,v 
value  for  that  point  on  the  plot.  This  means  that  Runs  b,  c,  e,  and  /  can  be  considered 
always  positive  and  Run  a  can  be  considered  negative.  Although  Run  c  does  not  have 
enough  time  to  settle  in,  it  clearly  reaches  a  positive  value. 


85 


2 

3 

o 

£ 


jo 

3h 

CD 

£ 

O 

3 

'o 

a 

cd 


o3 


X) 

o3 

X> 

O 

Dh 

Oh 

CD 

CD 

-+-? 

cd 

#CD 

3 

3 

m 

f-H 

cd 

X) 

3 

3 

3 

f-H 

CD 

> 

CD 

H 


jo 

3h 

c5 

cd 

CO 

3 

£ 

CD 


£ 

CD 

CD 

CO 

£ 

O 


.1  ^ 

'3  CD 

3  > 

•  i-H 

ci  o 

CN  ^ 
CD 

S  £ 

fcuO  O 

fa  'm 


86 


The  solid  black  dots  indicate  the  average  (over  10  trials)  CitV  value  at  each  cycle  of  oscillation.  The  scatter  plots  (red)  show  the 
spread  of  the  10  individual  values  at  each  cycle.  The  standard  deviation  decreases  over  time,  indicating  that  lock-in  is  reached 
after  some  initial  start-up  phase. 
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Figure  2-10:  Cm  contours  plotted  as  a  function  of  amplitude  ratio  and  reduced  fre¬ 
quency  for  the  circular  cylinder,  elliptical  cylinder,  and  whisker  (from  left  to  right). 
The  zero  Clv  contour  is  copied  over  and  shown  in  black.  The  contour  levels  are  in 
increments  of  1  and  the  3.5  contour  is  included  as  well.  Asterisks  indicate  where 
experiments  were  performed. 


Effective  Added  Mass  ( Cm ) 


Figure  2-10  shows  the  effective  added  mass  coefficient  for  the  circular  cylinder  to 
range  from  slightly  negative  to  3  within  the  excitation  region  (as  denoted  by  the 
zero  contour  of  Cl>v),  which  matches  the  expected  values  in  literature  [49].  Those 
values  deviate  from  the  potential  flow  solution  of  1,  which  means  that  the  wake  forces 
the  body  into  lock-in  over  a  wide  range  of  frequencies.  The  elliptical  cylinder  has  a 
potential  flow  added  mass  value  of  4,  but  sees  values  ranging  from  0  to  9.  It  reaches 
an  even  larger  range  (0  to  9)  because  of  the  doubly  sized  dimension  that  moves  in  the 
cross  flow  direction.  So  in  this  case,  the  wake  has  a  particularly  strong  effect  on  the 
VIV  response.  The  shed  vortices  are  coherent  and  therefore  force  lock-in  over  a  range 
of  frequencies,  despite  the  reduced  amplitude  of  response.  The  whisker,  however,  sees 
a  small  Cm  range:  3.3  to  3.6,  indicating  that  the  wake  plays  a  much  smaller  role  in 
causing  lock-in.  The  vortices  are  less  coherent,  so  the  wake  can  not  easily  force  the 
model  to  lock-in.  This  indicates  that  the  whisker  is  able  to  reduce  the  bandwidth  of 
response.  Also,  this  shows  that  the  added  mass  is  about  15%  larger  for  the  whisker 
than  the  cylinder,  giving  a  partial  explanation  for  why  the  whisker  geometry  is  so 
effective  at  reducing  vibrations. 
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Figure  2-11:  Vorticity  contours  of  a  2-D  slice  of  the  wake  behind  the  elliptical  cylinder 
(left)  and  whisker  (right)  models  undergoing  steady  tow.  A  standard  von  Karman 
vortex  street  appears  behind  the  ellipse,  but  the  whisker  decorrelates  the  vortex 
shedding,  causing  the  vortex  street  to  appear  much  farther  downstream,  where  the 
forcing  on  the  body  is  reduced. 

Damping  (£) 

The  circular  cylinder  sees  values  of  0  to  0.4  in  the  excitation  region,  the  elliptical 
cylinder  reaches  1.6,  and  the  whisker  goes  to  0.1.  These  values  also  illustrate  the 
difference  between  the  elliptical  cylinder  and  the  whisker.  Although  the  elliptical 
cylinder  reduces  the  amplitude  of  vibration,  lock-in  is  still  guaranteed  to  happen 
even  if  the  damping  reaches  the  high  value  of  1.6  at  low  amplitudes.  The  whisker  on 
the  other  hand,  only  guarantees  lock-in  if  the  system  damping  is  less  than  0.1  at  low 
amplitudes. 

2.3.2  PIV 

Steady  Flow 

First,  the  wakes  behind  the  whisker  and  elliptical  cylinder  undergoing  towing  motion 
with  no  vibrations  are  observed  using  PIV.  These  wakes  are  compared  at  a  single 
point  in  time  in  Figure  2-11.  The  whisker  wake  presents  a  von  Karman  street  farther 
downstream  than  the  ellipse.  Also,  the  vortices  in  the  whisker  wake  are  less  coherent. 
These  indicate  reduced  forcing  on  the  whisker  and  agree  with  the  conclusion  presented 
in  [58].  The  flow  patterns  match  well  with  those  shown  behind  a  bumpy  cylinder  by 


[99], 
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Figure  2-12:  Vorticity  contours  of  a  2-D  slice  of  the  wake  behind  the  elliptical  cylinder 
(left)  and  whisker  (right)  models  undergoing  the  forced  vibration  run  corresponding 
to  their  maximum  VIV  amplitude.  A  standard  von  Karman  vortex  street  appears 
behind  both  models. 

Forced  Vibrations 

The  wakes  behind  both  models  are  shown  in  Figure  2-12  under  conditions  corre¬ 
sponding  to  each  model’s  respective  maximum  VIV  amplitude  (elliptical  cylinder: 
A/d*  =  0.175,  fd*/U  =  0.27,  whisker:  A/d*  =  0.075,  fd*/U  =  0.33).  In  these  cases, 
both  models  form  a  von  Karman  vortex  street.  As  in  the  steady  case,  the  whisker’s 
wake  is  less  coherent  than  that  of  the  elliptical  cylinder.  The  difference,  however, 
in  coherence  and  formation  length  is  less  dramatic,  highlighting  the  importance  of 
including  motions. 

2.4  Results  (ct  >  0°) 

The  same  method  of  forced  vibrations  was  conducted  on  the  whisker  and  ellipti¬ 
cal  cylinder  models,  while  facing  the  flow  at  an  angle  ( a  >  0°).  Experiments  were 
conducted  in  increments  of  15°  up  through  90°.  Due  to  size  limitations  in  the  exper¬ 
imental  setup,  smaller  models  were  used  for  the  tests  at  a  >  45°. 

2.4.1  Force  Measurements 

Drag  ( Cd ) 

The  drag  on  the  whisker  at  different  angles  of  attack  is  shown  in  Figure  2-13.  During 
a  zero-amplitude  vibration  case,  the  drag  coefficient  is  approximately  1.2,  1.5,  and  2 
for  the  15,  30,  and  45°  cases,  respectively.  The  60,  75,  and  90°  cases  are  about  2.5. 
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The  elliptical  cylinder  responds  similarly  for  the  high  angle  cases  (45°  and  up). 


Both  models  respond  similarly  to  each  other,  which  is  at  values  higher  than  that 
of  the  cylinder.  The  benefit  of  using  the  whisker  geometry  to  reduce  drag  is  therefore 
lost  when  the  whisker  does  not  face  the  flow  in  the  streamlined  position. 


Energy  Transfer  ( C^v ) 

CijV  measurements  shown  in  Figure  2-14  reveal  that  vibrations  in  both  the  whisker 
and  elliptical  cylinder  increase  with  angle  of  attack.  The  whisker’s  amplitude  at 
a  =  0°  is  about  half  of  the  elliptical  cylinder’s  amplitude.  Up  through  a  =  30°, 
the  two  models  continue  to  respond  with  low  vibration  amplitude  ( A/d *  ~  0.2  or 
less).  About  45°,  the  amplitude  ratio  jumps  up  and  continues  to  increase  through 
a  =  90°,  where  both  models  respond  with  large  amplitude  (A/d*  ~  1.5).  This  value 
is  about  twice  as  high  as  that  experienced  by  a  circular  cylinder  (from  Figure  2-8: 
A/d*  ~  0.7).  The  whisker’s  undulations  are  too  small  relative  to  the  bluff  dimension 
to  aid  desynchronization  of  vortex  shedding  at  high  angles  of  attack.  Figure  2-15 
displays  this  trend  in  15°  increments. 

Use  of  the  whisker  geometry  as  a  vibration  suppression  device  in  scenarios  where 
flow  approaches  the  body  from  a  range  of  directions  is  therefore  unsuitable.  However, 
the  wide  spectrum  of  vibrational  response  (factor  of  ~  20)  may  be  advantageous  for 
the  seal  and  can  be  exploited  in  an  engineering  sensor. 


Effective  Added  Mass  ( Cm ) 

The  effective  added  mass  coefficient  Cm  of  the  whisker  and  ellipse  at  different  angles 
of  attack  is  shown  in  Figure  2-16.  The  results  indicate  that  the  models  experience 
VIV  over  a  wide  range  of  response.  At  90°,  the  whisker  and  elliptical  cylinder  have 
Cm  values  of  ~  0  to  — 2.  This  relatively  large  deviation  from  the  potential  flow  value 
of  1  shows  that  the  wake  adjusts  to  encourage  VIV. 
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Figure  2-13:  Drag  coefficient  during  forced  vibrations  for  (left  column)  ellipse  and 
(right  column)  whisker  at  (a,b)  15;  (c,d)  30;  (e,f)  45;  (g,h)  60;  (i,j)  75;  (k,l)  90°.  Con¬ 
tour  levels  for  15  and  30  °  are  in  increments  of  0.1.  Otherwise  they  are  in  increments 
of  0.5. 
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Figure  2-14:  C^v  contour  for  (left  column)  ellipse  and  (right  column)  whisker  at  (a,b) 
15;  (c,d)  30;  (e,f)  45;  (g,h)  60;  (i,j)  75;  (k,l)  90°.  Contour  levels  are  in  increments  of 

1. 
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Figure  2-15:  The  peak  amplitude  of  vibration  for  the  elliptical  cylinder  and  whisker 
increases  as  the  angle  of  attack  increases.  Re  =  5,300  —  7,030  from  a  =  0°  —  45°  and 
Re  =  2,530  -  3,210  from  a  =  60°  -  90°. 

2.4.2  PIV 

Steady  Flow 

Flow  visualizations  of  the  whisker  and  ellipse  being  towed  rigidly  down  a  tank  at 
various  angles  is  shown  in  Figure  2-17.  The  whisker  is  seen  to  desynchronize  the 
wake  more  than  the  the  ellipse,  especially  at  low  angles  of  attack  (0,15,30°).  The 
vortex  structures  right  next  to  the  body  is  elongated,  which  reduces  the  oscillatory 
effect  of  the  normal  alternating  vortex  shedding  pattern.  At  higher  angles  of  attack 
(45,60,90°),  the  wake  pattern  looks  similar  between  the  two  models,  however  the  wake 
width  is  narrower  for  the  whisker,  especially  in  the  90°  case. 


1-D  Free  Vibrations  -  Dye  Visualizations 

1-D  free  vibration  experiments  with  the  30x  scale  whisker  model  were  conducted  by 
my  colleagues  Yahya  Modarres-Sadeghi  and  Banafsheh  Seyed-Aghazadeh  at  UMass- 
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Figure  2-16:  Effective  added  mass  coefficient  during  forced  vibrations  for  (left  column) 
ellipse  and  (right  column)  whisker  at  (a,b)  15;  (c,d)  30;  (e,f)  45;  (g,h)  60;  (i,j)  75; 
(k,l)  90°.  Contour  levels  are  in  increments  of  0.5. 
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Figure  2-17:  PIV  of  the  wake  behind  an  elliptical  cylinder  (left)  and  whisker  (right) 
being  rigidly  towed  at  different  angles  of  attack.  Re  =  5,300. 
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Figure  2-18:  The  wake  shed  behind  the  freely  vibrating  whisker  for  the  (left)  45  and 
(right)  90°  cases.  (Experiments  conducted  by  B.  Seyed-Aghazadeh  and  Y.  Modarres- 
Sadeghi) 

Amherst  through  use  of  a  linear  air  bearing  stage.  Snapshots  from  dye  visualizations 
of  the  wakes  are  shown  in  Figure  2-18  for  the  whisker  at  a  =  45°  and  90°.  Similar 
to  the  steady  PIV  results  shown  in  Figure  2-17,  the  wake  under  these  high  angle 
conditions  is  characterized  by  coherent  vortices  and  a  large  wake  width. 


2.5  Discussion 

2.5.1  Whisker  in  Frontal  Flows 

The  harbor  seal  whisker  geometry  has  been  dynamically  tested,  confirming  the  state¬ 
ment  originally  put  forward  by  [58]  that  the  geometry  passively  reduces  VIV.  Exper¬ 
iments  of  1-D  forced  vibrations  are  used  to  find  the  zero  contour  of  CjA} .  The  peak 
value  of  this  contour  indicates  a  factor  of  4  and  10  lower  amplitudes  of  vibration  for 
an  elliptical  cylinder  and  whisker,  respectively,  in  comparison  with  a  baseline  circular 
cylinder  at  the  same  Reynolds  number.  There  is  some  variation  in  the  whisker’s  Clv 
value  over  time.  This  variation  is  a  result  of  the  whisker  requiring  a  longer  time  to 
reach  lock-in  than  the  cylinder  or  elliptical  cylinder. 

The  effective  added  mass  coefficient  and  damping  coefficient  indicate  that  the 
whisker  geometry  reduces  the  lock-in  region  compared  to  the  elliptical  cylinder.  For 
both  of  these  geometries,  the  still-water  effective  added  mass  in  the  y-direction  is  large 
because  of  vortex  shedding  along  the  body,  leading  to  much  smaller  forcing  terms  in 


phase  with  the  acceleration  of  the  body  and  therefore  smaller  variation  in  the  effective 
added  mass.  Both  shapes  exhibit  large  damping  and  small  amplitudes  of  motion  when 
excited  in  the  y-direction  because  of  their  elongated  body  shapes,  however  the  whisker 
has  a  much  smaller  bandwidth  for  excitation.  Flow  visualizations  show  that  lock-in 
is  reduced  as  a  result  of  the  desynchronized  wake  and  increased  formation  length. 


2.5.2  Whisker  at  Angles  of  Attack 

VIV  reduction  does  not  apply  when  flow  approaches  the  whisker  from  larger  angles 
of  attack.  Beyond  a  =  45°,  the  whisker  responds  similarly  to  a  smooth  elliptical 
cylinder,  which  is  with  an  amplitude  higher  than  that  of  a  circular  cylinder.  At 
a  =  90°,  the  whisker  vibrates  1.5  times  its  projected  diameter.  Flow  visualizations 
show  that  the  undulatory  features  that  desynchronize  vortex  shedding  in  the  whisker’s 
streamlined  configuration  are  too  small  relative  to  the  bluff  dimension  to  enable  the 
same  effect.  The  wake  looks  similar  to  that  of  the  elliptical  cylinder  at  the  a  >  45° 
cases,  which  has  strong  coherent  vortices. 

The  high  VIV  response  as  the  body  becomes  more  bluff  with  respect  to  the  flow 
removes  the  possibility  of  using  this  shape  as  a  vibration  suppression  device  in  flows 
with  varying  angle.  Designs  of  whisker-inspired  geometries  that  are  omnidirectional 
were  made  in  an  attempt  to  create  a  shape  that  experiences  low  VIV  regardless  of  the 
angle.  See  Appendix  B  for  details  of  this  study.  These  models  were,  however,  unsuc¬ 
cessful  at  re-creating  the  extreme  reduction  capabilities  of  the  whisker,  highlighting 
the  importance  of  the  elliptical  cross  section  in  reducing  VIV. 

Although  the  whisker  geometry  is  not  suitable  for  use  as  a  VIV  reduction  in  cases 
where  flow  will  approach  the  object  with  an  a  >  45°,  it  holds  great  promise  for  use  as 
a  flow  sensor.  The  VIV  reduction  in  frontal  flows  allows  easier  distinction  of  flows  that 
approach  with  a  side  component.  Both  seals  and  engineered  devices  could  use  this  for 
heightened  sensing  of  the  direction  of  oncoming  flow.  This  evidence  has  encouraged 
sensor  development,  which  is  described  in  the  first  half  of  Chapter  3. 
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2.5.3  Note  on  the  geometry 

This  data  is  intended  to  describe  the  VIV  properties  of  a  “whisker-inspired”  geometry. 
The  geometry  used  here  can  likely  be  further  optimized  for  VIV  reduction.  In  fact, 
it  was  noticed  that  the  morphology  of  the  harbor  seal  whisker  (in  both  real  whisker 
specimens  obtained  from  the  New  England  Aquarium  as  well  as  the  specimen  reported 
in  [45])  is  different  than  the  geometry  reported  in  [58],  which  was  replicated  in  this 
chapter.  The  undulation  wavelength  appears  to  be  nearly  double  the  wavelength 
used  here.  In  Chapter  4  both  the  current  geometry  and  one  that  captures  the  longer 
wavelength  will  be  tested. 
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Chapter  3 


Flow  Sensor  Design 


How  can  these  unique  vibration  properties  be  realized  in  a 

physical  prototype? 

3.1  Introduction 

It  is  desired  to  develop  a  sensor  that  mimics  the  capability  of  harbor  seals  to  de¬ 
tect  vortex  wakes.  Since  the  mechanism  by  which  the  seals  do  this  is  unknown,  we 
approach  this  goal  in  steps.  First,  a  device  that  conducts  basic  steady  flow  measure¬ 
ments  is  created.  It  is  designed  to  measure  flow  velocity  and  exploit  the  discovered 
seal  whisker  vibration  properties  (Chapter  2)  for  heightened  flow  direction  sensing. 
Second,  the  design  is  refined  and  adjusted  for  use  as  a  platform  to  systematically 
explore  wake-detection  capabilities.  That  detailed  study  is  then  described  in  Chapter 
4. 

The  current  work  utilizes  the  discovered  mechanism  of  frontal-flow  vibration  re¬ 
duction  to  achieve  clear  drag-based  readings  for  correlation  to  oncoming  flow  velocity. 
Towing  tank  tests  are  used  to  measure  the  sensor’s  range  of  vibrational  response  over 
various  flow  directions.  Field  experiments  test  the  sensor’s  usability  at  sea.  Finally, 
a  refined  device  is  put  forth  as  a  platform  to  conduct  tests  on  wake  detection. 
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3.1.1  Previous  Work 


Flow  Velocity  Sensors 

The  first  step  seeks  to  create  a  device  that  measures  flow  velocity.  A  variety  of 
such  sensors  already  exist.  These  span  the  use  of  acoustic-,  magnetic-,  mechanical-, 
and  pressure-based  sensing  techniques.  AUVs  generally  utilize  Doppler  velocity  logs 
(DVLs),  which  measure  velocity  based  on  travel  time  of  acoustic  pings.  These  active 
devices  demand  a  significant  portion  of  the  vehicle’s  battery  power.  Experiments  with 
real  seals  suggest  that  they  may  use  their  vibrissae  to  detect  and  maintain  their  speed 
[111].  A  whiskerlike  sensor,  which  is  mechanical-based  as  in  the  real  seals,  could  serve 
as  a  low-power,  low-cost  alternative  to  the  standard  devices  in  certain  scenarios. 

Previous  whiskerlike  or  hairlike  flow  sensors  have  been  developed  [6,  32,  52,  127]. 
The  “whiskers”  used  have  generally  been  limited  to  cylindrical  rods.  The  drag-based 
response  of  these  rods  can  indeed  be  used  to  determine  flow  velocity  and  direction. 

Vortex  Wake  Detection  Sensors 

An  important  limitation  to  the  existing  whiskerlike  and  hairlike  flow  sensors  is  that 
their  uniform  circular  cross  section  precludes  them  from  detecting  unsteady  flow  fea¬ 
tures.  Since  such  bluff  bodies  experience  significant  vortex-induced  vibrations  above 
Re  ~  50,  the  baseline  signal  they  emit  will  include  strong  oscillations,  which  can 
reach  amplitudes  of  about  one  cylinder  diameter  [49].  Any  additional  response  upon 
encountering  features  such  as  vortex  wakes  may  then  be  indistinguishable  from  the 
baseline  vibrations. 

Existing  flow  sensors  could  be  modified  to  directly  target  measurement  of  vortex 
wakes  by  actively  emitting  signals,  such  as  acoustic-based  ones.  Passive  measurement 
options  do  exist  through  instruments  such  as  hot-wire  anemometers  and  pressure 
probes.  Such  devices  have  often  been  used  in  laboratory  settings  to  measure  flow 
velocity  and  can  be  used  to  detect  vortices  [36,  113,  126]. 

The  remaining  barrier  is  in  conducting  passive  detection  while  moving  forward  at 
significant  speed.  Devices  reaching  into  the  flow  while  mounted  externally  to  an  AUV 
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would  need  to  function  despite  the  presence  of  steady  oncoming  flow,  which  can  lead 
to  VI V  on  the  sensor.  This  work  therefore  seeks  to  develop  an  instrument  for  passive 
vortex  wake  detection  via  high  wake- “signal” -to-baseline- “noise”  response. 

3.2  Flow  Sensor  -  Version  1 

The  design,  calibration,  and  testing  of  the  Version  1  device  was  first  presented  in  [10] 
and  is  described  again  here. 

3.2.1  Design 

This  device  was  designed  to  operate  with  a  similar  sensing  mechanism  to  that  used  by 
the  real  seals.  No  nerves  are  present  inside  the  length  of  the  real  whiskers.  Instead,  the 
sensing  takes  place  first  in  the  seals’  cheek  as  the  whiskers  bend  in  the  flow.  Similarly, 
the  sensor  here  was  designed  to  measure  deflections  at  the  base  of  the  whisker. 

The  complete  sensor  design  is  shown  in  Figure  3-1.  The  whisker  section  is  the 
same  model  used  in  Chapter  2  for  the  larger  angle  of  attack  cases.  This  part  was 
fabricated  via  stereolithography  at  a  scale  of  thirty  times  the  real  whisker  diameter 
(whose  dimensions  are  taken  from  [58]),  giving  it  an  average  crossflow  diameter  (d)  of 
1.59cm,  when  facing  the  flow  in  the  streamlined  direction.  The  length  of  the  whisker 
section  was  chosen  as  27.5 cm,  which  encompasses  five  full  wavelengths  of  undulations. 
The  whisker  protrudes  from  a  housing,  which  contains  Flexpoint  Bend  Sensors  that 
change  resistance  based  on  the  amount  they  are  deflected.  The  underlying  Bend 
Sensor  configuration  is  based  on  previous  work  done  by  a  colleague,  whose  design 
details  are  outlined  in  [130]. 

A  closer  view  of  the  Bend  Sensors  is  given  in  Figure  3-2a.  The  four  are  rigidly 
mounted  at  one  end  in  a  holder  and  allowed  to  bend  at  the  other  end  as  a  square  plate 
pushes  on  them.  That  plate  is  attached  to  the  end  of  a  rod  that  is  connected  to  the 
whisker.  The  starting  position  of  the  plate  is  such  that  it  is  just  touching  each  of  the 
Bend  Sensors,  allowing  any  small  deflection  to  produce  a  change  in  voltage  output. 
The  sequence  is:  Deflection  of  the  whisker  pushes  on  the  plate,  which  pushes  on  the 
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Figure  3-1:  The  whisker  (yellow)  mounted  to  a  pressure  housing,  which  contains  four 
Bend  Sensors  (dark  orange). 


Figure  3-2:  (a)  A  square  plate  (in  black)  at  the  end  of  a  rod  connected  to  the  whisker 
pushes  on  the  four  Bend  Sensors  as  the  whisker  deflects,  (b)  A  silicone  diaphragm 
(in  orange)  between  the  whisker  and  a  retaining  ring  allows  flexible  motion. 


Bend  Sensors,  which  change  resistance  and  thereby  provide  deflection  information. 
This  section  describes  how  the  magnitude  of  the  output  can  then  be  translated  to 
direction  and  velocity  of  the  oncoming  flow. 

An  interesting  design  parameter  was  encountered  in  requiring  the  whisker  to  flex¬ 
ibly  move  in  response  to  the  oncoming  flow  but  simultaneously  maintain  a  watertight 
seal  with  the  housing.  The  chosen  solution  was  to  use  a  flexible  base  with  the  dual 
function  of  allowing  movement  of  the  whisker  and  serving  as  a  waterproof  barrier  to 
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the  pressure  housing.  A  silicone  diaphragm  (Shore  50A  hardness)  was  sandwiched  be¬ 
tween  the  end  of  the  whisker  and  a  plastic  piece  (Figure  3-2b)  via  bolts.  The  whisker 
is  thus  firmly  attached  to  the  rest  of  the  system,  but  the  diaphragm  allows  motion. 
The  lack  of  holes  from  outside  to  inside  the  housing  allows  for  a  watertight  seal. 

This  flexible  membrane  will  deform  under  water  pressure  when  the  sensor  is  sub¬ 
merged.  Such  a  shape  change  would  cause  additional  changes  in  readings  from  the 
Bend  Sensors.  As  such,  it  was  chosen  to  fill  the  inside  of  the  housing  with  oil,  provid¬ 
ing  a  back  pressure  to  keep  the  diaphragm  in  a  stationary  position  until  the  whisker 
deflects. 

Electronics 

A  separate  housing  is  used  to  store  the  electronics  since  the  main  housing  is  filled 
with  oil.  A  standard  pressure  housing  was  built  to  hold  an  amplification  circuit,  a 
data  logger,  and  a  power  supply.  Wires  from  the  Bend  Sensors  go  from  the  main 
housing,  through  a  Subconn  connector  and  cable  to  the  electronics  housing,  to  the 
amplification  circuit,  and  then  into  the  data  logger  (Figure  3-3).  The  DATAQ  DI-710 
data  logger  records  data  onto  an  SD  Card  at  a  sampling  frequency  of  1440 Hz.  A 
battery  pack  is  used  to  power  the  data  logger. 

A  circuit  that  provides  filtering  and  amplification  to  the  Bend  Sensors  was  designed 
and  assembled.  The  diagram  is  outlined  in  Figure  3-4.  LM324  op-amps  are  used  in  a 
slightly  expanded  version  of  a  voltage  divider.  Resistor  values  were  chosen  to  closely 
match  the  base  resistance  value  of  the  corresponding  Bend  Sensor.  Small  capacitors 
are  included  in  a  RC  low-pass  filter  configuration. 


Figure  3-3:  The  electronics  that  are  stored  inside  a  separate  housing:  circuit  and  data 
logger.  The  battery  pack  is  not  pictured. 
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Figure  3-4:  Circuit  diagram.  The  Bend  Sensor  signal  lines  are  filtered  and  amplified. 

3.2.2  Calibration 

The  voltage  output  from  the  Bend  Sensors  was  related  to  a  physical  deflection  length. 
This  was  done  by  manually  deflecting  the  whisker  in  equal  increments  and  matching 
the  output  voltage  to  the  corresponding  tip  deflection  distance  determined  from  the 
adjacent  ruler.  Figure  3-5  shows  the  setup.  The  sensor  was  then  rotated  90°  and 
the  procedure  was  repeated  three  more  times  to  collect  data  for  all  four  of  the  Bend 
Sensors. 


Figure  3-5:  Setup  for  bench-top  calibration.  Deflection  of  the  whisker  was  conducted 
in  increments,  using  the  ruler  on  the  left  as  a  reference. 

Figure  3-6  shows  the  voltage  output  for  each  of  the  four  Bend  Sensors  over  the 
incremental  deflections.  Trendlines  were  best  fit  to  these  curves  and  the  corresponding 
equations  are  shown  in  Equation  3.1.  V  is  voltage  (Volts)  and  dt  is  tip  displacement 
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Figure  3-6:  Voltage  output  versus  whisker  tip  displacement  for  each  of  the  four  Bend 
Sensors. 


(mm).  Subscript  numbers  indicate  which  Bend  Sensor  the  variable  corresponds  to. 

Vi  =  0.00084  +  0.0015dti  +  1.4972 
V2  =  0.0034 dt2  +  1.2651 

(3.1) 

V3  =  0.00134  +  0.0008dt3  +  1.3683 
V4  =  0.0065 dt4  +  1.3292 

The  coefficient  of  determination  ( R 2)  for  each  of  these  equations  is  0.997,  0.993, 
0.999,  and  0.989  in  sequential  order.  Sensors  2  and  4  can  be  characterized  by  a  linear 
relationship  between  voltage  and  displacement,  but  Sensors  1  and  3  are  best  described 
with  a  quadratic  one.  This  variation  from  Bend  Sensor  to  Bend  Sensor  emphasizes 
the  need  for  quantifying  the  individual  responses. 
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3.2.3  Towing  Tests 


The  whisker  response  (vibration  and  deflection)  to  various  flow  conditions  was  quan¬ 
tified  and  collected.  This  was  done  by  conducting  towing  tests  in  the  large  tank  of 
the  MIT  Towing  Tank  Lab,  where  the  sensor  was  mounted  to  the  carriage  with  the 
whisker  model  vertically  submerged  (Figure  3-7).  It  was  towed  at  different  velocities 
(U  =  0.2,  0.4,  0.6,  0.8 m/s)  and  different  angles  of  attack  (a  =  0, 30,  45,  60, 90°).  As  de¬ 
scribed  in  the  previous  chapters,  zero  angle  of  attack  refers  to  the  orientation  where 
the  longer  diameter  is  exactly  in-line  with  the  oncoming  flow.  The  corresponding 
Reynolds  number  is  Re  =  3, 180  to  12,  720  based  on  the  average  crossflow  diameter 
(d)  at  a  =  0°.  Figure  3-7  also  displays  the  orientation  of  the  Bend  Sensors.  Note 
that  as  the  angle  of  attack  changes,  the  sensors  which  are  in  the  crossflow  and  inline 
direction  also  change. 

Plots  of  voltage  versus  time  were  generated  for  each  run  and  compiled  for  later 
use  as  a  lookup  table  to  invert  output  voltage  to  the  corresponding  flow  information. 
A  selection  of  these  runs  are  shown  in  Figure  3-8.  Each  towing  run  shows  a  change  in 
voltage  once  the  carriage  begins  moving  forward,  indicating  a  Bend  Sensor  deflection 
in  that  corresponding  direction.  Then  oscillations  develop  in  the  crossflow  Bend 
Sensors,  and  finally  the  voltage  ramps  back  to  the  initial  value  as  the  carriage  stops. 

As  discovered  in  Chapter  2,  the  whisker  geometry  vibrates  preferentially  with 
direction.  As  the  angle  of  attack  increases,  the  crossflow  vibrations  increase.  At 
a  =  0°,  the  voltage  of  the  crossflow  sensors  (1  and  2,  for  this  angle)  oscillated  with 
extremely  small  amplitude  (Figure  3-8a).  The  crossflow  sensors  in  the  90°  case  (3 
and  4),  however,  show  significant  oscillations  (Figure  3-8f).  The  inline  sensors  for  0° 
and  90°  experience  a  mean  deflection  but  very  little  oscillations  (Figures  3-8b  and 
3-8e).  For  angles  between  those  two  extremes,  such  as  45°,  the  crossflow  oscillations 
are  captured  in  a  combination  of  all  four  Bend  Sensors  (Figures  3-8c  and  3-8d). 
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Figure  3-7:  Towing  test  setup.  (Left)  The  sensor  is  mounted  to  the  carriage,  and  a 
data  logger  records  voltage  output.  Note:  The  wake  generator  (circular  cylinder)  pic¬ 
tured  was  not  used  for  the  experiments  shown  in  this  Chapter.  (Right)  Schematic  of 
the  inside  of  the  sensor  housing  (top  view) .  The  Bend  Sensors  (purple)  are  numbered 
with  this  convention  throughout.  Note:  For  a  =  0°,  Bend  Sensors  1  and  2  are  in  the 
crossflow  direction,  but  as  a  is  increased  to  90°,  Bend  Sensors  3  and  4  become  fully 
in  the  crossflow  direction. 


Sensing  Procedures 

The  results  from  the  calibration  runs  are  then  organized  so  that  a  given  sensor  reading 
can  be  converted  to  the  corresponding  flow  velocity  and  direction.  This  is  done  by 
first  extracting  three  main  features  from  the  calibration  data:  deflection  in  the  drag 
direction  and  vibration  amplitude  in  both  axes.  These  three  variables  are  converted 
from  voltage  to  a  physical  length  of  motion  via  Equation  3.1. 

The  oncoming  flow  angle  can  be  determined  by  comparing  vibration  magnitudes 
in  the  two  axes.  The  Bend  Sensor  that  lies  in  the  crossflow  direction  will  vibrate  with 
higher  amplitude  than  the  adjacent  Bend  Sensor.  Figure  3-9a  depicts  this  trend  of 
diverging  values  between  the  inline  and  crossflow  motions  as  the  angle  of  attack  in- 
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Figure  3-8:  Time  traces  of  the  Bend  Sensor  outputs.  Black  vertical  lines  indicate  start 
and  end  of  the  carriage  motion.  (Top  row)  a  =  0°:  minimal  change  occurs.  (Middle 
row)  a  =  45°:  ramps  and  oscillations  in  voltage  are  seen.  (Bottom  row):  a  =  90°: 
larger  oscillations  occur  in  the  crossflow  (Sensors  3,4)  and  larger  ramps  occur  inline 

(Sensors  1,2).  Sensors  1  and  4  are  shifted  by  a  constant  value  to  consolidate  plots. 
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creases.  Note  that  the  values  presented  here  are  not  universal  values  for  the  whisker’s 
response.  The  exact  values  will  change  with  the  mass-damping  (m*Q  value  of  the 
setup.  A  more  quantitatively  accurate  comparison  is  found  through  the  forced  vibra¬ 
tion  data  shown  in  Chapter  2.  The  plots  here  are  intended  merely  to  confirm  the 
trends. 

A  comparison  of  the  two  amplitudes  therefore  provides  the  first  piece  of  informa¬ 
tion:  angle  of  oncoming  flow.  However,  it  is  important  to  point  out  that  the  vibration 
amplitudes  vary  as  the  flow  velocity  changes.  The  divergence  between  in-line  and 
crossflow  amplitudes  remains  across  all  velocities  tested  here. 

The  amplitude  values  shown  here  were  determined  from  the  corresponding  Bend 
Sensors  for  that  particular  angle  of  attack  (i.e.,  Sensor  2  for  crossflow  vibration  of 
a  =  0°  and  Sensor  4  for  crossflow  vibration  of  a  =  90°).  For  calibration  runs  con¬ 
ducted  at  a  =  30°, 45°  and  60°,  no  Bend  Sensor  was  aligned  exactly  with  either 
crossflow  or  in-line  flow  directions.  Measurements  in  the  in-line  direction  were  there¬ 
fore  approximated  by  A(SensorX )  =  ThTg^Tl)  and  measurements  in  the  crossflow 
direction  were  approximated  by  A(SensorX)  =  •  A(SensorX)  is  the  am¬ 

plitude  of  motion  that  would  be  experienced  by  a  Bend  Sensor  appropriately  aligned 
with  the  corresponding  angle  of  attack  (aj.  This  geometric  relation  converts  the  mo¬ 
tion  (A(Sensorl)  or  A(SensorA))  measured  by  Sensors  1  or  4,  whose  locations  are 
shown  in  Figure  3-7,  to  an  estimated  value. 

Next,  the  sensor  can  be  used  to  determine  the  flow  velocity  (77).  As  shown  in 
Figure  3-9b,  the  inline  deflection  magnitude  (found  from  the  aft  Bend  Sensor  in  each 
case)  increases  with  velocity.  It  is  also  dependent  on  angle  of  attack,  because  as  a 
increases,  the  area  facing  the  flow  becomes  less  streamlined.  As  such,  higher-angle 
cases  display  a  higher  in-line  deflection.  So  once  the  angle  of  attack  is  determined, 
the  output  voltage  in  the  in-line  direction  can  be  directly  converted  to  flow  velocity. 

Accuracy 

Based  on  the  bench-top  calibration,  the  sensor  has  a  resolution  of  approximately 
1  mm  deflection  at  the  whisker  base.  Correspondingly,  it  can  safely  distinguish  angle 
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Figure  3-9:  (a)  Nondimensional  amplitude  of  inline  and  crossflow  oscillations  for  the 
whisker  at  various  angles  of  attack  to  the  oncoming  flow.  Each  amplitude  (A)  is 
divided  by  the  average  cross-stream  diameter  at  a  =  0°  :  d  =  1.59 cm.  (b)  Deflection 
in  the  drag  (inline)  direction  vs.  velocity  (Sensor  3  and  Sensor  1  for  a  =  0°  and  90°, 
respectively) .  The  next  design  iteration  will  include  more  sensitive  measuring  devices 
so  that  the  low-speed,  low-angle  of  attack  results  can  be  better  distinguished  above 
the  inherent  system  noise. 


of  attack  with  30°  accuracy.  Note  that  the  results  presented  here  are  not  intended  to 
depict  a  thorough  study  of  the  VIV  response  of  the  whisker  geometry.  That  study  was 
shown  in  Chapter  2,  in  which  measurement  devices  with  higher  resolution  are  used 
and  where  boundary  conditions  and  mass-damping  are  carefully  considered.  Rather, 
this  work  presents  a  device  that  captures  the  width  of  the  vibration  response  range 
and  then  utilizes  it  for  heightened  sensing. 


3.2.4  Research  Cruise  Deployment 

To  determine  the  sensor’s  performance  and  practicality  as  an  ocean-ready  instrument, 
it  was  tested  at  sea.  The  sensor  was  deployed  from  the  R/V  Wecoma  on  September 
24,  2011  during  the  UNOLS  Chief  Scientist  Training  Cruise  near  Astoria  Canyon, 
off  the  coast  of  Oregon  (46. 22° A)  125.22 °W).  The  sensor  was  mounted  to  the  frame 
holding  the  CTD  rosette,  as  shown  in  Figure  3-10.  Both  pressure  housings  were 
strapped  to  the  outer  frame  and  the  whisker  segment  protruded  into  the  flow,  largely 
unobstructed.  The  instrument  was  lowered  to  a  depth  of  20m  (at  1405  GMT),  held 
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in  place  for  10  minutes,  and  then  brought  back  to  the  surface  (at  1415  GMT). 


Figure  3-10:  Field  test  setup.  The  sensor  is  mounted  to  the  CTD  rosette  in  a  manner 
that  allows  the  whisker  to  encounter  unobstructed  flow.  The  view  is  from  the  port 
side  of  the  ship.  (Photo  courtesy  of  Paul  Suprenand) 

The  sensor  successfully  collected  data  and  maintained  the  watertight  seal.  The 
velocity  measurements  from  the  ship- mounted  ADCP  (RDI  300 kHz  Ocean  Mariner, 
sampling  frequency  =  1/120 Hz)  are  taken  as  the  results  to  be  compared  against. 
That  data,  shown  in  Figure  3-1  la,  was  processed  through  UHDAS  CODAS.  Data 
from  two  of  the  Bend  Sensors  during  the  10-minute  run  was  sent  through  a  low-pass 
Butterworth  filter  and  are  shown  in  Figure  3-llb.  Sensor  3  shows  oscillations  and 
various  ramps  in  deflection.  Sensor  2,  shows  deflection,  but  oscillations  are  minimal 
and  are  mostly  a  result  of  noise.  This  suggests  that  Sensor  2  is  in  the  mean  flow’s  inline 
direction  and  Sensor  3  in  the  crossflow  direction.  Sensor  3  was  roughly  aligned  with 
the  ship’s  bow,  and  the  ship’s  heading  during  the  time  of  deployment  ranged  from 
-120°  to  -140°.  The  ADCP  data:  (+)  u  and  (-)  v  values  indicates  a  southeastwardly 
flow.  This  is  also  what  is  the  whisker  sensor  data  indicates,  as  Sensor  2  should 
therefore  roughly  be  in-line  with  the  southeast  direction  (Figure  3-12). 

The  amplitude  of  the  in-line  deflection  can  be  used  to  determine  the  flow  velocity. 


115 


Figure  3-11:  (a)  Velocity  measurements  at  20m  depth  from  the  ADCP  and  (b)  voltage 
output  from  the  whisker- inspired  sensor. 


N 


Figure  3-12:  Orientation  of  the  ship  and  sensor.  The  ship  had  a  heading  of  -120°  to 
-140°  during  the  lOrnin  deployment.  The  arrow  indicates  the  general  flow  direction 
shown  by  the  ADCP.  The  numbers  indicate  the  orientation  of  the  corresponding  Bend 
Sensors. 

At  the  beginning  of  the  run,  the  value  is  around  A/d  =  0.05,  which  from  Figure  3-9b 
indicates  ~  0.2 m/s,  which  matches  the  same  rough  value  expected  from  Figure  3-lla. 
Some  drift  appears  towards  the  end  of  the  run,  however,  so  this  value  may  be  of  low 
accuracy. 

Discrepancies  between  the  two  sensors  can  be  attributed  to  a  few  factors.  The 
motion  of  the  CTD  rosette  frame,  which  is  hung  on  a  relatively  thin  cable,  does 
not  match  the  motion  of  the  ship  exactly.  Also,  the  flow  field  being  measured  was 
disturbed  to  some  degree  by  the  arrival  of  the  large  frame  directly  preceding  the 
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measurements.  Finally,  the  whisker  sensor  does  not  account  for  w  (vertical)  velocity. 
Future  tests  should  hold  the  sensor  in  position  longer  to  allow  effects  from  preceding 
motions  to  die  out,  and  they  should  include  a  gyro  and  accelerometers  on  the  sen¬ 
sor/frame  to  account  for  any  motion  of  the  sensor  itself.  Nonetheless,  these  results 
serve  as  preliminary  evidence  that  the  sensor  can  be  used  to  detect  orientation  and 
velocity  in  the  field. 
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3.3  Flow  Sensor:  Version  2 


3.3.1  Design 

A  design  iteration  on  the  first  prototype  was  implemented  in  order  to  1)  increase 
device  performance  and  2)  enable  systematic  studying  of  its  flow-induced  vibration 
characteristics.  The  Bend  Sensors  used  in  Version  1  were  found  to  plastically  deform 
after  prolonged  use.  Strain  gauges  were  chosen  as  a  superior  alternative  by  providing 
low  hysteresis  and  higher  resolution.  The  design  was  also  changed  to  allow  direct 
measurement  of  the  response  in  both  axes.  This  version  can  therefore  be  employed 
as  a  free  vibration  setup  with  measurements  in  both  degrees  of  freedom. 

The  new  design  is  shown  in  Figure  3-13.  It  includes  the  same  30x-scale  plastic 
whisker  model  (average  crossflow  diameter  d  =  1.59 cm,  length  l  =  27.5cm)  mounted 
in  a  cantilevered  configuration  as  before.  Now  the  model  is  connected  to  a  four-armed 
flexing  plate  with  strain  gauges  (Omega,  Model  KFH-03-120).  Similar  cruciform 
flexure-based  designs  have  been  used  in  other  types  of  sensors  [9,  28,  101]. 

As  the  whisker  and  flow  field  interact,  the  four  arms  bend  accordingly,  allowing 
free  vibration  in  both  in-line  (x)  and  crossflow  (y)  directions.  A  voltage  signal  is  read 
from  the  strain  gauges,  and  this  is  then  calibrated  to  the  deflection  at  the  model  tip. 
This  system  is  designed  for  detailed  towing  tank  tests,  rather  than  use  in  the  field. 
To  this  end,  it  is  not  designed  to  be  submerged.  The  strain  gauges  and  flexure  remain 
above  the  waterline  during  tests.  After  the  strain  gauges  are  applied  to  the  plate,  the 
entire  plate  is  painted  with  Rust-Oleum  to  mitigate  rusting. 

Figure  3-14(a)  depicts  the  flexure’s  response  to  deflection  of  the  whisker.  As  the 
rigid  model  moves  in  response  to  the  flow,  the  flexure  bends,  inducing  strain  on  its 
arms.  One  strain  gauge  is  placed  on  the  top  side  of  one  arm  in  each  axis.  These 
two  measurements  can  then  be  used  to  determine  how  much  the  whisker  is  moving 
in  each  respective  axis. 

A  top-down  view  of  the  flexure  device  is  provided  in  Figure  3-14(b).  A  four¬ 
armed  plate  of  1095  spring  steel  (thickness  =  0.38  mm)  was  fabricated  using  a  water 
jet.  Holes  in  the  middle  allow  connection  to  the  base  of  the  whisker  model.  The 
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Figure  3-13:  Isometric  view  of  the  wake  sensor,  which  consists  of  a  plastic  whisker 
model  mounted  on  a  four-armed  flexing  plate.  Strain  gauges  on  these  arms  are 
mounted  on  these  arms  for  measurement  of  the  whisker’s  deflections. 


end  of  each  arm  is  mounted  to  a  plastic  pin  and  bushing.  These  bushing  assemblies 
are  designed  to  reduce  coupling  between  the  x  and  y  axes.  As  the  plates  in  one 
axis  deflect,  the  bushings  on  the  off-axis  rotate,  preventing  the  off-axis  plates  from 
deflecting  unnecessarily. 


(a)  (b) 


Figure  3-14:  (a)  A  Solidworks  simulation  of  the  strain  induced  on  the  flexure  as  the 
whisker  deflects,  (b)  Close-up  of  the  flexure  device.  A  spring  steel  plate  made  using 
the  water  jet  is  attached  to  a  bushing  (clear)  and  pin  (gray)  at  the  end  of  each  arm. 
This  design  reduces  coupling  between  the  two  axes  by  allowing  rotation,  rather  than 
deflection  in  the  off-axis. 
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The  components  used  for  data-collection  are  shown  in  Figure  3-15.  The  strain 
gauges  are  arranged  in  a  Wheatstone  bridge  through  the  use  of  resistors  soldered  to 
an  external  board.  Their  output  is  acquired  with  a  PhidgetBridge  1046  Input  Board 
(sampling  frequency  =  60  Hz,  gain  =  16),  which  applies  basic  filtering  and  sends 
voltage  data  to  a  lab  desktop  via  USB. 


Figure  3-15:  The  data-collection  electronics.  The  strain  gauges  are  conFigured  in  a 
Wheatstone  bridge,  using  additional  resistors.  A  PhidgetBridge  is  used  to  filter  and 
transmit  the  voltages  to  a  computer  through  USB. 


3.3.2  Calibration 

The  flexure’s  output  voltage  was  calibrated  to  an  equivalent  whisker-tip  deflection. 
Bench  top  measurements  were  conducted  for  this  calibration.  The  sensor  was  held 
in  place  in  air  and  deflected  in  5 mm  increments.  The  output  voltage  (V)  at  each 
deflection  was  recorded  and  compiled.  Figure  3-16  shows  the  result  in  the  x  direction. 
The  calibration  curve  is  taken  to  be  the  fit  through  these  points.  Both  axes  were 
found  to  have  linear  calibration  curves,  with  the  coefficient  of  determination  (R2) 
value  being  0.98  and  0.97  in  x  and  y  respectively. 

The  natural  frequency  (/„)  of  the  whisker  sensor  was  found  by  conducting  pluck 
tests  and  calculating  the  zero-crossing  frequency  of  the  signals.  Figure  3-16  shows  the 
time  series  reading  of  the  pluck  test  in  air  for  the  y  direction.  In  air,  the  natural  fre- 
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Figure  3-16:  (Left)  A  linear  fit  was  placed  through  data  collected  for  different  tip 
deflections  of  the  whisker  in  x.  (Right)  A  pluck  test  in  air  was  used  to  find  the 
natural  frequency  of  the  whisker  in  y.  The  values  here  (Vyn,Vxn)  are  centered  around 
zero  by  subtracting  the  mean  of  the  signal. 

quency  in  both  directions  is  fn  =  5 Hz  in  both  axes.  In  water,  the  natural  frequencies 
are  fnx,w  =  3Hz  in-line  and  fny,w  =  2.2 Hz  in  crossflow. 

3.3.3  Coupling 

Excessive  friction  at  the  bushings’  contact  points  can  lead  to  a  coupled  deflection 
between  the  x  and  y  axes.  In  the  first  fabrication  attempt  of  this  design,  tolerance 
in  the  machining  of  the  pin/bushing  assemblies  was  not  tight  enough.  As  the  model 
deflected  solely  in  the  x  axis,  the  y  axis  also  bent,  and  vice  versa.  Pins  with  a  tighter 
tolerance  were  fabricated  and  found  to  cause  much  less  friction  in  the  system. 

To  carefully  test  whether  this  coupling  was  indeed  a  result  of  friction,  screws  were 
placed  through  the  bushings  to  completely  restrict  their  rotation.  Figure  3-17  shows 
a  picture  of  the  assembly  in  its  normal  (low-friction)  configuration  and  then  with 
screws  inserted  through  the  bushings. 

Pluck  tests  in  water  display  the  increased  coupling  when  the  bushings  are  re¬ 
stricted  in  rotation.  Figure  3-18  shows  the  power  density  (normalized  by  the  maxi¬ 
mum)  for  the  two  configurations.  In  the  normal  configuration,  an  excitation  in  the  y 
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Figure  3-17:  Picture  of  sensor  assembly  in  the  (left)  normal  configuration  and  (right) 
with  the  addition  of  set  screws  (circled  in  red)  to  restrict  bushing  rotation.  A  fourth 
set  screw  through  the  remaining  bushing  is  not  visible  from  this  view. 


direction  causes  an  oscillatory  response  almost  entirely  in  that  same  direction.  Min¬ 
imal  response,  and  therefore  minimal  power,  is  seen  in  the  x  signal.  On  the  other 
hand,  when  screws  are  placed  in  the  bushings,  an  excitation  purely  in  the  y  direction 
can  lead  to  significant  response  in  both  directions.  In  this  case,  the  power  in  x  is  even 
larger  than  that  in  y. 


P  luck  Test  in  Y  Pluckiest  in  Y 


Figure  3-18:  Power  density  of  signals  from  ^-direction  pluck  tests  in  water  for  the 
(left)  normal  configuration  and  (right)  configuration  with  no  bushing  rotation. 


The  VIV  characteristics  of  the  sensor  in  both  configurations  were  measured  using 
towing  tank  tests.  The  sensor  was  mounted  vertically  above  the  water  line  with  70% 
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of  the  whisker  section  protruding  into  the  water.  The  sensor  was  towed  at  several 
reduced  velocities  (Ur  =  ,  -  ,). 

v  Jny,wd/ 

The  deflection  amplitudes  in  the  crossflow  direction  are  plotted  in  Figure  3-19  for 
the  case  when  the  whisker  is  towed  with  no  coupling  and  when  it  is  towed  with  the 
set  screw  configuration  that  induces  coupling. 


Figure  3-19:  Vibrational  amplitude  in  both  axes  when  (left)  no  set  screws  were  used 
and  (right)  set  screws  were  used  to  restrict  bushing  rotation.  Ur  =  17. 


3.4  Conclusion 

Two  design  iterations  of  a  bio-inspired  flow  sensor  have  been  designed,  fabricated,  and 
characterized.  In  both  cases,  the  sensing  principle  is  based  on  deflection  measurements 
at  a  flexible  base,  similar  to  how  real  seals  collect  information  in  their  cheek  based 
on  the  whisker  motion.  Version  1  serves  as  a  proof-of-concept  that  the  whisker  VIV 
properties  can  be  realized  in  a  sensor.  Towing  tank  tests  and  field  tests  demonstrate 
its  ability  to  operate  as  a  low-cost,  low-power  flow  velocity  and  direction  sensor. 
Version  2  is  created  as  a  platform  for  systematic  studying  of  flow-induced  vibration 
properties,  with  the  goal  of  developing  wake  detection  capabilities. 
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3.4.1  Version  1 


The  first  version  shows  the  sensor’s  ability  to  use  drag-based  measurements  to  deter¬ 
mine  flow  velocity.  Readings  from  Bend  Sensors  connected  to  the  base  of  the  whisker 
are  calibrated  to  provide  flow  velocity  based  on  the  deflection  in  the  drag  direction. 
The  sensor  can  be  tuned  to  operate  optimally  for  a  desired  range  of  flow  speeds, 
simply  by  replacing  the  silicone  membrane  with  one  of  a  different  rigidity.  Also,  it 
captures  a  wide  range  of  vibrational  response  with  flow  angle,  as  predicted  in  Chapter 
2.  The  Bend  Sensor  readings  can  therefore  also  be  used  to  detect  the  direction  of 
oncoming  flow. 

The  sensor  was  designed  such  that  flexible  motion  is  possible  even  while  a  water¬ 
tight  seal  is  maintained.  Field  tests  confirm  that  the  device  can  be  used  at  depth. 
This  indicates  the  sensor’s  ability  to  be  used  for  measuring  currents  along  the  water 
column  or  while  mounted  on  an  AUV.  The  data  logging  is  done  automatically  once 
the  system  is  activated,  allowing  this  sensor  to  be  used  in  a  stand-alone  configuration. 

This  device  presents  itself  as  a  low-power  alternative  to  standard  velocity  sensors, 
such  as  DVLs,  as  it  does  not  require  active  emission  of  signals  in  order  to  conduct 
detection.  Also,  those  devices  may  cost  tens  of  thousands  of  dollars,  whereas  this 
costs  on  the  order  of  hundreds  of  dollars  or  less. 

3.4.2  Version  2 

The  second  version  was  designed  with  the  goal  of  creating  a  system  that  can  conduct 
vortex-wake  detection.  Existing  sensors  either  cannot  conduct  passive  detection  or 
else  can  not  do  so  without  experiencing  significant  background  “noise”  while  moving 
at  high  speeds.  A  sensor  designed  equivalently,  but  with  a  circular  rod  in  place  of 
the  harbor  seal  whisker  geometry,  would  have  difficulty  distinguishing  self-induced 
vibrations  from  the  vibrations  induced  by  oncoming  wake  features. 

This  device  was  therefore  designed  to  again  capture  the  VIV  reduction  properties 
of  the  whisker  geometry  and  now  to  also  allow  for  control  of  important  parameters 
for  wake- induced  vibration  studies  (i.e.  mass-damping  and  natural  frequencies).  The 
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improved  design  is  a  flexure-based  system  with  strain  gauges.  Simply  changing  the 
thickness  of  the  flexure  allows  for  variation  of  the  aforementioned  properties. 

Importantly,  contributions  from  the  in-line  and  crossflow  directions  of  motion  are 
now  separated.  The  two  axes  are  still  physically  connected,  however,  which  can  lead 
to  coupling  if  excessive  friction  is  present.  That  coupling  can  lead  to  a  dramatic 
increase  in  the  VIV  response.  The  reason  for  this  may  be  that  the  coupling  attempts 
to  induce  a  torsional  motion,  which  is  what  was  observed  in  a  separate  experiment 
conducted  with  a  flexible  whisker  model  (Appendix  C).  A  torsional  response  may  be 
a  result  of  the  elliptical  cross  section  [137]. 

This  new  system  doubles  as  a  free  vibration  experiment  setup,  on  which  detailed 
wake- induced  vibration  studies  are  conducted  in  Chapter  4.  Those  experiments  are 
used  to  study  and  develop  the  wake  detection  capabilities  of  the  sensor. 

3.4.3  Future  Work 

Room  for  improvement  remains  in  the  design,  particularly  in  terms  of  size  reduction 
and  increased  sensitivity.  The  designs  presented  here  are  intended  to  serve  as  the 
groundwork  for  what  can  be  done  with  seal-whiskerlike  sensors.  Future  integration 
into  an  AUV  sensor  suite  will  require  miniaturization  of  the  device  and  a  waterproof 
design.  Further  testing  of  alternative  sensing  devices  is  also  merited.  Optical  fibers, 
for  example,  could  prove  superior  to  strain  gauges  in  this  system. 

For  future  work  with  this  device  as  a  velocity  and  direction  sensor  (as  in  Version 
1),  it  is  recommended  that  a  row  of  whiskers  be  used.  Each  whisker  should  be  oriented 
at  a  different  angle  to  the  oncoming  flow,  allowing  for  accurate  quantification  of  the 
velocity  and  direction. 

For  future  work  with  this  device  as  a  wake  sensor,  please  see  Chapter  4.  The  goal 
is  for  this  device  to  enable  targeted  sensing  for  AUVs.  The  ability  to  detect  and  track 
a  vortex  wake  would  open  doors  for  the  oceanographic  community  to  use  robotic 
systems  to  home  in  on  hydrothermal  vents  or  enable  long-term  behavioral  studies  by 
following  animal  schools  for  extended  periods  of  time. 
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Chapter  4 


Whisker  in  Vortex  Wake:  Free 
Vibrations 


What  is  the  seals 7  mechanism  of  wake  detection? 

4.1  Motivation 

In  [58]  and  in  Chapter  2  of  this  thesis,  flow  simulations  showed  that  the  whisker 
geometry’s  undulations  break  the  spatial  coherence  of  vortex  formation  and  form 
streamwise  vorticity,  hence  reducing  substantially  the  unsteady  fluid  forces  on  the 
whisker.  This  result  implies  that  the  seal  can  move  forward  at  steady  speed  with 
small  self-induced  noise  from  VIV.  The  question  remains  on  how  the  vibrissae  detect 
unsteady  velocities,  especially  the  minute  velocities  found  in  wakes  30sec  after  a  body 
has  passed  [138]. 

To  answer  this  question  we  turn  to  the  problem  of  wake-induced  vibrations  (WIV), 
which  occur  when  a  flexibly  mounted  body  is  placed  in  the  wake  of  an  upstream  body 
so  it  is  subject  to  the  unsteady  forcing  of  its  vortical  structures.  Most  research  has 
been  conducted  in  studying  two  interacting  circular  cylinders  of  equal  or  comparable 
diameter  [65,  102,  5,  4]  and  it  is  shown  that  the  upstream  cylinder  has  significant  effect 
on  the  amplitude  and  frequency  of  motion  of  the  downstream  cylinder  even  when  it  is 
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placed  at  25  diameters  apart.  It  was  also  found  that  the  vortical  wake  of  the  upstream 
cylinder  causes  a  restoring  force  on  the  downstream  cylinder,  which  is  equivalent  to  a 
linear  spring,  hence  affecting  its  natural  frequency  [3].  In  the  few  studies  conducted 
with  interfering  cylinders  of  unequal  diameter  [115,  82,  40,  41,  39],  it  was  found  that 
the  upstream  cylinder  also  affects  the  response  of  the  downstream  cylinder,  reducing 
its  amplitude  relative  to  the  single  cylinder  case,  while  the  frequency  of  oscillation  of 
the  downstream  cylinder  was  found  to  be  either  equal  to,  or  twice  the  value  of  the 
frequency  of  the  upstream  cylinder.  Only  one  published  study  [91]  has  begun  to  show 
the  harbor  seal  whiskers’  response  under  WIV  conditions. 

In  this  chapter  we  study  a  flexibly  mounted  model  of  a  seal  whisker  placed  within 
the  wake  of  a  larger  upstream  circular  cylinder.  Strain  measurements  and  dye  visu¬ 
alizations  are  used  to  elucidate  the  whisker’s  mechanism  of  detecting  the  features  of 
the  unsteady  flow.  The  effects  of  the  separation  distance  and  size  difference  between 
the  two  objects  are  studied.  Experiments  where  the  whisker  model  is  replaced  first  by 
a  circular  cylinder  and  then  by  an  elliptical  cylinder  are  also  conducted  to  compare 
their  responses. 


4.2  Methodologies 

The  whisker  model  was  mounted  in  a  vertical  position  in  the  MIT  Towing  Tank, 
fixed  at  its  top  on  a  spring-like  structure,  which  allowed  it  to  vibrate  in  crossflow 
and  in-line  vibrations,  as  explained  in  detail  in  the  next  section.  The  models’  VIV 
responses  were  first  tested  by  towing  them  in  calm  water.  Next,  a  hydrodynamic 
wake  was  generated  upstream  of  the  whisker  model  through  the  use  of  a  vertically 
mounted  cylinder,  which  sheds  a  von  Karman  vortex  street  (Figure  4-1). 

In  the  case  of  the  live  animal,  the  whisker  diameter,  which  is  typically  about  0.5 
to  1  mm,  is  expected  to  be  significantly  smaller  than  the  dimension  of  the  dominant 
wake  vortices,  which  scale  with  the  dimension  of  the  upstream  body.  The  diameter 
of  the  upstream  cylinder  was,  therefore,  chosen  to  be  larger  than  that  of  the  whisker 
model  ( dcyi  =  {2.5,4, 11  }dw).  The  two  objects  were  towed  simultaneously  along  the 
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y  dcy|  =  {2.5,4,ll}dw 


dx  =  {8,17, 

top  view  side  view 

Figure  4-1:  (Left)  Top  view;  and  (Right)  side  view  of  the  experimental  setup.  A 
hydrodynamic  trail  was  generated  in  front  of  the  whisker  through  the  use  of  a  cylinder, 
which  shed  a  Karman  vortex  street.  Three  different  downstream  distances  ( dx )  were 
tested  for  each  of  the  three  models,  ranging  from  the  near  wake  to  the  far  wake. 
Three  different  cylinder  sizes  (dcyi)  were  also  tested.  Simultaneous  measurement 
of  the  oncoming  wake  frequency  and  whisker  response  was  done  through  use  of  a 
pressure  sensor  placed  near  the  edge  of  the  wake.  A  mariotte  bottle  supplied  dye  at 
a  constant  flow  rate.  A  camera  underneath  the  model  captured  a  bottom-up  view  of 
the  wake-whisker  interaction. 


tank  at  velocity  U,  and  the  whisker’s  amplitude  and  frequency  of  vibration  were 
measured.  The  value  of  U  covered  a  range  from  0.1  to  1.3 m/s. 

To  assess  the  effect  of  the  distance  between  the  two  objects,  the  tests  were  con¬ 
ducted  at  different  downstream  distances  ( dx  =  {8, 17, 110, 160}dw),  representing 
cases  of  close  proximity,  moderately  large  distance,  and  large  distance.  These  tests 
were  repeated  with  the  whisker  model  replaced  by  a  smooth  elliptical  (2:1  diameter 
ratio)  cylinder  model  and  then  a  circular  cylinder  model. 


4.2.1  Note  on  the  geometry 

Two  versions  of  the  whisker  geometry  are  tested  in  the  present  study.  The  geometry 
of  the  first  whisker  was  taken  directly  from  that  provided  in  [58]  (as  in  Chapters 
2  and  3)  and  scaled  up  30x.  It  is  referred  to  here  as  “whisker  #1”  and  it  has  an 
average  diameter  in  the  crossflow  direction  dw  =  1.59  cm.  A  second  whisker  geometry 
was  designed  (“whisker  #2”).  It  more  closely  matches  the  shape  measured  on  real 
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whisker  #1  whisker  #2 


Figure  4-2:  Two  variations  in  the  harbor  seal  whisker  geometry  are  captured  through 
the  depicted  parameters,  which  are  shown  here  at  the  scale  of  the  real  seal.  The 
inclination  angles  are  17.6°  and  15.27°  for  the  top  and  bottom  cross-sections  shown, 
respectively.  The  geometry  from  Hanke  et  al  (2010)  (“whisker  #1”)  has  a  A / Dm  ~  2, 
whereas  the  new  geometry  (“whisker  #2”)  has  a  A / Dm  ~  5.  The  models  used  in  this 
chapter  are  scaled  up  from  the  dimensions  shown  here  30  and  20  times,  respectively. 


harbor  seal  whisker  specimens  both  that  the  authors  received  from  the  New  England 
Aquarium  and  as  described  in  [45] .  This  model  was  enlarged  20x  from  the  real  whisker 
scale,  providing  it  a  diameter  dw  =  1.06cm.  The  only  nondimensional  difference 
between  the  two  geometries  is  the  wavelength-to-diameter  ratio:  A / Dm  «  2  in  the 
first  and  A / Dm  ~  5  in  the  second,  where  Dm  is  the  mean  overall  diameter. 

4.2.2  Strain  Measurements:  Free  Vibrations 

Apparatus  (Flow  Sensor  Version  2) 

As  detailed  in  Chapter  3,  a  whisker  sensor  was  designed  to  freely  vibrate  by  mounting 
the  base  of  the  model  on  a  flexing  plate  (1095  spring  steel,  thickness  =  0.38 mm)  with 
four  arms,  which  allows  the  model  to  freely  vibrate  in  the  in-line  (a:)  and  crossflow  (y) 
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Figure  4-3:  The  model  whisker  (“whisker  #1”  pictured  here)  was  mounted  on  a  set 
of  flexing  plates  that  allow  free  vibration  in  in-line  and  crossflow  directions.  Bushings 
are  used  to  reduce  coupling  between  the  axes. 

directions  (Figure  4-3).  It  was  therefore  in  a  cantilevered  configuration,  resembling 
that  on  the  seal,  although  the  main  body  of  the  whisker  remains  rigid  in  the  exper¬ 
iment.  This  same  device  is  used  throughout  this  chapter  to  conduct  free  vibration 
experiments. 

Strain  gauges  (Omega,  Model  KFH-03-120)  on  the  flexing  plates  were  arranged  in 
a  Wheatstone  bridge  to  output  a  signal  proportional  to  the  deflection  of  the  plates. 
Strain  signals  were  converted  to  a  voltage  once  acquired  through  a  PhidgetBridge 
1046  USB  Input  Board  (sampling  frequency  =  60 Hz,  gain  =  16).  The  measurements 
collected  were  calibrated  to  the  amount  of  deflection  at  the  tip  of  the  model.  Calibra¬ 
tion  curves  relating  voltage  and  tip  deflection  were  found  to  be  linear  for  both  axes 
on  all  models,  with  the  lowest  coefficient  of  determination  ( R 2)  value  being  0.953. 

The  natural  frequency  (fn)  of  the  models  was  found  by  conducting  pluck  tests 
in  water  and  calculating  the  zero-crossing  frequency  of  the  signals.  The  whisker  #  1 
and  elliptical  cylinder  models  have  a  natural  frequency  of  fn>x  =  3.1  Hz  in-line  and 
fn^y  =  2.2 Hz  in  crossflow.  The  whisker  #2  model  was  designed  to  have  similar  mass  to 
these  two,  allowing  it  to  have  similar  natural  frequencies  (fnjX  =  3Hz,  /n>2/  =  1 ,8Hz). 
Keeping  the  mass  of  the  two  whisker  models  similar  (for  ease  of  experimental  setup) 
meant  that  they  needed  to  be  at  different  scales  in  order  to  maintain  the  same  number 
of  undulations  on  each  model.  This  is  the  reason  that  the  whisker  #1  model  was  scaled 
up  30  times  whereas  the  whisker  #2  model  was  scaled  up  20  times.  The  circular 
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cylinder  model  has  a  smaller  mass  relative  to  the  other  models,  giving  it  a  lower 
frequency  (fn  =  1.8 Hz),  which  also  necessitated  the  use  of  a  slightly  thinner  flexure 
(1095  spring  steel,  thickness  =  0.2 mm)  to  cover  an  appropriate  reduced  velocity 
(U*  =  U/(fndw ))  range. 

All  four  models  have  a  mass  ratio  (to*  =  rn/rn(j)  of  1.4,  where  to  is  the  model 
mass  and  rn(j:  =  pir^f-span  is  the  displaced  fluid  mass.  The  whisker  #2  model  has 
a  submerged  span  of  L/dw  =  30  and  the  other  three  models  have  L/dw  =  12.  The 
Reynolds  number  based  on  the  diameter,  Re  =  Udw/v ,  ranges  from  1,060  to  20,670, 
so  it  is  generally  larger  than  the  Reynolds  number  of  seal  whiskers  which  is  around 
1,000.  The  higher  Reynolds  number  in  the  experiments  allows  easier  measurement  of 
forces  and  motions. 

4.2.3  Dye  Visualizations:  Free  Vibrations 

Apparatus 

A  mariotte  bottle  setup  was  constructed  in  order  to  introduce  dye  into  the  flow  for 
visualizations.  This  system  is  made  of  a  sealed  container  of  dye  with  an  air  inlet 
on  the  top  and  a  dye  outlet  on  the  bottom  (Figure  4-1).  The  air  inlet  ensures  that 
the  pressure  inside  remains  constant  even  as  the  reservoir  level  lowers.  The  dye  can 
therefore  be  delivered  at  a  constant  rate.  The  tracer  was  made  by  combining  8  parts 
water,  5  parts  industrial  food  coloring,  and  1  part  alcohol,  which  results  in  a  neutrally 
buoyant  liquid. 

The  dye  was  introduced  into  the  flow  between  the  wake  generator  and  the  whisker. 
A  tube  from  the  mariotte  bottle  extends  into  the  wake  generator  (hollow  PVC  pipe) 
and  dye  exits  from  a  small  hole  on  the  downstream  side  of  the  cylinder.  The  hole 
aligns  with  the  65%  span  level  of  the  whisker.  An  underwater  camera  captured  videos 
of  the  wake-whisker  interaction.  A  GoPro  Hero  3  Black  Edition  was  mounted  below 
the  whisker,  which  provided  a  bottom- up  view.  It  was  positioned  sufficiently  below 
and  downstream  of  the  model  to  avoid  disturbing  the  flow  seen  by  the  model.  Videos 
of  the  models’  wakes  were  captured  with  720p  resolution  and  at  120  frames  per  second. 
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Videos  of  the  cylinder  wake-whisker  interaction  were  captured  at  1440p  resolution  and 
at  48  frames  per  second. 

4.2.4  Wake  Frequency  Measurements 

Finally,  to  compare  the  whisker’s  response  to  the  specific  wake  being  generated,  infor¬ 
mation  on  the  upstream  wake  was  measured  in  real  time  during  the  towing  tests.  A 
pressure  sensor  (Freescale  MPXV7007,  on-chip  conditioning)  was  placed  at  the  same 
downstream  location  as  the  whisker  (Figure  4-1).  To  avoid  interfering  with  the  flow 
encountered  by  the  whisker,  the  pressure  sensor  was  placed  at  a  lateral  location  closer 
to  the  edge  of  the  wake  (at  an  angle  of  /3  =  60  —  75°  between  the  upstream  cylinder 
and  the  whisker) .  Plastic  tubing  was  connected  to  the  sensor  inlet  and  placed  into  the 
wake,  such  that  the  sensor  avoided  direct  contact  with  water.  To  reduce  flow-induced 
vibrations  of  the  tube  itself,  it  was  attached  to  the  leading  edge  of  a  small  airfoil  wing, 
which  was  rigidly  mounted  to  the  carriage.  A  National  Instruments  USB-6211DAQ 
board  was  used  to  collect  the  analog  signal  (sampling  frequency  =  1000772;).  The 
vortex  formation  frequency  of  the  wake  was  extracted  from  the  pressure  signal. 

4.2.5  Analysis 

The  strain  signals  were  recorded  as  function  of  time  for  each  run.  A  bandpass  Butter- 
worth  filter  (pass  range:  0.4  to  20Hz)  was  applied.  To  exclude  transient  phenomena, 
only  a  middle  segment  of  these  signals  was  considered  further.  The  vibration  ampli¬ 
tude  for  that  run  was  determined  by  taking  the  root-mean-square  (RMS)  amplitude 
of  that  segment.  The  value  was  then  nondimensionalized  by  dw.  The  crossflow  fre¬ 
quency  of  vibration  was  determined  by  taking  a  fast  Fourier  transform  (FFT)  of 
the  signal,  selecting  the  peak  response  frequency,  and  nondimensionalizing  it  by  the 
model’s  natural  frequency  in  water:  f*  =  fy/fn,y 

For  the  experiments  which  utilize  simultaneous  measurements  of  wake  pressure 
and  strain  on  the  whisker  sensor,  the  response  was  analyzed  over  a  more  fine  time 
scale.  First,  the  two  signals  were  re-sampled  to  the  same  frequency  and  the  following 
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signal  conditioning  was  applied  to  both  signals.  A  bandpass  Butterworth  filter  was 
used  to  pass  frequencies  between  1/2  and  2x  the  expected  wake  frequency,  based 
on  equation  1.2.  Additional  Butterworth  filters  were  used  to  notch  out  the  natu¬ 
ral  frequencies  of  the  foil  (~  1, 10,20 Hz).  Cross-correlation  of  the  two  signals  was 
conducted  to  quantify  any  similarity. 

This  correlation  was  first  conducted  by  dividing  the  signals  by  their  RMS  value 
in  order  to  normalize  their  magnitudes.  An  initial  cross-correlation  was  carried  out 
to  determine  the  value  of  any  time  delay  between  the  two  signals.  The  signals  were 
shifted  by  a  delay  corresponding  to  the  peak  in  this  initial  cross-correlation  result. 
A  second  cross-correlation  was  then  used  to  determine  the  level  of  synchronization 
between  the  two  signals.  This  value  was  normalized  by  the  magnitude  to  reach  a 
cross-correlation  coefficient  (c),  which  is  the  variable  presented  in  the  following  plots. 


4.3  Results 

4.3.1  VIV  response 

First,  the  vortex-induced  vibrations  of  each  model  were  measured  as  the  models  were 
towed  through  calm  water.  The  whiskers  and  elliptical  cylinder  were  oriented  with 
their  long  axis  aligned  with  the  flow  direction.  Figure  4-4  shows  the  nondimensional 
tip  deflection  for  all  three  models  in  the  crossflow  direction  (A*  =  Ay/dw)  and  in 
the  in-line  direction  (A*  =  Ax/dw)  versus  the  nondimensional  towing  velocity  (U*  = 
U/(fn,ydw )).  The  nondimensional  frequency  of  crossflow  oscillation  (/*  =  fy/fn,y)  is 
also  plotted. 

[37]  tested  a  cylinder  in  a  similar  configuration  as  in  our  experiments,  viz.  pivoted 
at  the  top,  freely  vibrating  in  the  in-line  and  crossflow  directions  and  found  a  peak 
amplitude  of  vibration  of  1.5 d  at  the  tip.  That  value  is  close  to  what  is  seen  for  the 
cylinder  here,  providing  some  validation  for  this  new  setup. 

The  same  trends  seen  in  forced  vibration  experiments  (Chapter  2)  are  found  here: 
The  whisker  at  0°  experiences  low  VIV  compared  to  a  circular  cylinder.  The  oscil- 
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Figure  4-4:  VI V  response  of  each  model  positioned  in  the  streamlined  direction  (a  = 
0°).  Bars  depict  data  spread  over  three  repeat  runs.  (Left)  Both  cross-flow  and 
in-line  amplitudes  of  the  whisker  and  elliptical  cylinder  are  significantly  smaller  than 
those  of  the  circular  cylinder.  The  frequencies  at  which  the  whisker  and  ellipse  vibrate 
remain  close  to  their  natural  frequency.  The  diagonal  dotted  line  indicates  a  frequency 
that  matches  the  Strouhal  relation.  (Right)  Both  whiskers  and  the  elliptical  cylinder 
vibrate  with  A*  <  0.2.  Whisker  #2  has  the  lowest  VIV  amplitude. 


lation  frequency  for  the  elliptical  cylinder  and  whisker  remain  close  to  their  natural 
frequency  across  all  speeds  tested.  It  should  be  noted  that  in  previous  seal  whisker 
vibration  studies  ([58]  and  Chapter  2),  the  additional  degree-of-freedom  in  the  in-line 
direction  was  not  included.  Motion  in  the  in-line  axis  can  induce  additional  vibrations 
[17].  In  spite  of  this,  the  whisker’s  VIV  is  seen  to  remain  low. 

Figure  4-4(b)  shows  a  zoomed  in  view  of  the  cross-flow  response.  The  elliptical 
cylinder  exhibits  similar,  slightly  higher  vibrations  than  the  whisker.  The  whisker 
#2  model  has  an  even  lower  response  than  the  original  model.  These  two  geometry 
variations  may  align  with  the  two  minimums  in  response  seen  in  the  Lam  &  Lin 
(2009). 

Dye  visualizations  of  the  0°  configuration  illuminate  wake  patterns  that  explain 
the  recorded  deflection  measurements  (Figure  4-5).  A  circular  cylinder  sheds  large 
coherent  vortices  in  a  2S  pattern.  The  elliptical  cylinder  also  sheds  coherent  vortices 
but  with  a  narrower  wake.  The  whisker  #1  model  displays  a  less  coherent  vortex 
wake  that  forms  farther  downstream.  The  whisker  #2  model  creates  a  similar  wake 
(Figure  4-6  ). 
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Figure  4-5:  Dye  visualizations  of  the  wake  behind  a  circular  cylinder,  elliptical  cylin¬ 
der,  and  whisker  (from  left  to  right)  at  a  =  0°.  The  first  two  objects  shed  regular 
2S  patterns,  whereas  the  whisker  #1  model  sheds  an  incoherent  wake.  Images  were 
taken  during  runs  corresponding  to  each  model’s  peak  amplitude:  cylinder  at  U*  =  10, 
ellipse  at  U*  —  28,  and  whisker  #1  at  U*  =  28. 


Figure  4-6:  Dye  visualization  of  the  wake  behind  the  whisker  A  2  model.  Similar  to 
the  whisker  #1  model,  it  sheds  an  incoherent  wake.  The  image  was  taken  during  the 
run  corresponding  to  the  model’s  peak  amplitude,  which  is  at  U*  =  32. 


Because  of  the  slightly  superior  VIV  reduction  performance  of  the  second  whisker 
model,  the  remainder  of  this  paper  focuses  on  it.  If  the  model  number  is  not  explicitly 
mentioned,  it  should  be  taken  to  refer  to  the  whisker  #2  model. 


4.3.2  Whisker  response  within  a  wake  ( dcyi  =  2.5 dw,dx  =  8 dw) 

We  then  consider  the  case  where  the  whisker  is  at  a  relatively  close  distance  (dx  =  8 dw) 
from  the  upstream  cylinder  {dcyi  =  2.5 dw).  The  whisker  is  facing  the  flow  in  its 
streamlined  direction  (a  =  0°).  The  whisker’s  nondimensional  deflection  at  the  tip 
is  shown  in  Figure  4-7(a).  The  amplitude  of  vibration  in  this  case  peaks  at  >  1.5 dw. 
This  value  has  to  be  contrasted  with  the  whisker’s  vibration  amplitude  as  it  moves 
through  open  water  (depicted  with  the  “x”  ’s) ,  which  is  significantly  lower  than  that 
-  peaking  at  ~  0.05 dw.  Hence,  for  this  distance  from  ( dx )  and  size  ratio  with  the 
upstream  cylinder  (dcyi/dw),  the  whisker  can  vibrate  with  30x  larger  amplitude  when 
it  encounters  the  wake  than  it  does  in  open  water.  As  above,  the  symbols  mark  the 
average  of  three  repeat  runs  and  the  bars  denote  the  data  spread. 

Also  of  note  is  the  comparison  of  the  whisker’s  frequency  of  response  between  the 
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two  cases,  shown  in  Figure  4-7(b),  where  the  nondimensional  frequency  of  vibration 
(/*)  is  plotted  versus  U*.  The  frequencies  plotted  here  were  determined  by  selecting 
the  peak  in  the  FFT  of  the  corresponding  signal’s  time  trace.  As  the  whisker  moves 
through  open  water,  it  oscillates  at  a  relatively  constant  frequency.  However,  when 
the  whisker  encounters  the  oncoming  vortical  wake,  it  oscillates  at  the  vortex  shedding 
frequency  of  the  upstream  cylinder  (/„).  This  frequency  is  measured  with  a  pressure 
sensor  during  the  same  runs,  nondimensionalized  by  fnyj,  and  plotted  with  the  red 
bars  (again  showing  the  data  spread  over  three  runs) .  Hence  the  whisker  synchronizes 
with  the  oncoming  wake  dominant  frequency,  which  is  distinct  from  its  VIV  frequency 
in  open  water. 

These  are  two  strong  indicators  for  wake  detection  that  the  seal  can  potentially 
use:  The  whiskers  which  normally  vibrate  at  low  amplitude  when  in  open  water, 
vibrate  with  more  than  an  order-of-magnitude  larger  amplitude  when  in  the  wake  of 
a  body;  more  importantly,  the  frequency  of  vibration  synchronizes  with  the  dominant 
frequency  in  the  wake. 


Figure  4-7:  For  a  whisker  within  the  wake  of  an  upstream  cylinder  {dcyi  =  2.5 dw) 
at  a  distance  dx  =  8 dw:  (a)  The  whisker  vibrates  with  far  higher  amplitude  in  the 
wake  than  it  does  on  its  own.  (b)  The  peak  frequency  of  the  whisker’s  vibration 
synchronizes  with  the  oncoming  wake  frequency.  The  diamonds  and  “x”s  are  the 
average  over  three  repeat  runs.  The  bars  indicate  the  data  spread  of  the  measured 
wake  frequency  over  the  same  three  runs. 
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Figure  4-8:  Schematic  of  the  camera  perspective  of  the  dye  runs.  The  bottom-up 
view  shows  vortices  being  shed  by  the  wake  generator  on  the  right  and  impinging  on 
a  cross-section  of  the  whisker  on  the  left. 

4.3.3  Dye  visualizations  in  wake 

Flow  visualizations  of  the  whisker  interacting  with  the  wake  confirms  and  explains 
these  two  findings.  Similar  parameters  displayed  above  are  used  here:  dcyi  =  2.5 dw,  dx  = 
17 dw.  Videos  are  taken  with  a  bottom-up  perspective  depicted  in  Figure  4-8.  Dye 
exits  the  wake  generator,  off  the  right  side  of  images,  is  entrained  by  the  vortices  shed 
by  the  cylinder  and  then  impinges  on  the  whisker.  The  whisker  is  painted  white,  and 
the  cross-section  which  aligns  with  the  dye  plane  is  outlined.  A  sequence  of  snapshots 
from  the  video  are  shown  in  Figure  4-9. 

The  motion  of  the  whisker  is  driven  by  the  alternating  pressure  field  of  the  wake, 
forcing  it  to  oscillate  at  the  wake  frequency.  It  is  remarkable  that  the  trajectory  of 
the  whisker  resembles  the  slaloming  motion  recorded  for  trout  swimming  behind  a 
circular  cylinder,  which  extracts  energy  from  the  oncoming  von  Karman  vortices  by 
also  synchronizing  with  the  frequency  of  the  vortices  [85] . 

4.3.4  Cross-correlation  analysis 

More  detailed  analysis  of  the  response  between  the  whisker  and  the  upstream  wake 
was  conducted  by  cross-correlating  the  strain  signal  at  the  whisker  base  with  the 
pressure  signal  in  the  wake.  A  cross-correlation  value  of  1  indicates  perfect  synchro¬ 
nization  between  the  whisker  and  the  wake  frequencies,  and  a  value  of  0  indicates  no 
synchronization. 

Figure  4-10  shows  a  sample  result  from  this  analysis.  On  top  are  the  two  signals 
measured  over  the  course  of  one  run.  The  signals  were  normalized  and  aligned  in 
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Figure  4-9:  Sequence  of  dye  visualizations  (left)  and  schematics  (right)  showing  the 
whisker’s  free  vibration  path  through  the  vortex  wake.  The  frequency  of  the  whisker’s 
motion  aligns  with  the  vortex  wake  frequency,  as  the  whisker  slaloms  between  the 
vortices.  (U*  =  16,  Re  =  3, 180) 
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Figure  4-10:  Sample  result  of  the  cross-correlation  analysis.  (Top)  The  normalized 
signals  of  strain  and  pressure  measured  on  the  whisker  and  in  the  wake,  respectively. 
(Bottom)  The  cross-correlation  coefficient  at  each  oscillation  cycle.  (dx  =  17 dw,  dcyi  = 
Adw,  U*  =  27) 

time  via  an  initial  cross-correlation.  On  bottom  is  c,  the  cross-correlation  coefficient 
for  those  two  signals  at  each  oscillation  cycle.  In  this  case,  there  are  segments  of 
poor  correlation  near  the  beginning  and  end  but  a  middle  segment  of  the  run  reaches 
strong  correlation. 

This  type  of  analysis  is  repeated  for  the  following  experiments.  The  results  are 
presented  with  a  few  different  metrics  that  aim  to  highlight  the  whisker’s  ability  to 
detect  the  wake  under  a  variety  of  conditions. 

Effect  of  cylinder  size 

To  test  the  effect  of  the  size  of  the  wake  generator  on  the  detection  ability  of  the 
whisker,  three  different  size  cylinders  were  used  ( dcyi  =  {2.5,4,  ll}dw).  The  whisker 
was  placed  at  a  downstream  distance  based  on  cylinder  diameters  dx  =  15 dcyi.  This 
is  to  ensure  that  the  wake  would  have  evolved  equivalently  in  each  case. 

In  Figure  4- 11  (a),  the  maximum  cross-correlation  coefficient  per  run  is  plotted. 
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Figure  4-11:  Cross-correlation  metrics  of  the  whisker  vibration  with  the  wake  fre¬ 
quency  for  three  different  size  cylinders  at  dx  =  15 dcyi.  (Left)  All  cases  see  a  high 
maximum  cross-correlation  coefficient.  (Right)  The  percentage  of  the  run  where  the 
whisker  is  “synchronized”  generally  increases  with  cylinder  diameter. 


For  all  sizes  of  cylinders  tested,  there  is  a  high  (>  0.85)  maximum  value.  We  also 
define  a  sufficient  “synchronization”  as  when  the  instantaneous  cross-correlation  co¬ 
efficient  is  either  >  80%  or  <  —80%.  The  amount  of  time  that  the  whisker  spends 
being  “synchronized”  with  the  wake  is  seen  to  vary  with  the  size  of  the  cylinder,  and 
therefore  the  size  of  the  dominant  wake  features  (Figure  4- 11(b)). 

In  Figure  4-12  the  percent  of  the  run  that  is  “synchronized”  (as  averaged  over  all 
towing  speeds  tested)  is  seen  to  decrease  with  separation  distance  and  lateral  distance 
away  from  the  centerline.  This  synchronization,  however,  remains  present  in  all  cases 
shown  here,  indicating  that  the  seal  could  detect  the  presence  of  the  wake  from  a 
distance  downstream  and  also  near  the  edge  of  the  wake. 


4.3.5  Comparison  between  the  whisker  and  the  cylinder  and 
elliptical  geometries 

Experiments  were  carried  out  using  the  whisker,  elliptical  and  circular  cylinders 
at  three  downstream  distances  the  first  closer  to  the  upstream  cylinder  ( dcyi  = 
2.5dw,dx  =  8 dw),  then  moderately  far  away  ( dx  =  17 dw)  and  finally  in  the  far  wake 
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Figure  4-12:  Average  of  percent  synchronization  across  all  speeds  tested  ( dcyi  =  4 dw). 
The  trends  are  displayed  versus  separation  distance  both  in  the  downstream  ( dx )  and 
lateral  ( dy )  directions.  Experiments  with  a  lateral  offset  were  conducted  only  on  one 
side  of  the  cylinder  but  the  results  are  assumed  to  be  symmetric  and  are  therefore 
mirrored  here. 

0 dx  =  110c4). 

The  frequencies  of  response  are  shown  in  Figure  4-13,  with  dx  increasing  from  left 
to  right.  At  the  closest  distance,  all  geometries  synchronize  with  the  theoretical  wake 
frequency,  based  on  St  =  0.2.  As  the  distance  increases  to  dx  =  17 dw,  the  circular 
cylinder  deviates  from  the  wake  frequency,  but  the  elliptical  cylinder  and  whisker 
remain  synchronized.  When  the  models  encounter  the  far  wake  at  dx  =  110d„,, 
they  exhibit  different  trends.  The  whisker  locks  to  the  upstream  cylinder’s  Strouhal 
frequency  over  a  narrow  frequency  range.  It  should  be  noted  that  the  far  wake 
develops  new  vortical  patterns  that  are  different  than  the  originally  shed  patterns 
close  to  the  cylinder  [141]. 

The  amplitudes  of  vibration  for  the  three  models  are  shown  in  Figure  4-14,  with 
the  circular  cylinder  on  the  left,  elliptical  cylinder  in  the  middle,  and  whisker  on 
the  right.  Results  from  the  three  different  separation  distances  are  shown,  as  well 
as  the  response  in  open  water.  Except  at  the  closest  distance,  the  circular  cylinder 
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Figure  4-13:  The  frequencies  of  response  nondimensionalized  by  the  theoretical  fre¬ 
quency.  At  the  close  distance  (8dw),  all  four  models  vibrate  at  the  theoretical  wake 
frequency.  As  the  distance  increases  to  17 dw,  the  elliptical  cylinder  and  whisker  con¬ 
tinue  to  oscillate  at  the  wake  frequency,  but  the  cylinder  deviates  from  it. 

vibrates  at  roughly  similar  amplitudes  in  the  wake  as  it  does  in  open  water.  The 
elliptical  cylinder  and  whisker,  on  the  other  hand,  oscillate  substantially  more  as 
they  encounter  the  wake.  These  amplitudes  generally  decrease  as  the  downstream 
distance  increases. 

At  the  two  closest  distances,  all  three  shapes  the  cylinder,  the  elliptical  cylinder, 
and  the  whisker  show  signs  of  wake  induced  response,  which  is  similar  to  trends 
reported  in  the  literature  for  cylinders  in  tandem  at  similar  distances  dx/d  =  2  to  5.6 
[65,  5,  4,  3], 

Unlike  the  other  two  models,  the  whisker  continues  to  vibrate  with  amplitude 
significantly  higher  than  its  VIV  response  at  the  far  distances  {110, 160} .  The  wake- 
induced  vibration  (WIV)  amplitude  of  the  whisker  at  the  largest  distance  reaches  >  7 
times  higher  than  the  VIV  in  open  water,  indicating  the  whisker’s  superior  ability  to 
detect  a  wake  from  a  far  distance.  Testing  whether  the  effect  sustains  for  the  longer 
detection  distances  (O(104)  whisker  diameters)  stated  in  [25,  118,  138]  is  not  possible 
in  this  facility. 

4.3.6  Wake  Detectability 

The  WIV  frequencies  and  amplitudes  indicate  that  the  elliptical  cylinder  performs 
about  the  same  or  even  higher  than  the  whisker.  But  a  distinction  between  the  two 
geometries  is  highlighted  through  a  detectability  index ,  which  we  define  as  the  ratio 
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Figure  4-14:  Vibration  amplitudes  of  the  circular  cylinder,  elliptical  cylinder,  and 
whisker  (from  left  to  right,  respectively).  The  latter  two  vibrate  more  when  they  are 
in  the  wake  of  an  upstream  cylinder  (colored  lines)  than  they  do  in  open  water.  The 
circular  cylinder  vibrates  with  a  similar  amplitude  for  all  cases.  The  amplitudes  of 
response  decrease  as  the  distance  from  the  upstream  cylinder  increases.  The  whisker 
retains  an  amplitude  greater  than  its  VIV  amplitude  even  at  the  largest  distance. 


A  =  mm(^ro;afce ).  This  ratio  of  the  models’  vibration  amplitude  in  the  wake  and  their 
self-induced  vibrations  at  the  same  flow  velocity  is  chosen  as  a  metric  that  highlights 
the  object’s  ability  to  detect  the  presence  of  a  wake  above  its  VIV  response.  The 
minimum  ratio  over  all  U*  cases  is  chosen  to  be  conservative.  The  value  of  that 
variable  is  plotted  in  Figure  4-15.  The  whisker’s  detectability  is  higher  than  that 
of  the  elliptical  and  circular  cylinders.  This  inequality  remains  true  over  all  of  the 
distances  tested.  This  is  due  in  part  to  the  whisker  having  a  smaller  VIV  amplitude 
than  the  other  two  models  over  most  of  the  reduced  velocity  range. 


4.4  Discussion 

The  study  of  a  model  problem,  the  interaction  of  a  harbor  seal  whisker  model  with  the 
wake  of  an  upstream  cylinder  with  large  diameter,  provides  the  basic  fluid  mechanisms 
that  live  seals  may  employ  when  they  exhibit  their  outstanding  ability  to  detect  and 
track  the  wakes  of  prey  even  with  a  30sec  late  start.  In  open  water  the  whisker 
model  exhibits  very  low  vibration  compared  to  its  width  when  its  long  axis  is  aligned 
with  the  incoming  flow.  However,  once  it  enters  the  wake  of  an  upstream  cylinder, 
it  oscillates  with  large  amplitude  on  the  order  of  its  own  width.  More  important, 
its  frequency  coincides  with  the  Strouhal  frequency  of  the  upstream  cylinder,  hence 
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Figure  4-15:  The  detectability  of  the  whisker  is  higher  than  the  other  models.  As 
the  separation  distance  increases,  the  wake  becomes  more  difficult  to  detect,  but  the 
whisker  remains  the  highest. 

making  detection  of  an  upstream  wake  and  the  size  of  the  wake-generating  body 
possible. 

This  is  a  new  result,  which  is  different  from  the  case  of  two  interacting  cylinders 
of  comparable  diameter,  wherein  the  downstream  flexibly  mounted  cylinder  oscillates 
at  a  frequency  determined  by  its  mass  and  spring  properties,  the  value  of  the  added 
mass,  and  the  wake  spring  constant  of  the  wake  defined  in  [3].  In  the  case  of  a 
stationary  upstream  cylinder  with  significantly  larger  diameter  than  the  downstream 
body  (2.5x,  4x,  and  llx  in  our  experiments),  the  flexibly  mounted  downstream  body 
vibrates  freely  at  the  Strouhal  frequency  of  the  upstream  cylinder,  following  a  path 
that  slaloms  among  the  vortices  of  the  oncoming  wake.  This  path  leads  to  significant 
energy  extraction  as  also  noted  in  the  visualizations  of  live  trout  swimming  behind 
bluff  cylinders  [85]  and  the  upstream  swimming  of  anaesthisized  fish  and  flapping  foils 

[7]- 

While  a  seal  can  swim  in  open  flow  with  minimal  flow  noise,  as  shown  in  Chapter 
2,  once  it  enters  the  vortical  flow  of  the  wake,  its  amplitude  of  response  increases 
substantially  even  at  large  distances,  while  its  frequency  locks  at  the  frequency  of  the 
incoming  vortices  enabling  it  to  detect  the  size  of  the  upstream  body. 

In  our  study  we  employed  wakes  of  bluff  bodies  to  demonstrate  the  detection 
mechanisms  up  to  a  distance  of  160  whisker  diameters  downstream  of  the  body. 
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Far  wakes  are  known  to  reorganize  their  vortical  structures  creating  new  dominant 
frequencies  [141],  which  places  limits  on  the  detectability  of  the  specific  object  shape, 
but  the  presence  of  a  distinct  wake  is  still  detectable.  The  use  of  different  sized 
cylinders  as  the  wake  generator  indicates  that  the  whisker’s  detection  ability  tends  to 
increase  with  the  size  of  the  wake  features. 

The  study  of  a  circular  cylinder  and  a  2:1  elliptical  cylinder  serves  to  underline 
the  special  properties  of  the  whisker  shape.  The  circular  cylinder  responds  with  large 
amplitude  even  in  open  flow,  hence  making  detection  difficult. 
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Chapter  5 


Conclusions 

5.1  Overview 

Behavioral  experiments  with  harbor  seals  have  demonstrated  their  outstanding  ability 
to  detect  and  track  hydrodynamic  signatures  left  by  swimming  animals  and  moving 
objects.  Despite  having  their  auditory  and  visual  sensory  cues  blocked,  the  seals 
successfully  follow  the  paths  of  bodies  that  had  swum  ahead  of  them  several  seconds 
earlier.  The  seals’  vibrissae  (whiskers)  are  reported  to  be  sensitive  enough  to  detect 
even  minute  water  movements,  such  as  those  left  in  the  hydrodynamic  trails  of  fish. 

Various  animals  utilize  hairlike  structures  as  direct  touch  sensors.  For  aquatic 
animals  that  swim  forward  at  relatively  high  speeds,  however,  it  is  rather  surprising 
that  they  would  be  able  to  use  their  whiskers  to  sense  hydrodynamic  information.  The 
cross-sectional  shape  of  the  whiskers  is  close  in  form  to  a  bluff  (non-streamlined)  body, 
whose  wake  is  known  to  spontaneously  form  a  double  array  of  staggered  vortices,  the 
Karman  street.  The  unsteady  forces  caused  by  the  Karman  street  result  in  vortex- 
induced  vibrations  (VIV)  on  the  body,  with  an  amplitude  comparable  to  its  cross- 
sectional  dimension  and  with  high  frequency,  close  to  that  predicted  by  the  Strouhal 
law,  resulting  in  high  transverse  velocity.  Such  vibrations  would  act  as  strong  noise 
for  the  sensor,  since  it  would  be  difficult  to  discriminate  between  the  wake-induced 
fluctuations  and  external  flow  unsteadiness. 

The  phenomenon  of  VIV  is  ubiquitous  for  bluff  bodies  and  very  complex,  as  the 
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body  vibrations  caused  by  the  Karman  street  alter,  in  turn,  the  vortical  structures 
in  the  wake  and  hence  the  resulting  forces  and  motions.  That  establishes  a  feed¬ 
back  mechanism  that  leads  to  a  stable  or  unstable  dynamic  equilibrium,  which  varies 
significantly  as  function  of  the  basic  parameters  of  the  problem. 

In  this  thesis,  the  harbor  seal  whisker’s  ability  to  reduce  VIV  is  quantified  and 
explained.  This  feature  is  realized  in  a  physical  device.  Lastly,  a  mechanism  of  vortex 
wake  detection  is  elucidated.  A  single  whisker  is  sufficient  to  detect  a  wake  and  is 
shown  to  do  so  passively. 


5.2  Principal  Contributions 

5.2.1  Forced  vibrations  on  a  whisker  model 

Experiments  with  a  scaled  harbor  seal  whisker  model  undergoing  forced  vibrations 
demonstrate  the  undulatory,  elliptical  geometry’s  ability  to  passively  reduce  VIV. 


Low  Civ  =  0  contour 

The  lift  coefficient  in  phase  with  velocity  {C^v)  was  measured  over  a  range  of  imposed 
frequencies  and  amplitudes.  The  zero  contour  was  found  to  be  very  low  ( A/d  <  0.1), 
indicating  a  dramatic  reduction  in  VIV  amplitude  from  that  of  a  circular  cylinder  at 
the  same  Reynolds  number  ( A/d  ~  0.7). 


Small  bandwidth  of  response 

The  added  mass  coefficient  reveals  that  the  whisker’s  bandwidth  of  excitation  is  dra¬ 
matically  reduced.  Cm  on  both  circular  and  elliptical  cylinders  deviate  from  the 
potential  flow  solution,  as  the  vortex  shedding  forces  the  body  to  lock-in  over  a  wide 
range  of  exception.  The  whisker,  however,  breaks  down  the  correlation  along  the 
model  span,  leading  to  a  narrower  range  of  excitation. 
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Decorrelated  Karman  street 


Wake  visualizations  of  the  whisker  geometry  help  explain  this  dramatic  reduction  of 
VIV.  The  shed  vortices  are  decorrelated,  resulting  in  a  jumbled  wake,  which  exerts 
low  forcing  back  onto  the  model.  A  Karman  street  does  form,  but  does  so  farther 
downstream  than  a  similar  model  with  smooth  elliptical  cross  section. 

Low  drag 

The  whisker  geometry  is  also  seen  to  exhibit  low  drag:  Cd  ~  0.7,  which  is  even  lower 
than  that  of  the  elliptical  cylinder.  Thus,  the  whisker  is  able  to  attain  both  low 
drag  and  low  vibration  amplitude  a  feat  not  achieved  simultaneously  in  most  VIV 
suppression  devices. 

Wide  range  of  response  with  angle  of  attack 

The  VIV  response  of  the  whisker  geometry  changes  significantly  with  the  angle  at 
which  it  faces  the  flow.  The  amplitude  of  response  increases  ~  20x  as  the  angle 
of  attack  (a)  increases  to  90°.  The  elliptical  cross  section  of  the  geometry  is  seen 
to  dominate  the  response  at  higher  angles,  as  the  undulatory  features  are  no  longer 
large  enough  to  have  a  strong  effect. 

This  reduces  the  usability  of  this  geometry  as  a  VIV  suppression  device  in  situa¬ 
tions  where  oncoming  flow  may  approach  from  a  range  of  directions.  It  can,  however, 
be  used  to  sensitively  detect  direction  of  steady  oncoming  flow. 

5.2.2  Design  of  whisker  sensor 

The  unique  vibration  properties  discovered  in  the  whisker  morphology  were  brought 
into  use  as  a  sensor  through  two  successive  designs. 

Flow  velocity  and  direction  sensor 

Drag-based  measurements  are  used  to  sense  oncoming  flow  velocity.  The  wide  range  of 
response  based  on  angle  of  attack  allows  this  to  also  be  used  as  a  flow  direction  sensor. 
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Figure  5-1:  Pictures  of  the  whisker-inspired  flow  velocity  and  direction  sensor. 


Figure  5-2:  Picture  of  the  re-designed  whisker  sensor,  which  was  made  for  systematic 
testing  in  a  vortex  wake. 

This  is,  however,  not  the  feature  that  enables  wake  detection,  which  is  described 
through  the  next  round  of  experiment  results. 

Do  seals  actually  use  their  vibrissae  as  velocity  sensors?  Perhaps  they  default  to 
other  sensory  systems,  as  in  [111],  in  which  seals  could  maintain  constant  speed  both 
with  and  without  presence  of  its  vibrissae.  In  that  case,  the  seals  may  have  used 
visual  cues  instead.  Despite  a  potential  preference  for  using  other  sensory  systems 
to  determine  speed,  the  discovered  whisker  vibration  properties  still  present  useful 
sensing  strategies  for  engineering  applications. 

Vortex  wake  sensor 

An  improved  design  is  put  forward  for  testing  of  vortex-wake  detection.  A  flexure- 
based  system  allows  for  free  vibration  of  the  whisker  model.  Strain  gauges  are  used 
to  measure  deflection  at  the  base. 
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Figure  5-3:  Schematic  of  the  slaloming  trajectory  (dotted  line)  that  the  whisker  takes 
as  it  encounters  a  vortex  wake. 

5.2.3  Mechanism  of  wake  detection 

Experiments  with  the  whisker  sensor  undergoing  free  vibrations  are  used  to  describe 
the  vibrational  response  of  the  whisker  as  it  encounters  a  vortex  wake.  The  whisker 
exhibits  a  dramatic  difference  in  vibrational  response  to  steady,  open  flow  versus  a 
vortex  wake. 

High  vibration  amplitude  in  wake 

The  whisker  vibrates  with  much  higher  amplitude  in  a  vortex  wake  than  it  does  in 
steady  flow.  Depending  on  the  distance  downstream  and  the  size  ratio  difference 
between  the  upstream  body  and  the  whisker,  this  difference  can  span  more  than  an 
order  of  magnitude. 

Synchronization  with  wake  frequency 

Also  important  is  the  frequency  of  the  whisker’s  response.  Upon  entering  the  vortex 
wake,  the  whisker’s  frequency  of  vibration  synchronizes  with  the  vortex-shedding 
frequency  of  the  wake  generator.  This  is  consistent  across  several  distances  and  size 
ratios  tested. 

Energy  extraction 

The  underlying  mechanism  of  the  whisker’s  response  in  the  vortex  wake  resembles 
that  seen  in  trout  slaloming  behind  rocks.  In  both  cases,  the  downstream  body  is 
able  to  move  in  a  trajectory  that  allows  for  extract  energy  from  the  oncoming  wake. 
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5.3  General  Implications 


5.3.1  Sensory  systems  employed 

It  is  unknown  what  information  real  seals  rely  on  in  actual  prey  pursuit  conditions. 
Several  studies  provide  evidence  for  sensory  information  coming  from  outside  of  the 
vibrissal  system.  For  example,  bioluminescence  produced  either  by  the  prey  or  the 
surrounding  plankton  could  conceivably  aid  a  visual  search  [63].  Also,  seals  could 
use  backlit  light  from  the  surface  to  aid  visual  sensing  [63].  In  one  experiment  [108], 
the  seals’  ability  to  catch  live  trout  did  not  appear  to  be  significantly  affected  by  the 
removal  of  vibrissae.  The  experimental  setup,  however,  was  such  that  vision  may  have 
been  sufficient.  [96]  suggests  that  echolocation  may  be  more  useful.  Off  the  coast  of 
California,  seals  were  instrumented  with  accelerometers  on  their  jaws.  They  measured 
the  number  of  jaw  motion  events  per  day  to  be  much  larger  than  the  number  of  pings 
that  toothed  whales  send  out.  So  they  suggest  that  the  acoustic  method  may  be 
more  effective  than  vibrissae.  [33]  proposes  that  the  vibrissae  on  Arctic  seals  do  serve 
a  functional  sensory  role,  as  the  animals  had  a  difficult  time  finding  ice  holes  when 
their  vibrissae  were  restricted.  Or  perhaps  the  vibrissae  can  be  used  for  orienting  the 
seals’  mouths  with  the  prey  in  the  last  few  inches,  since  the  snout  may  obstruct  the 
seals’  view  of  the  fish  at  that  point  [108]. 

Certainly  the  seals’  choice  of  sensory  system  is  situation-dependent.  And  the 
animals  likely  utilize  a  combination  of  sensory  channels  to  successfully  complete  the 
complex  behaviors  required  to  capture  prey.  The  mechanisms  elucidated  in  this  thesis 
help  explain  how  seals  can  conduct  wake  detection  from  long  range.  More  remains 
to  be  explored  in  determining  additional  ways  that  vibrissae  can  be  used,  either  as  a 
supplement  or  substitute  for  other  types  of  sensory  information. 

5.3.2  Comparison  with  other  species 

As  noted  in  [56],  ecological  factors  may  be  correlated  with  the  varying  evolution  of 
vibrissal  morphology  across  different  species.  For  example,  the  monk  seal  is  the  only 
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phocid  species  that  exclusively  inhabits  tropical  and  subtropical  waters,  meaning  that 
it  should  be  able  to  more  heavily  rely  on  vision  in  those  clearer  waters.  Therefore 
its  vibrissal  profile  perhaps  did  not  need  to  evolve  from  smooth  to  beaded.  Bearded 
seals,  which  also  have  a  smooth  profile,  are  benthic  foragers  and  therefore  likely  need 
their  whiskers  more  in  a  tactile  role.  Gray  seals  do  have  a  beaded  profile,  but  it  is 
less  pronounced.  They  are  largely  benthic  or  demersal  feeders,  so  they  also  may  not 
need  to  rely  on  their  vibrissae  for  wake  tracking. 

The  use  of  vibrissae  as  hydrodynamic  receptors  may  not  be  uncommon  at  all  for 
marine  mammals.  Only  toothed  whales  (dolphins,  sperm  whales,  etc.)  have  sonar 
systems.  All  others  have  well-developed  vibrissae.  Even  those  with  whiskers  that 
lack  undulatory  features,  but  are  flattened  in  profile,  may  be  able  to  employ  the  same 
mechanism  of  detection,  since  the  data  in  this  thesis  shows  the  same  major  trends  in 
response  of  the  harbor  seal  whisker  geometry  compared  to  a  smooth  elliptical  cylinder. 

5.4  Future  Work 

Certainly  this  is  merely  the  beginning  of  exploring  the  space  where  flow-induced 
vibrations  and  hydrodynamic  perception  meet.  The  following  are  a  few  areas  of 
immediate  interest  for  future  studies. 

5.4.1  Optimizing  the  whisker  geometry 

The  two  whisker  geometries  tested  in  this  thesis  were  either  constructed  directly  from 
provided  measurements  [58]  or  made  with  a  single  modification  of  a  single  parameter. 
A  thorough  optimization  of  the  whisker  geometry  could  lead  to  further  improved 
VI V  reduction  properties.  Appendix  A  provides  images  of  real  whiskers,  which  are 
intended  to  aid  such  a  study  for  any  interested  parties. 

It  is  interesting  to  note  that  the  two  geometries  tested  in  this  thesis  have  undula¬ 
tion  wavelengths  which  align  fairly  well  with  the  minima  in  vibration  response  seen 
in  wavy  cylinders,  as  described  in  Section  1.2.3.  Perhaps  some  harbor  seal  vibrissae 
have  evolved  to  a  nondimensional  wavelength  of  A / Dm  =  2  and  more  have  evolved 
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to  the  better  value  of  X/Dm  =  6  -  the  two  local  minima  in  responses  reported  in 
[80].  Also,  engineering  solutions  appear  to  have  converged  with  biology,  as  the  opti¬ 
mal  wavelength  to  diameter  ratio  for  strakes  is  about  5.  Optimization  of  the  whisker 
geometry  would  help  determine  what  wavelength  provides  the  global  minimum  in 
vibration  response. 

An  attempt  to  design  a  “whiskerlike”  geometry  with  low  VIV  for  all  angles  of 
oncoming  flow  is  described  in  Appendix  B.  Further  attempts  can  be  conducted  to 
achieve  a  reduction  greater  than  the  ~  20%  found  there.  The  elliptical  cross  section 
is  seen,  however,  to  be  the  dominant  feature  responsible  for  VIV  reduction,  so  it  is 
not  believed  that  further  optimization  would  achieve  significantly  better  results. 

One  key  parameter  that  differs  between  the  morphology  of  real  whiskers  and  the 
models  tested  in  this  thesis  is  that  of  curvature.  Real  whiskers  are  curved  along  the 
span,  and  it  would  be  good  to  determine  how  much  the  VIV  response  of  the  harbor 
seal  whisker  changes  with  the  addition  of  curvature.  Footage  of  swimming  harbor 
seals  by  [95]  led  those  authors  to  observe  that  the  major  curvature  of  the  hair  shaft 
was  generally  in  the  downstream  (caudal)  direction.  In  that  case,  the  curvature  would 
likely  not  dramatically  increase  the  VIV  response. 

5.4.2  Determining  the  limits  of  detection 

Distance  and  time  limits 

Chapter  4  showed  the  whisker’s  ability  to  detect  the  presence  of  a  vortex  wake  at  least 
160  whisker  diameters  away.  This  was  done  with  a  relatively  small  wake-generating 
body,  so  it  is  expected  that  a  larger  object  (which  would  more  closely  match  the  size 
ratio  between  a  real  whisker  and  a  fish),  would  leave  a  stronger  wake  and  therefore 
enable  detection  from  an  even  further  distance.  The  time  delays  reported  in  the 
behavioral  experiments  are  ~  30sec  long,  which  is  equivalent  to  tens  of  meters.  It  is 
of  interest  to  quantify  the  limit  of  detection,  in  terms  of  distance/time. 

An  order-of-magnitude  argument  can  be  made  for  this  limit,  utilizing  the  known 
sensitivity  of  the  whiskers  and  the  known  evolution  of  velocities  in  the  wake  of  bluff 
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bodies.  Based  on  a  dipole  experiment,  harbor  seals  can  detect  flow  speeds  as  low 
as  0.25 mm/s  [24],  The  relationship  of  the  wake-averaged  defect  velocity  (Uq/U) 
behind  a  bluff  body  versus  downstream  distance  (x/d)  is  shown  in  [123].  Given  a 
fish  swimming  at  U  =  O(10_1)m/s,  the  seal  would  theoretically  be  able  to  detect 
Uq/U  =  0(10-1)/0(10-4)  =  0(10_3)m/s.  Based  on  the  data  compiled  in  [123],  this 
value  equates  to  O(103)  (using  a  conservative  estimate  of  the  trend  for  non-stratified 
wakes).  This  means  that  a  small  fish  of  size  0(10) cm  could  be  detected  from  0(100) 
meters  away. 

Other  types  of  wakes 

In  this  thesis,  a  von  Karman  street  was  employed  as  a  representative  wake.  An 
efficiently  swimming  fish,  however,  will  leave  a  reverse  von  Karman  street,  which 
has  vortices  rotating  in  the  opposite  direction  and  produces  thrust.  It  should  be 
determined  whether  or  not  the  whisker  response  is  the  same  in  that  wake.  Also,  to 
more  closely  represent  the  real  scenario,  the  upstream  disturbance  should  be  oriented 
perpendicular  to  the  whisker  length.  I  postulate  that  these  two  changes  will  not 
induce  dramatic  changes  in  the  response,  as  the  dominant  wake  features  will  be  the 
same. 

Additional  variations  in  the  wake,  such  as  bifurcations,  could  be  tested  for  de¬ 
tectability  by  the  whisker.  Fish  leave  a  range  of  hydrodynamic  trails  depending  on 
their  kinematics,  and  therefore  do  not  always  leave  a  perfect  Karman  street.  Such 
tests  may  help  uncover  the  limits  of  detection. 

5.4.3  Implementing  a  control  loop 

The  device  developed  here  is  able  to  passively  detect  the  presence  of  a  vortex  wake. 
For  robotics  applications,  the  next  technological  hurdle  is  to  conduct  vortex- wake 
tracking,  which  requires  actuation  and  control,  but  would  enable  localization  of  a 
moving  source.  The  first  step  could  be  to  design  a  basic  control  system  to  detect 
the  location  of  a  wake  coming  from  an  unknown  but  stationary  position.  Then  the 
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control  could  be  enhanced  to  detect  the  heading  of  the  wake  and  follow  it  in  real  time 
as  the  wake- generating  body  turns. 

5.4.4  Determining  the  function  of  the  whisker  array 

In  this  thesis,  I  show  that  a  single  whisker  is  sufficient  to  perform  wake  detection. 
To  carry  out  more  complex  tasks,  however,  multiple  whiskers  may  be  required  or 
beneficial.  Studies  of  the  interaction  between  the  whiskers  may  uncover  additional 
insight  into  how  multiple  sensory  inputs  can  be  used  to  conduct  such  tasks. 

The  shear  number  of  whiskers  (~  40  on  each  side  of  the  muzzle)  suggests  that  there 
exists  some  beneficial  function  to  having  an  array.  It  is  known  that  multiple  whiskers 
can  be  used  in  concert  to  aid  discrimination  of  size  and  shape  in  tactile  sensing,  but 
the  hydrodynamic  function  of  the  vibrissal  array  is  a  subject  to  be  explored. 

The  arrangement  of  mystacial  vibrissae  in  the  array  across  the  muzzle  is  known: 
There  are  several  horizontal  rows  parallel  to  the  mouth  line,  the  vibrissae  are  generally 
longer  and  thicker  in  the  posterior  and  ventral  directions,  and  the  lower  rows  tend 
to  have  more  vibrissae.  Muscles  in  the  muzzle  allow  the  seal  to  position  the  whiskers 
in  a  variety  of  arrangements.  The  position  of  individual  whiskers  during  real  prey 
pursuit  is  unknown.  If  the  videos  provided  in  [95],  which  show  seals  swimming  in 
lab  conditions,  depict  the  same  whisker  positions  during  wake  tracking,  the  following 
are  a  few  notable  features  and  corresponding  hypotheses  about  the  hydrodynamic 
implications  of  this  arrangement.  Screenshots  from  the  videos  are  provided  in  Figure 
5-4. 

•  The  “fan”  arrangement  of  the  array  places  the  whiskers  at  many  different  angles 
to  the  horizontal  plane.  This  is  likely  done  to  cover  a  wide  area  of  the  flow.  The 
row  or  individual  whisker  that  vibrates  the  most  can  indicate  which  direction 
the  wake  is  coming  from. 

•  It  appears  that  despite  being  protracted  into  a  “fan” ,  all  whiskers  remain  facing 
the  flow  in  their  streamlined  orientation  (a  ~  0°),  indicating  that  they  will  all 
experience  low  VIV. 
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(a)  (b) 


Figure  5-4:  (a)  Front  and  (b)  side  views  of  the  arrangement  of  whiskers  as  a  harbor 
seal  swims  forward.  Screenshots  taken  from  supplementary  video  in  [95]. 

•  Along  each  row,  the  whiskers  increase  in  length.  This  could  be  used  to  estimate 
the  size  of  a  flow  feature,  such  as  a  vortex  which  may  hit  a  short  whisker  but 
not  the  tip  of  a  long  one. 

•  Alternatively,  the  range  in  lengths  of  the  whiskers  could  exist  to  help  determine 
frequency  information.  As  whiskers  with  different  lengths  will  have  different 
natural  frequencies,  each  has  a  different  range  that  it  will  respond  best  to. 

•  The  curvature  of  the  whiskers  within  one  row  generally  increases  from  anterior 
to  posterior.  This  could  possibly  serve  to  reduce  wake  effects  generated  by  the 
upstream  whiskers  felt  on  the  downstream  ones. 

•  The  most  anterior  whiskers,  which  are  the  shortest,  may  not  play  a  major  role 
in  hydrodynamic  sensing,  as  they  do  not  reach  far  out  into  the  flow.  Perhaps 
they  are  mostly  used  for  fine  tactile  sensing.  As  seen  in  [19],  a  sea  lion  first  used 
the  longer  vibrissae  to  locate  an  object  and  then  the  shorter  ones  on  top  of  its 
snout  to  collect  finer  details  about  the  shape. 

More  information  about  the  whiskers’  response  under  real  prey  pursuit  conditions 
remains  to  be  uncovered.  A  fine- resolution  view  of  the  whiskers  encountering  various 
hydrodynamic  trails  will  greatly  aid  our  understanding  of  how  the  array  may  be 
employed  for  sensing. 
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Appendix  A 


Whisker  Image  Database 


This  appendix  compiles  several  pictures  of  real  harbor  seal  whiskers.  These  images 
are  intended  to  serve  as  a  tool  for  further  studies  of  the  whisker  shape,  which  may 
include  the  goal  of  optimization  for  VIV  suppression. 

Figures  A-l  and  A-2  are  pictures  of  five  different  whiskers  shed  by  harbor  seals 
at  the  New  England  Aquarium.  The  whiskers  were  mounted  vertically  in  a  water 
tank.  Images  were  taken  using  digital  holography  (DH)  with  the  help  of  Yi  Liu  in 
the  MIT  Optics  Lab.  Back-propagation  was  used  to  focus  the  images  at  the  correct 
plane.  Each  whisker  was  captured  in  multiple  sections  due  to  the  long  whisker  length 
relative  to  the  camera  field  of  view. 

Figure  A-3  is  an  image  from  [95],  which  shows  cross  sections  of  different  seal 
species’  whiskers  based  on  CT  scans.  All  species  show  flattening  towards  the  tip. 
The  sea  lion  whisker  is  essentially  an  elliptical  cylinder  the  whole  way  down  the  span 
and  the  latter  two  show  the  undulatory  features  along  the  span  by  alternating  from 
a  roughly  circular  to  elliptical  cross  section. 

Images  and  dimensions  are  also  provided  in  [58],  but  they  are  not  replicated  here 
since  there  may  be  an  error  in  the  wavelength  reported.  Finally,  some  dimensions  are 
given  in  [46,  45],  which  are  copied  over  in  Figures  A-4,A-5,  and  A-6. 
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Figure  A-l:  Digital  holographic  images  of  whiskers  from  a  harbor  seal  at  the  New 
England  Aquarium.  All  are  in  the  90°  configuration.  The  camera  field  of  view  and 
therefore  the  dimensions  of  each  DH  image  box  are  1.23  x  1.23cm. 
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Figure  A-2:  Digital  holographic  images  of  whiskers  from  a  harbor  seal  at  the  New 
England  Aquarium.  All  are  in  the  90°  configuration. The  camera  field  of  view  and 
therefore  the  dimensions  of  each  DH  image  box  are  1.23  x  1.23cm. 
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Figure  A-3:  From  CT  scans,  these  are  reconstructions  of  (A)  sea  lion  (B)  elephant  seal 
and  (C)  harbor  seal  vibrissae.  Cross  sections  are  enlarged  and  shown  at  6  different 
locations  along  the  span.  Figure  from  [95] 


P.  groenlandicus 

C.  cristata 

H.  grypus 

Number  of  Vibrissae 

33-41  (38  ±  4) 

24-39  (32  ±  11) 

25 

Overall  length  (mm) 

39-3  ±  17.3 

33-4  ±  16.8 

37.8  ±  13.1 

Total  bumps 

12  ±  5 

10  ±  6 

9  ±  3 

Peak-to-peak 

Proximal 

2.65  ±  0.42 

2.24  ±0.13 

3.50  ±0.40 

distance  (mm) 

Middle 

3.00  ±  0.42 

2.50  ±  0.48 

3.90  ±  0.37 

Distal 

2.70  ±  0.55 

2.26  ±  0.24 

3.80  ±  0.57 

Tip 

2.50  ±  0.50 

2.00  ±  0.34 

3.90  ±  0.49 

Crest  width  (mm) 

Proximal 

0.63  ±  0.07 

0.58  ±  0.08 

0.94  ±0.17 

Middle 

0.61  ±  0.04 

0.57  ±  0.09 

0.97  ±0.19 

Distal 

0.45  ±  0.08 

0.45  ±  0.06 

0.89  ±0.18 

Tip 

0.34  ±  0.06 

0.31  ±0.06 

0.91  ±0.15 

Trough  width  (mm) 

Proximal 

0.47  ±  0.06 

0.46  ±  0.09 

0.80  ±0.16 

Middle 

0.43  ±  0.05 

0.42  ±  0.07 

0.76  ±0.15 

Distal 

0.33  ±  0.07 

0.36  ±  0.04 

0.72  ±0.14 

Tip 

0.25  ±  0.06 

0.25  ±  0.06 

0.74  ±0.14 

Figure  A-4:  Table  showing  measured  geometric  parameters  of  the  vibrissae  of  these 
three  species:  harp,  hooded,  and  grey  from  left  to  right.  Figure  from  [46]. 
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Crest  Width  (mm) 


Top  Peak-to-Peak 
Distance  (mm) 

Bottom  Peak-to- 
Peak  Distance  (mm) 

Crest  Width  (mm) 

Trough  Width 
(mm) 

Crest  Width/ 
Trough  Width 

Total  Length  (mm) 

#  of  Beads/cm 

Harp 

3.97  ±0.70 

3.94±0.69 

0.88  ±0.1 4 

0.69±0.13 

1.28 

80.6±10.8ab 

2.2±0.4bc 

Harbor 

3.27±0.39 

3.26  ±0.40 

0.92±0.13 

0.73  ±0.1 2 

1.26 

86.5±10.5ab 

2.3±0.4abc 

Ringed 

3.56  ±0.73 

3.53  ±0.72 

0.70  ±0.21 

0.49±0.21 

1.44 

79.6±10.9ab 

2.5±0.5ab 

Spotted 

4.01  ±0.63 

3.99  ±0.63 

1.04  ±0.1 5 

0.85  ±0.1 5 

1.22 

91.2±7.9a 

1.9±0.2C 

Gray 

3.43  ±0.38 

3.41  ±0.39 

0.76  ±0.1 3 

0.63  ±0.11 

1.21 

73.8±12.5b 

2.7±0.5a 

Bearded 

3.64  ±0.66 

3.62  ±0.65 

1.18±0.27 

1.18±0.27 

0.99 

77.2±9.5ab 

0 

Fur  Seals 

3.64  ±0.66 

3.62  ±0.65 

0.96±0.12 

0.96±0.12 

1.01 

83.0±20.3ab 

0 

Sea  Lions 

3.64±0.66 

3.62  ±0.65 

0.84  ±0.1 9 

0.84±0.19 

1.00 

90.7±11.3a 

0 

Mean  ±  SD  values  for  each  species  are  given.  Significant  differences  between  species  for  total  length  and  number  of  beads  per  cm  are  indicated  by  different  letters. 
Significant  differences  between  species  for  the  other  traditional  morphometric  measurements  are  given  in  Figure  3.  Species  with  different  letters  are  significantly 
different  from  one  another. 
doi:1 0.1 371  /journal.pone.0034481  .t002 


Figure  A-5:  Table  showing  measured  geometric  parameters  of  the  vibrissae  of  multiple 
species.  Figure  from  [45]. 


Figure  A-6:  (a)  Crest  width  and  (b)  trough  width  measured  on  vibrissae  of  different 
species.  The  diameter  decreases  towards  the  tip.  Figure  from  [46]. 
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Appendix  B 


Omnidirectional  Whisker-Inspired 
Geometries 


Can  a  “whiskerlike”  geometry  be  designed  to  reduce  VIV 
regardless  of  oncoming  flow  angle? 

B.l  Experiment  Motivation 

A  wide  range  of  vibrational  response  over  all  angles  of  attack  can  be  useful  for  a 
sensor,  but  it  is  a  poor  feature  for  a  VIV  reduction  mechanism,  particularly  in  the 
ocean  where  currents  may  come  from  multiple  directions.  Therefore,  it  is  desired  to 
know  if  a  “whiskerlike”  geometry  can  be  created  to  experience  low  VIV  over  a  wide 
range  of  flow  directions. 

B.2  Shape  Design 

Seven  geometries  were  designed,  based  on  various  features  of  the  original  whisker 
shape.  They  were  fabricated  into  plastic  via  stereolithography,  a  layer-by-layer  pro¬ 
cess.  Each  model  has  the  same  length  ratio  (L/d)  as  the  original  whisker  and  is  made 
from  the  same  material.  Their  mass  ratios  range  from  m*  =  0.9  to  1.8. 
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(a) 


(b) 


Figure  B-l:  Screenshots  of  Model  1  and  Model  2,  which  were  designed  keeping  all 
dimensions  the  same  as  the  original  whisker  except  that  one  dimension  of  each  cross 
section  was  expanded  or  shrunk  to  make  it  circular. 


B.2.1  Direct  Conversion  to  Circular  cross  sections 

The  first  two  (Model  1  and  Model  2)  take  the  cross-section  dimensions  provided  in  [58] 
and  change  one  of  the  two  diameters  such  that  each  cross  section  becomes  circular. 
The  same  wavelength  as  the  original  whisker  is  maintained.  See  Figure  B-la  and 
B-lb. 


B.2.2  “Spiral-like”  Shapes 

The  next  four  (Models  3,  4,  5,  and  6)  are  created  via  parametric  equations.  The 
equation  for  a  standard  spiral  (Equation  B.l)  is  repeated  multiple  times,  each  time 
with  slight  offset  in  the  spanwise  (z)  direction  and  in  amplitude.  Put  together,  these 
“strips”  of  helices  form  a  spiral  body  with  undulating  amplitude.  The  base  helix 
has  amplitude  (a)  of  0.475mm,  the  smallest  radius  of  the  whisker.  As  shown  in 
Equation  B.2,  over  each  one  of  n  iterations,  the  amplitude  is  increased  by  a  constant 
amount  (A a)  up  to  the  maximum  radius  of  0.595mm,  following  the  dimensions  of 
the  whisker  given  in  [58].  The  2  offset  is  also  increased  by  a  constant  amount  (half 
wavelength/number  of  iterations)  in  each  iteration.  The  same  process  is  then  repeated 
with  decreasing  amplitude  to  create  strips  of  helices  back  down  to  the  base  spiral 
radius  but  with  increasing  offset  in  the  2  direction  (Equation  B.3).  This  strip  method 
is  depicted  visually  in  Figure  B-2. 
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A  min 


A 


max 


Figure  B-2:  Spiral  strips  combined  together.  The  blue  and  red  lines  are  the  initial  and 
final  spirals.  Wavelength  (A)  and  the  two  amplitudes  (Amin  and  Amax )  are  labeled. 


x  =  acos(u>t ) 
y  =  a  sm(ut) 

z  =  t  (B.l) 

a  =  a  +  Aa 

z  =  z  T  — 
n 

a  =  a  —  A  a 

z  =  z  +  —  (B.3) 

n 

Model  3  follows  the  same  wavelength  as  that  of  the  previous  two  models  and  the 
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Figure  B-3:  All  the  models  created  via  parametric  equations.  From  the  left:  Model 
5,  Model  4,  Model  3,  Model  6. 

original  whisker  dimension.  Model  4  uses  half  of  the  wavelength.  Model  5  was  created 
using  the  same  parametric  equations  as  above,  but  now  with  a  more  extreme  change  in 
diameter.  The  overall  minimum  and  maximum  radii  of  the  actual  whisker  (0.24mm 
and  0.595mm)  were  used  as  the  base  and  top  radii  of  the  spiral  body,  in  order  to 
create  a  model  with  more  pronounced  undulations.  Model  6  was  created  using  the 
original  whisker  dimensions,  but  with  a  sinusoidal  change  in  amplitude  rather  than  a 
linear  one. 

The  coordinates  generated  from  these  equations  are  then  meshed  together  via 
MeshPlot,  imported  into  SolidWorks,  scaled  up  to  match  the  mean  cross-stream  di¬ 
mension  of  the  previous  models  (~  1  in),  and  fabricated.  Figure  B-3  shows  all  the 
models  created  via  the  parametric  equations. 

B.2.3  Twisted  Ellipse 

The  final  model  generated  is  that  of  a  twisted  elliptical  cylinder.  This  geometry  was 
created  by  taking  the  average  elliptical  cross  section  of  the  whisker  geometry  and 
twisting  it  along  the  spanwise  axis.  The  amount  of  twist  applied  was  chosen  such 


SolidWorks  Educational  Edition. 
For  Instructional  Use  Only. 


Figure  B-4:  The  twisted  elliptical  cylinder,  as  generated  in  Solidworks 

that  the  wavelength  of  undulations  matched  that  of  the  whisker.  This  value  is  given 
by  (f)  =  2.95 z/d,  where  z  is  the  span  wise  position.  A  picture  of  the  model  is  shown 
in  Figure  B-4. 

Work  by  [76]  numerically  tests  the  forces  felt  on  a  stationary  body  with  a  similar 
shape.  Their  geometry  was  also  generated  by  twisting  an  elliptical  cylinder,  but  it 
has  slightly  different  parameters:  A/davg  =  0.1,  X/davg  =  tt/3  =  1.05,  <fi  =  2 z/d, 
Re  =  3,  900,  where  A  is  the  amplitude  of  undulation.  In  my  case,  those  parameters 
are  A1/davg  =  0.67,  A2/davg  =  0.33,  A /davg  =  1.14,  <fi  =  2.95 z/d,  and  Re  =  10,800. 
These  geometrical  parameters  do  not  differ  dramatically,  so  the  results  from  the  two 
models  can  be  broadly  compared. 


B.3  Forced  Vibrations:  Force  Measurements 

B.3.1  Apparatus 

Forces  were  measured  as  the  models  underwent  a  range  of  imposed  motions.  The 
hydrodynamic  forces  are  used  to  calculate  nondimensional  coefficients  that  describe 
the  real  vibrational  response.  The  same  linear  motor  and  force  sensor  apparatus  used 
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CLv  CLv 


Figure  B-5:  CijV  contours  of  Models  1  and  2. 


in  Chapter  2  was  employed  for  the  present  experiments. 


B.4  Results 
B.4.1  C,v 

The  zero  C/v  contour  is  used  to  describe  the  VIV  response  of  the  geometry  under  zero 
damping  conditions.  The  peak  therefore  indicates  the  peak  expected  VIV  amplitude. 

Models  1  and  2,  which  are  shown  in  Figure  B-5  have  a  peak  amplitude  of  ~  0.65. 
Both  of  those  first  two  models  thus  present  a  reduction  in  VIV  amplitude  of  about 
20%  compared  to  that  of  the  smooth  circular  cylinder  shown  in  Figure  2-8. 

Results  from  the  models  created  via  parametric  equations  are  shown  in  Figure 
B-6.  Similar  to  above,  the  peaks  fall  between  A/d**  =  0.65  to  0.8. 

Finally,  the  C%,  contour  for  Model  7,  the  twisted  ellipse,  is  shown  in  Figure  B-7. 
The  peak  amplitude  is  found  to  be  A/d*  ~  0.6. 

The  zero  contour  peaks  for  all  of  the  models  tested  in  this  study  are  similar  to  that 
of  a  circular  cylinder,  based  on  literature  [49]  and  my  results  shown  in  Figure  2-8. 
The  maximum  reduction  seen  is  in  the  twisted  ellipse,  with  a  reduction  of  ~  15%. 
None  of  them  achieve  the  significant  VIV  reduction  seen  in  the  whisker  or  smooth 
elliptical  cylinder.  Also,  these  results  show  a  weaker  reduction  compared  to  that 
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Amplitude  Ratio  (Y/d*)  Amplitude  Ratio  (Y/d*) 


CLv 


Figure  B-6:  C^v  contours  of  Models  3,  4,  5,  6. 


CLv:0  deg 


Reduced  Frequency  (fd*/U) 


Figure  B-7:  Ci}V  contours  of  Model  7. 
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in  [76],  which  reports  a  reduction  in  mean  drag  up  to  14%  and  reduction  in  root 
mean  square  ( rms )  of  fluctuating  lift  up  to  74%.  This  large  difference  highlights  the 
importance  of  including  dynamic  motion  in  testing,  as  noted  in  [99].  Allowing  body 
motions  to  develop  can  strongly  affect  the  resulting  drag  and  lift  forces  experienced. 
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Appendix  C 


Effect  of  Flexibility  on  VIV 


How  does  flexibility  affect  the  whiskers ;  vibrations  ? 

C.l  Motivation 

Real  pinniped  vibrissae  are  somewhat  flexible.  All  of  the  experiments  previously 
described  in  this  thesis  were  conducted  on  rigid  models.  It  is  desired  to  know  how 
flexibility  may  change  the  vibration  properties  of  the  harbor  seal  whisker  geometry. 
Flexible  circular  cylinders  are  known  to  experience  different  VIV  responses  than  rigid 
ones,  including  features  such  as  vibrating  along  the  span. 

C.2  Experiment  Setup 

C.2.1  Model  Fabrication 

The  following  is  an  overview  of  the  experiment  design  and  fabrication.  The  full  details 
can  be  found  in  Christopher  Gerber’s  Bachelor’s  Thesis  [44], 

The  experimental  setup  was  designed  around  the  following  parameters:  compliant 
material,  adjustable  natural  frequency,  and  minimal  deformation  as  a  result  of  that 
adjustment.  The  chosen  solution  was  to  cast  a  long,  flexible  rubber  model  with  two 
separate  Kevlar  strings  running  down  the  center.  The  high  strength  of  the  Kevlar 
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(a)  (b) 


Figure  C-l:  (a)  One  half  of  the  whisker  mold,  which  was  made  from  six  blocks  of 
machinable  wax  bolted  together.  A  CNC  tool  path  for  the  whisker  geometry  was 
created  in  MasterCAM  from  Solidworks  and  then  routed  out.  (b)  The  mount  to 
which  the  whisker  was  attached.  The  model  was  strung  in  the  empty  space  pictured 
at  the  bottom.  The  two  struts  have  an  airfoil  profile  to  minimize  drag.  The  tube  on 
top  attaches  directly  to  the  front  of  the  Tow  Tank  carriage. 


allowed  the  model  to  be  held  taut  in  a  base  position,  and  turnbuckles  were  used  to 
tension  it  to  a  desired  level.  The  flexibility  of  the  rubber  allows  vibrations  to  develop 
along  the  body. 

The  whisker  geometry  dimensions  given  in  [58]  were  used  directly,  and  scaled  up 
lOOx,  giving  an  average  crossflow  diameter  of  d  =  5.3 cm.  A  mold  with  the  negative 
of  this  geometry  was  made  out  of  machinable  wax  using  a  CNC  mill  (Figure  C-la).  A 
two-part  rubber  (Smooth-On  PMC  724)  was  cast  into  the  mold  to  create  the  whisker 
model. 

The  length  of  the  model  (6  ft)  was  chosen  such  that  it  would  fit  an  existing  mount 
(Figure  C-lb)  that  was  used  for  previous  studies  in  the  MIT  Towing  Tank.  The 
diameter  was  chosen  such  that  the  aspect  ratio  would  be  at  least  20  to  allow  for 
potential  motion  along  the  span.  Also,  it  was  chosen  such  that  embedded  sensors 
(0.5m  cubes)  would  be  small  compared  to  the  model,  to  add  minimal  rigidity  to  the 
otherwise  flexible  structure. 

Accelerometers  (Kistler  8688A5,  +/  —  5 g)  were  placed  inside  the  model  and  used 
to  measure  model  vibration  without  disturbing  the  flow,  as  a  string  potentiometer 
would,  for  example.  Three  of  them  were  embedded  in  the  center  of  the  model  during 
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(a)  (b) 


Figure  C-2:  (Left)  Accelerometer  held  in  place  via  strings  and  a  thin  piece  of  plastic. 
(Right)  Three  accelerometers  and  two  Kevlar  strings  inside  of  the  mold  just  before 
casting. 

casting.  They  were  placed  evenly  along  the  span  with  one  in  the  middle.  They  were 
held  in  place  via  strings  strewn  across  the  diameter  of  the  mold  and  a  thin  piece  of 
plastic  in  the  direction  of  the  other  diameter  (Figure  C-2a).  Figure  C-2b  shows  all 
the  items  (three  accelerometers  and  two  Kevlar  strings)  inside  of  the  mold  before  the 
rubber  was  poured  in. 

C.2.2  Mounting  Scheme 

After  experimentation  with  various  boundary  conditions,  it  was  determined  that 
“clamped-clamped”  was  the  best  option.  Because  of  the  weight  of  the  model,  the 
model  tended  to  sag  along  the  span.  The  strings  alone  were  insufficient  to  keep  the 
body  level  along  its  span,  so  fixed  boundary  conditions  helped  reduce  this  effect.  The 
clamps  were  used  to  support  a  majority  of  the  weight  and  the  strings  provided  the 
additional  tension  needed  to  pull  the  body  to  a  horizontal  position. 

A  schematic  of  the  setup  is  shown  in  Figure  C-3a,  and  Figure  C-3b  shows  the 
entire  system  mounted  in  the  Towing  Tank. 

C.2.3  Modeling 

The  expected  natural  frequencies  of  the  structure  can  be  calculated  using  beam  theory 
with  the  appropriate  boundary  conditions  and  assuming  the  body  is  in  a  vacuum. 
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Figure  C-3:  (a)  Schematic  of  the  setup.  The  whisker  is  held  with  clamps  on  the  end. 
Turnbuckles  are  used  to  apply  an  initial  tension  to  the  Kevlar  strings.  The  spring 
and  pulley  system  on  both  ends  allows  for  additional  stretching  as  vibrations  occur 
along  the  body,  (b)  Picture  of  the  assembled  system  mounted  in  the  Towing  Tank. 


The  structure  can  be  modeled  as  a  beam  and  as  a  string.  The  natural  frequency  is 
calculated  for  both  cases  and  combined  to  determine  the  overall  natural  frequency,  as 
shown  in  equations  C.l.  This  system  is  seen  to  be  tension-dominated.  The  calculation 
is  repeated  for  the  first  three  modes  of  vibration  (n  —  1,  2,  3).  A  pluck  test  was  done 
on  the  structure  while  in  air  and  showed  the  first  mode  to  vibrate  at  fn:air  =  2.54Hz. 

fbeam  =  =  [°’2>  °-9,  2.0\Hz  (C.la) 

f string  =  n\\j =  [2-7,  5.5,  8.2 }H Z  (C.lb) 

ftotai  =  yj faring  +  fL*n  =  [2-73,  5.53, 8.43] (C.lc) 


The  expected  natural  frequency  in  water  can  be  calculated  using  the  expected 
value  in  air  and  incorporating  the  appropriate  added  mass  and  mass  ratio  values,  as 
in  equation  C.2.  In  this  case,  the  body  is  approximated  as  a  uniform  ellipse  with  a  2:1 
diameter  ratio,  giving  it  an  added  mass  coefficient  Cm  of  4.  The  mass  ratio  m*  =  1.38 
was  determined  from  the  Smooth-On  PMC-724  material’s  data  sheet  value  for  specific 
gravity.  A  crossflow  direction  pluck  test  on  the  model  in  water  showed  an  oscillation 
frequency  of  fn, water  =  1.416 Hz,  which  matches  the  expected  value  very  well. 


176 


Go  Pro  Hero  3 


Mourn 


Spring 


Pulley 


Tumboclde 


Kevlar  String 


Marker 


Whisker 


Figure  C-4:  Schematic  of  the  setup,  which  now  includes  a  GoPro  camera  to  capture 
motion  of  markers  on  top  of  the  whisker  at  center  span. 


(C.2) 


C.2.4  Motion  Tracking 

Matthew  Hildner  developed  a  second  means  of  measuring  vibrational  response.  Opti¬ 
cal  measurements  of  the  whisker’s  motion  were  collected  through  use  of  an  underwater 
camera  and  markers  on  the  whisker  surface.  This  method  allows  for  quantification  of 
motion  in  all  three  degrees  of  freedom. 

A  schematic  of  the  modified  setup  is  shown  in  Figure  C-4.  A  GoPro  Hero  3  is 
mounted  ~  12m  above  the  whisker  to  minimize  flow  disturbance  seen  by  the  whisker. 
It  is  oriented  with  a  30°  angle  to  the  markers  so  as  to  capture  all  three  degrees  of 
motion  simultaneously. 

Camera  Calibration 

Calibration  of  the  camera  from  pixels  to  dimensional  length  was  carried  out.  Because 
of  the  orientation  angle  of  the  camera,  the  distance  in  the  image  of  the  crossflow 
direction  (depicted  with  the  marked  pole  in  Figure  C-5)  changes  exponentially. 

Video  Processing 

An  object  tracking  algorithm  was  developed  in  MATLAB  based  on  a  few  existing 
toolboxes.  First,  the  GoPro  Cineform  software  was  used  to  adjust  the  color  of  the 
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Figure  C-5:  A  calibration  image.  In  this  case  a  pole  marked  with  1  in  increments  is 
used  to  determine  the  pixel  to  crossflow  distance  conversion. 


(a)  (b) 

Figure  C-6:  A  few  of  the  color  adjustments  done  to  better  highlight  the  markers. 


images.  This  helped  isolate  the  green  thumb  tack  markers  from  the  white  whisker 
body  and  blue  tank.  All  objects  besides  the  markers  were  converted  to  a  red  (Figure 
C-6a).  Then  a  series  of  transformations  were  done  in  MATLAB,  leaving  a  clear 
contrast  of  white  markers  against  a  black  background  (Figure  C-6b). 

Next,  MATLAB’s  Blob  Analyzer  was  used  to  identify  the  groupings  of  white 
pixels  and  track  their  location  over  time.  A  Kalman  filter  predicts  the  location  of  the 
markers  for  the  next  frame,  which  helps  reduce  incorrect  tracking. 


C.3  Results 

C.3.1  Accelerometer 

Data  from  the  middle  accelerometer  in  the  crossflow  direction  was  collected.  This 
signal  was  integrated  two  times,  detrended,  and  filtered  with  a  low-pass  Butterworth 
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Figure  C-7:  From  top  to  bottom:  crossflow  acceleration,  velocity,  and  displacement 
of  the  whisker  as  measured  from  the  middle  accelerometer  in  the  U*  =  4  case. 

filter  to  estimate  the  displacement.  An  example  plot  of  the  progression  from  acceler¬ 
ation  to  velocity  to  displacement  is  shown  in  Figure  C-7. 

C.3.2  Motion  Tracking  of  Model 

The  markers  were  tracked  over  the  course  of  various  towing  tank  runs,  their  posi¬ 
tions  were  converted  to  dimensional  lengths,  and  their  amplitudes  were  calculated 
in  terms  of  root- mean-square  ( rms ).  The  results  are  compiled  in  Figure  C-8a  along 
with  the  results  from  the  accelerometer.  The  dominant  frequencies  were  measured 
using  both  techniques  by  selecting  the  peak  response  in  the  FFT  of  the  signal  and 
nondimensionalizing  by  the  crossflow  natural  frequency  in  air  (Figure  C-8b). 

The  two  measurement  techniques  line  up  reasonably  well  with  each  other.  Both 
show  a  higher  peak  value  than  those  measured  in  the  rigid  whisker  experiments.  A 
high  response  is  seen  to  start  to  develop  at  high  velocities.  This  is,  however,  not 
believed  to  necessarily  indicate  galloping  because  the  model  was  observed  to  undergo 
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Figure  C-8:  The  nondimensional  crossflow  (a)  amplitude  and  (b)  frequency  measured 
both  by  the  accelerometers  and  by  the  motion  tracking  software. 


Figure  C-9:  (a)  Mean  and  (b)  rms  rotation  of  the  model  across  different  towing 
velocities.  Both  increase  dramatically  at  high  speeds. 


significant  mean  rotation  along  the  spanwise  axis  at  high  speeds.  The  model  could 
no  longer  be  considered  as  being  tested  at  a  =  0°. 

The  rotation  was  quantified  through  the  motion  tracking  software.  Both  mean 
rotation  and  rms  rotation  are  shown  in  Figure  C-9.  As  towing  speed  increases,  the 
mean  rotation  increases.  The  rms  rotation  is  non-zero  in  all  cases,  and  reaches  a 
peak  at  U  =  0.8m/ s. 
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C.4  Discussion  and  Future  Work 


Both  in-line  and  crossflow  vibrations  reach  higher  values  in  this  flexible  whisker  setup 
than  in  the  rigid  whisker  experiments.  Also,  a  significant  torsional  excitation  was 
measured,  especially  at  certain  speeds.  The  excitation  of  this  mode  may  be  driven 
by  the  elliptical  cross  section  of  the  geometry  [137] . 

It  is  possible  that  torsional  vibration  is  responsible  for  inducing  the  additional 
vibrations  seen  in  the  crossflow  and  in-line  directions.  To  investigate  this,  we  are 
modifying  the  freely  vibrating  rigid-whisker  system  to  include  a  torsional  degree  of 
freedom. 
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